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Sorption of Nitrogen and Water Vapor on Textile Fibers’ 
By John W. Rowen and R. L. Blaine 


Measurements were made of the adsorption of nitrogen and water vapor on six purified 


textile fibers and titanium dioxide. 


It was found that all of the fibers had a relatively low 


capacity for the adsorption of nitrogen as compared with the capacity for the adsorption of 


water vapor. 


0.98 square meter per gram for viscose rayon. 


The surface area values ranged from 0.31 square meter per gram for nylon to 


The values of the free surface energies of 


adsorption as calculated by the Gibbs adsorption equation were the same for wool, cotton, 


silk, and rayon fibers but differed for the two synthetic polymers, nylon and acetate rayon 


I. Introduction 


There is evidence that surface characteristics of 
textile fibers have some effect on the properties of 
the finished-textile fabrics. Examples of such 
evidence of surface properties are apparent in the 
phenomenon of water-repellency [1, 2],? in the 
loss of light-reflectance [3] and in the difference in 
moisture-sorption capacities of fabries [4, 5] after 
various treatments. Although the importance of 
surface characteristics and properties has long been 
recognized in the technology of colloids and other 
sciences, the relation of these surface properties to 
the properties of fabries is not well understood. 
10, 11, 12] have in- 


ierpreted the adsorption of water by some textile 


Several workers [6, 7, 8, 9, 


materials (such as cellulose, wool, nylon, ete.) as 
a phenomenon in multimolecular [13] adsorption. 
If one assumes that the adsorption is a multilayer 
phenomenon, one may calculate that the surface 
area accessible to the first layer is in the range 
from 100 to 200 m?/g. 
several orders of magnitude greater than the 


This range of values is 


value (0.60 m?/g) caleulated by Emmett and 
DeWitt [14] from nitrogen adsorption measure- 
On the other hand, 
this range of values is not appreciably greater 
than the value (53.0 m?/g) reported by Purves [9] 
et al., working with nitrogen on a swollen cotton 


ments on paper cellulose. 


cellulose. 

ln order to obtain reliable surface area measure- 
ments of textile 
measure the surface area under conditions en- 


materials it is necessary to 





sper will also appear in the December 1947 issue of the J. Ind. and 


in brackets indicate the literature references at the end of this 
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suring true physical adsorption. Brunauer, Em- 
mett, and Teller [13] and Harkins and Jura [15] 
have shown that the adsorption isotherms of 
195°C (the boiling point of nitrogen, 
physically 


nitrogen, at 
at which temperature nitrogen is 
adsorbed on solids) may be used to calculate the 
true surface areas of a great many crystalline and 
amorphous materials. 

Knowledge of the molecular area or surface 
area [16] is believed necessary in order to evaluate 
such properties as the free surface energy [17] of 
adsorption and the work of adhesion [18] be- 
It therefore 
seemed highly desirable that the surface areas of 
cotton, wool, 


tween the liquid and the solid phase. 


the more common textile fibers 
viscose, silk, nylon, and cellulose acetate—be 
measured. It also seemed desirable to obtain the 
water adsorption isotherms of the same samples 
used in the nitrogen adsorption experiments. The 
work presented in the paper represents a portion 
of the preliminary work of a larger program 
sponsored and supported by the Office of the 
Quartermaster General, War Department.’ It is 
intended to provide comparative data on the 
adsorption of the two vapors on six different 


purified textile materials. 


II. Materials and Their Preparation 


Six purified textile fibers and asample of titanium 
dioxide (anatase) were used in this study. The 
four fibers, cotton, wool, silk, and rayon are from 
the samples studied by Wiegerink [19]. The 
titanium dioxide was employed as a_ reference 
standard. The cotton was purified by extracting 


* Permission to publish this paper was granted by the Office of the Quarter- 
master General 
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with aleohol for 8 hours and with ether for 8 
hours, followed by four washings in distilled 
water at 50° to 60°C. The wool was also ex- 


tracted with ether and alcohol until the extracting 
liquid was free of residue traces. The wool was 
50 
The raw silk was degummed in soap solution and 
then extracted with aleohol and ether. It 


then thoroughly washed in distilled water. 


then washed at to 60°C in distilled water. 
was 
The 
viscose rayon was washed four times in a 1-percent 
50° to 60°C and 
then washed four times with distilled water at the 


water solution of ammonia at 


same temperature. 

The sample of nylon used in this study was 40 
denier, 13 filament varn, having one-half turn per 
inch of Z twist. 
KE. I. 


for 6 hours with ethyl ether. 


The sample was supplied by the 
duPont de Nemours Co. It was extracted 
This length of time 
proved sufficient for the complete removal of the 
0.4 percent of oil base lubricant applied during 
spinning. The acetate rayon was 150 denier, 46 
filament yarn having 2) turns per inch and was 
by the American Co. The 
sample was washed four times in distilled water 
at 50° and 60°C. 


supplied Viscose 


III. Test Methods 


The sorption of water by the textiles was de- 


termined gravimetrically using an analytical 
balance in conjunction with the conditioning 


Worthing- 
In this apparatus the partial pressure 


apparatus deseribed by Carson and 
ton [20]. 
of the water is maintained by means of saturated 
salt solutions, and the samples may be weighed 
without removal from the conditioned atmosphere. 
The moisture vapor pressures were obtained from 
wet and dry bulb temperatures and also by means 
The dry 

samples of yarn 
them 


of a Dunmore [21] electric hygrometer. 
weights of the 
by 


5-gram were 


obtained placing over phosphorous 


pentoxide until constant weight was obtained. 
This constant weight was attained in about 5 days. 
It was observed that no loss of weight occurred 
from the fifth to the thirtieth day. After this 
drying period, the fibers were in equilibrium with 
an atmosphere whose relative water vapor partial 
pressure was less than 0.01. After the dry equi- 
librium point had been attained, the moisture 
content of the sample was determined at a series 
The 


of increasing relative water vapor pressures. 
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rate of change of weight was plotted against tin, 
at each partia! pressure, and the equilibrium \ aly 
usually attained within 24 to 48 hours, was 
indicated by the rate curve. 

The adsorption of nitrogen was determined }) 
a method similar to the one used by Brunauer and 
Emmett [22]. The saturation vapor pressure o 
the nitrogen was determined by use of a vapo, 
the 
employed for the adsorption experiment. 


1. 
The 


made on 20- to 50-gram samples of the fibers 


pressure thermometer using purified — gas 


low temperature adsorption tests we 
The sample bulb was packed as tightly as possibl, 
for each determination to reduce the dead spa 
correction. For the experiments reported in this 
article, the textile fibers were evacuated for 16 1 
20 hours at room temperature (about 25° C 
mm of Hg. 
inary tests were made in which the temperatu 


toa 


pressure of about 10 Other prelim- 
of evacuation was higher and one experiment was 
the 

For the low temperature adsorption of 


made of cotton-fiber surface available to 
oxygen. 
TiO, the temperature of evacuation was 200° ( 
Helium was used in evaluating the dead space in 


the sample bulb. 


IV. Results and Discussion 


1. Adsorption Measurements 


As shown in figures 1, 2, and 3, the adsorptior 
isotherms of both water vapor and nitrogen wer 
found to be sigmoid-shaped, corresponding to th 
type II of the five types classified by Brunaue: 
|233]. 
istic of physical adsorption on rigid solids with 


This type of adsorption curve is character- 


pore diameters larger than the molecular diameters 
of the adsorbed molecules [24]. However, this 
type of adsorption curve also occurs with non- 
rigid solids in which the adsorption sites are no! 
restricted to the surface. 

It may be noted from figure 1 that the mil- 
moles of water vapor and nitrogen adsorbed by 
the Tid, differ only by a small percentage [25 
This is true despite the fact that the nitrogen was 
adsorbed after thorough evacuation, whereas water 
vapor was adsorbed from mixtures in air afte! 
On the other hand 
comparing the adsorption isotherms of the textile 
fibers in figures 2 and 3, and the data in table ! 
reveals that the order of magnitude of adsorption 


drying the sample over P,Os. 


capacity of nitrogen and water vapor are not the 
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sam All the fibers had relatively low capacity 





MILLIGRAMS OF NITROGEN VAPOR PER GRAM OF SOLID 


MILLIMOLES OF NITROGEN VAPOR PER GRAM OF SOLID x 107 


1.0 
for the adsorption of nitrogen as compared with 
the relatively high capacity for adsorption of water 








vapor. 
“ik . , “4 . 1 , Oo ol 1] 03 04 05 o6 7 
The wool and viscose rayon fibers had slightly 7 ’ - 
‘ iy: ™.. RELATIVE NITROGEN VAPOR PRESSURE P/P, 
greater nitrogen adsorption capacities than silk 
. ae Ficure 2 Adsorption isotherms nitr x 
and cotton, and the adsorption capacities of the scums 3.—Adsorption seothermes of nitrogen on siz textile 
fibe rs at 199° C. 


acetate and nylon fibers were the lowest. In the 














Th ¢ . . , Wool; X, viscose; silk , cotton; , acetate , nylon, 
0 adsorption of water vapor, the wool and viscose 
rayon fibers also had greater adsorption capacity 
it 
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TABLE 2.—-Surface areas of tertile fibers 


! 


than silk and cotton, whereas the acetate and 
nylon also had the lowest capacities. 
certain lack of 


shifting of relative position. 


There was, 
and 
Sufficient study has 
not been made to indicate the significance of these 
differences. 


however, a parallelism Area available 
to nitrogen 
molecule at 
195° C (meg 


Wool 
Viscose rayon 


Silk 


abulate 


2. Calculated Quantities able th 


(a) Surface Areas 


The surface areas of the fibers available to the 
nitrogen vapor were calculated according to the 
BET equation [26], commonly used for 
calculating the surface areas of solid materials. 


now 


In evaluating the quantities of gas required fora 
monolayer on each of the textile the 
customary method of plotting P/V (Py—P) versus 
P/P. was employed (26). 
in figure 4. 
are presented in table 2. 


fibers, 


The plot is presented 
The calculated surface area values 
All of the surface areas 





Cotton 
Acetate rayon 
Nylon 


Titanium dioxide (TiO, 


Purves’ experimental value suggests. The sur. 
thus appears 
vary from at least 0.60 to 53.0 m*/g, depending 


upon its history. The for cellulose 0 


face area of cotton cellulose 


value 


0.60 m?*/g is only slightly larger than values 
(0.1 
observed with the microscope and reported by 


0.3 m*/g) obtained by estimation of areas 





0.2 0.3 


RELATIVE NITROGEN VAPOR PRESSURE P/P, 


Figure 4. 
, Wool; X, viscose; 0, silk; 


were less than | m’/g. The value obtained for 


cotton cellulose —0.72 m?*/g—compares favorably 
with the 0.60 m*/g value reported by Emmett 
and DeWitt [14]. However, it differs markedly 
from the value of 53.0 [9] reported by Purves for 
a sample of cellulose. It is 
pointed out that the surface area available to the 
nitrogen vapor may be greatly dependent upon the 
and treatment 


swollen cotton 


source as the comparison with 
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BET plots used for obtaining surface areas accessible 


, cotton; V, 


to nitrogen. 


wetate; A, nylon. 


Stamm and Millett [6] and Howell and Jackso 
27). It is likely, therefore, that the surface area 
as determined with nitrogen, represents a lower! 
limit. 

In view of the large differences between th 
the 
isotherm and the nitrogen isotherm, it was © 
the “‘water’’ 
This caleulation was carried ou! 


adsorption capacities observed in wat 


interest to calculate areas of thi 


same samples. 
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ough the authors were aware of the possi- 
hat the data could be interpreted from 











nother point of view. The plot of P/V(P)»—P) 
sus P/P, is shown in figure 5, and the areas 
hat may be calculated from them, by assuming 
|| adsorbing sites were on a free surface, are 
abulated in table 3. It may be noted in this 
able that the caleulated surface areas of the 
ding 
. 
lu 
rei 
fon »2 Oo3 04 O05 
RELATIVE WATER VAPOR PRESSURE P/P, 
oF nE 5 BET plots used Jor obtaining surface areas 


accessible to water vapor. 


xtile fibers were from 150 to 200 times as great 
vhen determined by water adsorption as_ they 


vere When determined by nitrogen adsorption. 


BLE 3 Calculated surface areas accessible to water vapor 


and nitrogen vapor 





Area (1) avail- | Area (II) avail- Ratio of 
Mat ible to water ible to nitro area I 
. molecules at gen molecul to area 
20” at — 195 
4 
yt 0. 06 21 
t 24 us 208 
140 76 is4 
OS 72 1M”) 
\ x8 is 1M 
40.0 31 14 
xid rio 7.0 7.9 0.9 


This apparent discrepancy may be explained in a 


mber of ways: 

The sites that adsorb water are not restricted 
lo free surface. 
2. If the water adsorbing sites are restricted to a 
surface, there may be an additional internal sur- 
lace specific to certain adsorbates as well as an 
al surface. 
Che internal surface within the fibrous struc- 
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ture exists only in the presence of a swelling agent 
such as water. 

4. The smaller diameter and the polarity of the 
water molecule enables it to penetrate into capil- 
laries not accessible to the nitrogen molecule. 

Stamm and Millett [6], after a careful considera- 
tion of various adsorption data on cellulose (in- 
cluding selective adsorption from solution and 
heats of swelling) concluded that all the estimated 
values in the literature for the surface areas of 
cellulose fell into two groups. The first group 
was about 0.2 m?*/g—‘‘the microscopically visible 
surface’’—and the other group was about 300 
m?/g¢—‘the of the transient capillary 
structure created within the cell walls by the 


surface 


swelling agents’’. 

It is pointed out by the present authors that all 
four of the above factors may operate to account 
discrepancy. There is insuf- 
for evaluating the 


for the apparent 
ficient evidence at 
relative merit of the four hypotheses suggested 


present 


above. 

Plotting the nitrogen and water vapor adsorp- 
tion data according to the method suggested by 
Harkins and Jura [28] resulted in a line or series 
of lines concave to the Y axis rather than convex 
as usually found for other materials. In the two 


instances somewhat lower surface area values 


were obtained using this method of plotting. 


(b) Free Energy of Adsorption 


In view of the possibility that we are here 
dealing with an internal as well as an external 
surface, it seemed of interest to attempt to treat 
both systems and 
to apply the relationship of J. Willard Gibbs [29]. 
This equation is applicable to interfaces between 


systems as heterogeneous 


solid surfaces and vapors; it was written by 


Gibbs as follows: 


o=e,—tn,— uw I,;— JT 1) 


Where: 


o==the work spent in forming a square centi- 
meter of clean new surface, 

e,— surface total energy, 

t=temperature 

n=surface entropy, 

u, chemical potential of the solid at bound- 
ary, 

uwo=chemical potential of the vapor at bound- 


ary, 
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,=surface concentration of the solid, and 
r,=surface concentration of the vapor. 


Under isothermal conditions [30] and under 
conditions in which [, is set at zero (by choosing 
the Gibbs boundary so that the surface concen- 


. . ° ° 
tration of the solid is zero) eq (1) reduces to: 
de ld uo. (2) 


Under appropriate conditions dy, may be replaced 
by its equivalent d (RT In P) where P is the vapor 
pressure, and we have: 


do/d(RT \n P) rs, (3) 


e °P 
[do RT| rdlnP, 


and 


| do RT/MS { "gdm P. 


Where Gibbs’ surface concentration 
pressed [32, 33] as the number of moles of vapor 


rs is eX- 


adsorbed per square centimeter: 
ro=q/MS. 


q=grams of absorbate, 

M=its molecular weight, and 

S=total surface area 

As o is the work required to form a square 

centimeter of surface, it is also the free surface 
energy, therefore, covering the surface with a 
vapor lowers the value of ¢ This change of the 
free surface energy, which accompanies adsorp- 
tion is referred to here as the free energy change 
of adsorption and designated AF as follows: 


AF RT/MS | “q/p dp (7) 


The AF of eq (7) evaluated after eX- 
perimental determination of S, gq and P. The 
integral fq/p dp can be obtained by graphical 
This is done by estimating the area 


may be 


integration. 
under the curve that is obtained when one plots 
values of gp vs P. Plots of the data of figures 
2 and 3 in this manner are shown in figures 6 
and 7. As the plots do not lend 


themselves easily to extrapolation, the method 


(See also [S]). 


used by Boyd and Livingston [32] was employed 


in obtaining the areas. The AF values thus ob- 
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02 03 o¢ 05 o6 
RELATIVE NITROGEN VAPOR PRESSURE P/P, 


Figure 6.— Plot of nitrogen adsorption data in accord 
with a modified Gibbs equation. 


©, Wool; X, viscose; 0, silk; ©, cotton; V, acetate; A, nylon 


(GM WATER VAPOR) 











o6 
RELATIVE WATER VAPOR PRESSURE P/P, 


_ 


Figure 7.— Plot of water adsorption data in accord 


with a modified Gibbs equation. 


>, Wool; X, viscose; 0, silk; O, cotton; A, nylon. 
tained were in turn plotted against the relat 
vapor pressure, P/P,, and the resulting curve 
are shown in figure 8. 

It is noted that the free energy of adsorpty 
of water is the same at all relative vapor pressur 
for wool, silk, viscose, and cotton. The free « 
ergies of the two synthetic polymers, nylon, a 
acetate, did not fall on the common curve 
was the case in the surface area measurement 


\s 


the free energy of adsorption of these two fiber 


Wit! 


was in a separate class (see table 4). 
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02 03 04 O05 06 OF O8 o9 10 
RELATIVE VAPOR PRESSURE P/P, 
RE 8.— Free surface energy of adsorption as a function 
of relative vapor pressure, 
, Wool; X, viscose; 0, silk; O, cotton 
xeeption of nylon and cellulose acetate, whose 
fiber-water contact angles are greater than zero, 
Mhe free energy of adsorption appears to be in- 
sensitive to the chemical nature of the adsorbent. 
It is possible that the bonding between the water 
nud the four fibrous materials may be a hydrogen 
If such bonds existed in the four 
ases, the average value of a variety of hydrogen 


ond type. 


bond energies might very well be the same for the 


our materials. 


Free energy of adsorption of water and nitrogen 


on textile polymeric materials 
W ater Nitrogen 


P) Po=10 


V. Summary 


he sorption isotherms of water and nitrogen 


important textile materials have been 


ned. 
he fiber surface available to the 


areas 
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nitrogen have been computed with the aid of the 
theory of multilayer adsorption. The surface 
areas of the textile fibers were found to vary 
between 0.31 and 0.96 m*/g. These values are of 
the same order of magnitude as the reported areas 
estimated with the aid of the microscope. 

3. The water and nitrogen isotherms on samples 
of the same textile fibers showed marked differ- 
ences in sorptive capacities. 

4. The free energy of adsorption of nitrogen on 
the six fibrous materials have been calculated. 

5. The implications of these calculations are 


discussed. 


Acknowledgment is made for some of the 


laboratory measurements made by Lee A. Dunlap 


and H. J. Valis. 
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Part of the Journal of Research of the National Bureau of Standards 
Sieve Test of Metal Powders 
By Rolla E. Pollard 
Reproducible results were obtained in sieve tests of sponge iron, electrolytic iron, 
electrolytic copper, and nickel when certain variables affecting the sieving characteristics 
of the powders were eliminated or controlled. One of these was a cumulative sampling 
error resulting from repeated riffle cutting of limited powder supplies. Another was the 
effect of exposure of the powder to humid atmospheres. The effects of both variables were 
demonstrated in tests with sponge iron. Variations of considerable magnitude were observed 
when either variable was uncontrolled. 
Significant variations also were noted when the same powders were sieved with different 
sets of certified sieves. These differences were reduced considerably in magnitude when 
comparisons were made on the basis of the size of the average opening as determined during 
the certification tests instead of the nominal sieve opening. 
Development of more effective methods for controlling these variables or evaluating 
their effects would improve the reproducibility of sieve tests. 
I. Introduction [1],' by the National Bureau of Standards [2], by 
the American Society for Testing Materials [3], 
ln most arts and industries producing or using and by the American Standards Association {4}. 
Swdered material, close control of the = size A calibration service is available at the National 
stribution of particles is necessary. Sieving is Bureau of Standards where, for a small fee, sieves 
» most rapid and convenient method of size may be tested and, if acceptable, certified to 
alysis, but it is limited to particles that are conform to the requirements of these specifications. 
tained on a 325 mesh screen, having a nominal Standard methods have been adopted by the 
waning of 44 uw. Manufacturing difficulties have American Society for Testing Materials and other 
evented production of sieves of smaller mesh, technical organizations for sieve analyses of 
nce the size distribution of finer materials must particular materials such as ores [4], refractories 
determined by other means. [5], fine and coarse aggregates [6], pigments [7], 
lwo standard sieve series, with specified powdered coal [8], soap [9], cement [10], and roofing 
mensions and permissible variations, have been materials {11, 12]. Three of these specifications, 
dely used in this country for many years. The namely those for powdered coal [8] and roofing 
vier Standard Seale introduced in 1910 is based materials [11] and [12] include criteria for the 
the size of the average opening, 0.0029 in., in reproducibility of sieving tests by these standard 
1) mesh wire cloth. The openings of the other methods. All three require that duplicate deter- 
es in the series vary by the fixed ratio of the minations by the same operator, using the same 
are root of two. The U. S. Standard sieve sieves, shall check within 1 percent of the total 
es, proposed by the National Bureau of Stand- weight of the sample on all sieves. The standard 
rds in 1919, is based on a 1 mm (0.0394 in.) method for coal [8] requires that duplicate deter- 
pening in the No. 18 sieve. In general, the minations by different operators using different 
es vary in the same ratios as the Tyler sieves shall check within 3 percent on the mate- 
es and there is little difference between rial sieved on the No. 200 sieve. 
Orresponding sieves. The U. S. Standard Sieve spre, Fags - - _— 
is been adopted by the U.S. Government paper. 
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Tentative standard procedures for sieve analy- 
sis of granular metal powders have recently been 
proposed by the American Society for Testing 
Materials (B214—46T the Metal 
Association (5-46T Considerable information 
on sieving is available in the literature, but most 


and Powder 


of it is concerned with materials other than metal 
powders. Investigators generally agree that the 
sieving characteristics of different materials vary 
considerably. In particular, the rate of sieving, 
sharpness of separation, and average particle size 
of sieve fractions are influenced by the shape of 
the particles and their size distribution. The 
diversity of methods used in the manufacture of 
metal powders causes wide variations in particle 
shape and particle size distribution. Conse- 
quently different metal powders or powders of the 
same composition made by different processes, 
might be expected to display variations in sieving 
characteristics, 

It is generally agreed that satisfactorily repro- 
ducible sieving results can be obtained for a given 
material when a single set of sieves and a stand- 
ard method of testing are used. If the sieving 
conditions are not closely controlled, or if different 
sets of sieves are used, variations considerably 
greater than those cited as criteria for reproduci- 
For example, in cooper- 
Wig 


hod 


bility may be obtained. 
ative tests of standard samples of cement, 
and Pearson [13] found differences as high as 13 
percent in the values reported from 80 different 
laboratories, all using standard sieves. Similar 
differences have been reported for cooperative 
tests of foundry sand [14], powdered coal [15] and 
road building materials [16]. Weber and Moran 
[21] in tests with granular and powdered sodium 


bicarbonate obtained differences of as much as 12 


TABLE 1. 


Material Method of production 

Sponge iron Reduced inill seale pulverized, screened and 
annealed 

Electrolytic irot Electrodeposit in sheet form pulverized, screened 
annealed, and again screened 

Electrolytic coppe Electrolytic “Grade B’ 

Nickel Mechanical pulverization 

Zin Atomization 

rin cle 

Lead 

* Hall Standard Flowmeter: (-gram sample, 
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percent in sieving, values obtained with ¢ iffediore. we 
standard sieves in the same laboratory. time 

In beginning an investigation of the basic pfijses. T 
ciples of powder metallurgy at the Nati@were no 
Bureau of Standards, one of the first object™iers wer 
was to determine the effects of particle size one or t 
tribution on the properties of metal powders meas 
the compacts made from them. In the proces@fe useful 
accumulating supplies of sieved fractions 9 within | 
various powders for these studies, signific licroscoy 
variations in sieve analyses were noted wihoo jron 
samples of the same powder were sieved a mills 
different times with the same sieves. Variati@ ticle 
of considerable magnitude also were obtai se pat 
when different sets of certified sieves were hex 
for the same powder. ~ . : 

Therefore, further study of the problems « tig 
nected with the sieve analysis of these powid = 
seemed desirable. The results presented herg “ 
are those obtained in preliminary studies of | ticles © 
effect of variation in sieving time, the reprodug!! ‘ee 
bility of sieve analyses and the particle size dees freq! 
tribution within the sieve fractions. arbon 1 





II. Materials 







The character and properties of metal powde 
used in obtaining the data presented herein 
given in table 1. Four different methods of p 
duction are represented; reduction (iron), electr 
ysis (iron and copper), pulverization (nickel), a 
atomization (zine, tin, and lead). 

The size distribution of the sponge iron, elect: 
lytic iron, copper and nickel powders was fow 
that 


size could be obtained by sieving. 


to be such several fractions of substant 


These powde n 


Fineness (approximate Apparent poy 


density 

gm icc 
32 to 40 percent through 325 mesh 2. 51 
37 to 46 percent through 325 mesh 2. 28 
51 to 59 percent through 325 mesh 2.70 
68 to 73 percent through 325 mesh 3. 65 
99 percent through 325 mesh 2. 02 
8Y percent through 325 mesh 4.19 
84 percent through 325 mesh 5.61 
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. were suitable for studies of the effect of 


time and the reproducibility of sieve 
st The atomized powders on hand at the 
e not suitable for such studies because the 
were extremely fine and sieving produced 
or two small fractions. However, micro- 


measurements of accumulated fractions 
iseful information on particle size distribu- 
vithin fractions. 


revealed that the 





scopic examination 
iron, which had been prepared by reduc- 
of millseale, was composed of irregular plate- 
articles as shown in figure 1. The thickness 
bta se particles appeared to be considerably less 
' either of the other two dimensions as in- 
ed by the relatively small adjustment re- 
oe d to change the focus from the top to the 
‘he m of the particle. Many of the larger 
of yprticles consisted of several plates apparently 
rod together by oxide (fig. 2). The pearlite 
iz s frequently observed indicated the presence 
bon in the metal 








] Sponge iror powder particle s, as received, 


00. 
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The electrolytic iron was similar to the sponge 


iron in shape of particle and oxide content, but 
the plates were more distorted and there was no 
microscopic evidence of carbon in the iron (figs. 3 
and 4 

The electrolytic copper powde was composed 
of two distinetly different types of particles. One 
was plate-like, and the other was irregular and 
tree-like in shape, as shown in figure 5, 6, and 7. 
Many of the spaces between the branches in the 
tree-like particles appeared to be filled with light 
blue or dark blue oxide. 

The nickel particles (figs. 8, 9, and 10 


roughly rectangular or elliptical in shape 


were 
Plate- 
like particles were prevalent, and the larger 
particles often consisted of two or more plates 
apparently cemented by oxide. 

The zinc, tin, and lead powders, produced by 
atomization were composed of globular particles, 
all dimensions being of the same order of magni- 
tude. These particles appeared to be largely solid 
little 
(figs. 11 to 16, inclusive). 


metal with porosity and few inclusions 





Fieure 2 lreas of pearlite are visible in sections of 


sponge iron part cles, 


be composed of several plate emented together 


th oxide Etched with nital x AM 
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or particles retained on 140 Figure 5 Electrolytic copper particles retained on 





mesh sieve x 50. 





Fieure 4 Sections of electrolytic tron particles which shou Ficure 6 Sections of irregular, tree-like copper } 
aiartion ¢ ‘ res i of bd 
! . is a ese yy pu verization Some of the space between branches apparently are filled I 
Some partick ' tes held together by oxide Etched wit! , Etched with ammonium hydroxide and hydrogen perox 
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Figure 9.—Sections of plate-like nickel particles. 


Space s hetween branches of tree-like, electrolytic 
copper are filled with oxide some of the particles were distorted during pulverization and the interspaces 
apparently filled with oxide. Etched with nitric and acetic acids in acetone 


etchant as figure 6. 500 FAM) 


, - ’ Fiaure 10.—Section of nickel particle that apparentl 
Nickel powder particles retained on 140 mesh ; : ; PI — 

5”) a consisis of several plate s cemented togethe r with oxide. 
Sieve < dO. . f 


Same ctchant as figure ¥ < 500. 


ve Tests of Metal Powders 491 





\ hic hy were selec ted from the follow Ing SIX SleVves of 


the U.S. Standard Series 











Figure 14. The globula particles of tin contain fe 
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me h were se ‘lected from the following rsix SleVes of 


the U.S. Standard Series 


si) 247 212 
x) 19 ) 0s 7 
Ww) Ww ‘ 4 a 
2 4 11s " 





‘Two sets of sieves, designated as set I and set I] 
were used for most of the tests. The No. 230 sieve 
was omitted from both of these sets. In later 
tests set No. IL] was used, which omitted the 
No. 80 but included a No. 230 sieve. All sieves 
ead powder retained on were tested at the National Bureau of Standards 
~50 and certified to conform to ASTM Specification 
11-9 and Federal Specification RR-S366a. 

The sieves were brushed with a stiff bristle 






brush after each test. At frequent intervals, for 
; instance at the beginning of a series of tests or 
“ when a different material was to be sieved, the 
¥. sieves were cleaned in hot soapy water with stiff 
bristle brushes, rinsed with alcohol and dried 
an oven at 110°C 
The sampling procedure included the use of a 
riffle-type sample splitter, to reduce the entire 
supply of each metal powder (50 to 100 |b in 
most cases) to a “sample supply” that could be 
stored in one to five 1-pint Mason jars (3 to 15 Ib). 
Smaller samples for individual tests were taken 
from these jars, after thorough shaking, by means 
of a curved spatula (‘‘scoopula’’), which reached 


to the bottom of the container. 


IV. Time of Sieving 





ke 16. Lead particles are roughly spherical with few pe : : A . , 
lo determine the effect of shaking time on 


sieving, 25, 50, 100, and 150 gram samples of 


voids or inclusions. 


tehed with nitrie and ac cids in cerol. <5 
, sponge iron powder (3 samples of each weight 
Ill. Method of Testing were sieved for various periods of time, using 
: Sieve set L. Two series of tests were made, one in 
\ Tyler RO-Tap shaker provided with an March, the other in July 1946. The schedules of 
atic time switch was used in making the sieving time periods for both series are given in 

tests as specified in ASTM B-46T and table 2. 
t45T Standards. The machine was ad- To obtain information on possible differences 
1 to accommodate a set of five tested sieves, caused by sampling, each of the three samples of 
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TABLE 2 Schedules of sreving time for two serves of sieve 
tests of sponge iron—NSieve Set I 
Weigt Numi 
o ‘ ! ind evir 
my sami} 
SERIES I 
ninnul 
10. 20. 30. 40, 10 
“4 5, 10, 2), 30, 40, 60, 20 
TLD 5, 10, 20, 30, 40, 60, SO, 30 
mo ; . 10, 20, 30, 40, 60, SO, 100, 4 
SERIES II 
, 10, wD , 40, 60, 80, 100, 12, lo 
“a 5 10, 2 ), 40, G0. 80, 100, 120, 2p 
lw) 10, 20, 30, 40, 60. 80, 100, 120, 30 
IM” 5, 10, 20, 30, 40, 60, SO, 100, 120, 49 


selected 
After each period the 
weighed sieve fractions were thoroughly mixed and 
for the The total 
weight added to the top sieve at the beginning of 


each weight was resieved through the 
range of time periods. 
resieved next time period. 
each period was used in calculating the percentage 
by weight retained by the sieves during that 
period. 
The last 


schedule was selected as a duplication of one of the 


time period in each series of the 


earlier periods for the purpose of determining to 
what extent disintegration of particles had oc- 
curred as a result of grinding or abrasion during 
prolonged sieving. The total time of sieving, 
including the last or duplicate periods ranged 
from 110 to 505 minutes. 

The results of two series of sieve tests on sponge 
18, 19, and 20. The 
sieving rate was greatest during the initial period, 


iron are shown in figures 17, 


decreased rapidly for periods up to 15 to 30 
minutes, and then decreased slowly with continued 
sieving time. Appreciable amounts of material 
continued to pass through the sieves after sieving 
It is evident, there- 


fore, that for practical purposes at least, there is 


periods up to 120 minutes. 


no definite time of sieving at which complete 
separation can be said to have occurred. 

The period of high sieving rate or rapid change 
in weight distribution represents the time required 
for the definitely undersized particles to pass 
through the sieve openings. The period of slowly 
diminishing sieving rate represents the continued 
the so-called “diffieult”’ 


separation of particles 
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SIEVING TIME MINUTES 


Ficure 17. 
made at different times using sieving schedule giv 
table 2. 


Sieve tests of 25-gram samples of spong 


Series I was made in March, Series II in July. 
(particles whose controlling dimensions are clo 
to those of the sieve openings or which, becaw 
of their shape, pass through the openings \ 
difficulty) or the larg 
particles through oversize openings. 


continued passage of 


The values for the duplicate time periods, give] 


in table 3, indieate that the continuous chang: 
weight distribution during prolonged sieving \ 
not caused by disintegration of the parti 
The values agree within the range of differe! 
between samples, although the sieving tim 
tween duplicate periods was as much as 4} 


Evidet 


disintegration of these sponge-iron particles 


minutes (see sieving schedule, table 2). 


result of grinding or abrasion, during prolong 
sieving 


») 


The average total sieve losses during the su 


is negligible. 
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4c 60 80 
SIEVING TIME MINUTES 


URE 18 Sieve tests of 0 gram samples of sponge tron, 
schedules, for the 25, 50, 100, and 150 gram 
sumples, respectively, were 1.0, 1.5, 1.2, and 1.0 
percent for Series | and 2.0, 1.9, 1.3, and 0.8 for 
Series I. 


imulative losses over as many as 10 different 


Each of these values represents the ac- 


sieving periods, hence these losses are not con- 
sidered to be excessive. As the total losses de- 
creased appreciably with increase in weight of 
sample, particularly in Series Il where the sieving 
schedules were approximately equal, they probably 
resulted from repeated handling rather than from 


lormation ol appreciable amounts of dust. 


V. Reproducibility 


lhe curves of figures 17, 18, 19, and 20, which 
resent the change in sieving rate with increas- 
sieving time are noticeably different for Series 
nd Series II. 

ed for the pan fractions that are measures of 
The differ- 


This difference is particularly 
iaterial passing all of the sieves. 


Sieve Tests of Metal Powders 
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Figure 19.—WSieve lests of 100-gram samples of sponge iron. 
ence between the pan fractions of the two series, 
therefore, represents the algebraic sum of the 
differences between the corresponding sieve frac- 
tions. The values for the corresponding samples 
of the two series fall into definitely separate 
groups after the initial 5-minute sieving period, 
With increased sieving time the mean values for 
the two series are separated by an amount greater 
than the difference between individual samples of 
either series. 

The results obtained for Series | show differ- 
ences between individual samples of the same 
weight, sieved for the same length of time, rang- 
ing from 0 to a maximum of 1.9 percent of the 
total weight of the sample. In Series Il the 
differences between corresponding samples range 
The differ- 


ences between the average values for the two 


from 0 to a maximum of 2.4 percent. 


sevies vary from 0 to a maximum of 4.6 percent. 


It should be noted that these differences are con- 
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differences between the series 



























a | i © INDIVIDUAL SAMPLES SERIES | , two , are = mue| 

| ta. ae greater in magnitude and are consistently n tly 

a T T same direction. It must be concluded, the: efop 
that they were caused by other uncontrol} 
variables. 

40 One of these variables would be an addition 
sampling error resulting from the use of a sump 
splitter to cut the powder supply, original], 

35 lb, to the two different “sample supplies” of 3 
15 Ib from which the samples for the two ser 

8 59 oe were taken, This error, being produced by 

3 mechanical operation, would be expected to | 

y smaller than that resulting from the many 

5 ash | | method used in selecting samples from {| 

Hi 200 we sn “sample supply”’. 

5 Another possible variable is suggested by t! 

y 20 t t t t different seasons of the vear in which the tes fi 
were made, one in March, the other in July 194 
The physical condition or sieving characteristic. ‘ 

‘S| . 
of the powder may have been changed by exposy 
to different atmospheric conditions. A change vf 

at | | | | | size distribution, for instance, could be produc : 

iY j sii elle under oxidizing or corrosive conditions such . 
Si =t=-= === 4 would prevail with high humidities, by the forma. ]j- 

5 a { t tion of additional composite particles similar 1» 9" 
those previously described as being composed | |g 
two or more smaller particles cemented together | 

% 20 40 so _ 26 by oxide (see figs. 2, 4, 6,9, and 10). Variation | 

pe in humidity might also affect the sieving charact 
Ficure 20.— Sieve tests of 150-gram samples of sponge iron. istics of the powder by changing the tendency : 5, 10 
particles to agglomerate, or by modifying |! follow 

sistently in the same direction. In other words, static charges accumulated during handling | Ser 
there is no overlapping of results between the two sieving. 4 
series, and the maximum scatter of corresponding To test the effect of these and other possibl jj 
values for both series considered together (5.3 variables such as clogging of sieves or changes | Ser 
percent) is more than twice as great as that of sieve motion, additional sieve tests were made 
either series alone (1.9 to 2.4 percent). November and December 1946, and January a 

The close grouping of the values for the samples February 1947. For these tests sieve set No. Ill 
of each series and the fact that smooth curves was used. 
can be drawn through the points representing A new hard maple plug was placed in the RO 
them indicate that testing conditions were rela- Tap machine, and the sieve motion was check 
tively stable during each series of tests. The at frequent intervals over the 4-month period | Ser 
separation of the curves by an amount greater counting the number of hammer blows per minut’ | 
than the seatter between the samples indicates The rate of hammer blows was found to average | 
that testing conditions or the sieving character- 154 per minute with variations no greater that | 
istics of the powder were different for the two +1 blow per minute. This fixed the speed 0} 
series. The differences between the correspond- rotary motion of the sieves at 289+2 revolution Th 


ing fractions of the same series are comparatively 
small and can be attributed chiefly to the sam- 
pling error inherent in the method used in select- 


ry 
The 


ing the samples from the “sample supply”’. 
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per minute, 


Samples of sponge iron 25, 50, and 100 gram 


in weight, all taken from the same freshly rifll 


cut 


‘sample supply’, were sieves for periods | 
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TABLE 


Samples were land resieve 


samples 


d to 


i ffect of prolor ged sieving ot 


che 


Last 20-mitr 
period 


ean 


ples 


Dev 


tion 


fron 
iron 


" 
mre 


SERIES 


Sponge 


dules given in tabk 


samples 


i 


1 


SERIES I 


10, 20, 30, 40, and 60 minutes after the 


following treatments of the powder: 


series II] 


3 samples of each weight for each 


time period were oven dried for 
| hour at 110°C 
3 samples of each weight for each 


prior to sieving. 
Series IV 
time period were exposed for 64 
to 72 hours to a very humid at- 
vessel (a 


mosphere in a closed 


desiccator) over water in which 
wicks were dipped to increase the 
evaporating surface. 

The humidified samples of Series 
IV were mixed, 


dried for 110° C 


after sieving 


5 hours at and 


resieved. 


the of 


representing the pan fractions of the 


results curves 


ro DF 


are illustrated by 


vram samples. Similar results were obtained 
the 25 


unidified samples, Series IV, are consistently 


and 50 gram samples. The value for 
than those of the dried material, Series III, 
to 6.5 


nounts corresponding to 3.5 ) percent 
original weight of the sample 
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Ficure 21. Three ies of sieve tests of sponge tron made 
with 100-qram sam ples all from the 
the following treatments of the 
dried: Series T\ humidified; Series | . 


IV dried and resieved 


“sample s upply”’ 
il 


samples of Series 


Same 


powde Ne ries 


atter 


this difference (approximately one-fourth) was 

the 
dried and resieved as shown by the position of the 
It that the 


values for individual samples of Series II] are 


recovered when humidified samples were 


values for Series V. will be noted 


all quite close to the curve. The values for the 


497 











Series LV are more seattered, and the irregularities 
are matched by similar deviations in Series V. 
The relative humidity in the humidification 
chamber undoubtedly was very high, probably 
saturated at times due to overnight or week-end 
Exposure to this 
atmosphere for 64 to 72 hours produced no change 


changes in room temperature. 


in the appearance of the powder except for a slight 
tendency to cake and to retain the shape of the 
dish in which the samples were contained. The 
gain in weight of the 100 gram samples, due 
largely to the presence of absorbed or condensed 
was 0.4 to 1.5 0.8 
A large part of this weight increase was 


moisture, percent, average 
percent. 
lost when the humidified samples were dried prior 
to resieving. The oxide content of the humidified 
samples after drying and resieving, as measured 
by the percentage loss in weight in hydrogen 
according to Metal Powder Association Tentative 
Standard 2—-45T The hydrogen 
loss of samples taken from the original “sample 
supply” was 0.9 percent, indicating that the 
reducible oxide content of the powder was in- 
creased by 0.3 percent of the original weight of 
the sample during humidifying, sieving, drying 


, Was 1.2 percent. 


and resieving. 

All samples of the three series for each sieving 
period were taken from the same “‘sample supply”’ 
at the same time. The sieving for the three series 
was done concurrently, that is, all samples for 
each time period were sieved within the 72 to 96 
hours required to humidify, sieve, dry, and resieve 
the samples of Series [V and V. These precau- 
tions eliminated the consideration of such variables, 
as the sampling error resulting from riffle splitting 
clogging of sieves or changes in sieve motion as 
major factors in causing the differences observed 
between the three series of tests. It must be con- 
cluded, therefore, that these differences resulted 
only from the exposure of the powder to a humid 
atmosphere. 

Two additional series of sieve tests of sponge iron 
were made about 1 month after Series III, IV, and 
V were completed. The same sieves were used, 
but the samples, 25, 50, and 100 grams in weight, 
were taken from two different “sample supplies”’ 
riffle split within the same week. The samples 
were dried for 1 hour at 110° C and sieved for 
periods of 5, 10, 15, 20, 30, 40, and 60 minutes. 

The pan fractions of the 100 gram samples are 
shown in figure 22 in comparison with those of 
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Ficurer 22.— Three series of sieve tests of 100-qram sar p 
taken from three different “‘sample supplies” riffle cut | 
the same lot of sponge tron powder in the order Series 11] 


Series VI, Series VII. 


Series II]. The difference between each succeed- 
ing series, in the order in which they were riffle-cut 
amounts to about 1 percent of the original weigh 
of the sample. The sampling, handling, and 
sieving for all of these tests were done in a heated 
laboratory during the winter months when th: 
humidity and relatively 
Changes in sieving characteristics resulting from 


was low constant 
exposure to such atmospheres should be small 
It is believed, therefore, that the differences be- 
tween Series II], VI, and VII resulted chiefly 
from sampling errors produced by repeated riffle 
cutting. 

Variations similar to those obtained with sponge 
iron have been observed with electrolytic iron, 
electrolytic copper, and nickel in sieve tests made 
over extended periods of time. When a series of 
sieve tests of any of these powders was made on 
the same day or over a short period of time, the 
results for individual samples usually agreed within 


about 1 percent of the original weight of the 
sample. A similar series of tests made at anothe: 


time using the same sieves and the same source of 
powder, again usually showed good agreement be- 
tween the individual samples, but the average 
results often differed from those of the previous 
series by as much as 4 to 8 percent. 

The variations observed in the sieve tests of 
these powders, extending over periods of 1 week 
to 8 months, are given in table 4. The maximum 
differences between corresponding fractions o! 
electrolytic iron were as high as 5.6 percent with 
Sieve Set | (column 2) and 4.4 percent with Sieve 
Set 2 (column 4). With electrolytic copper, the 
differences were 7.8 percent (column 6) and 7.! 
percent (column 8), respectively. 
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\RLE 4 Nieve tests of ele ctrolytic iron, electrolytic copper, and nickel po vders shou ng variation of results 


£ 
} 
: 
“Twi 


{100-gram samples, sieved 30 min 





is? f 
ae 5 Weight pereent retained 
ee | 
- —4 Electrolytic iron Electrolytic copper Nickel 
—_—4 
Sieve Set I (2 Sieve Set II (20 Sieve Set I (33 Sieve Set IT (48 Sieve Set I (7 Sieve Set II (16 
§. ~—vs tests tests tests tests tests) tests 
‘\ 
Deviation Deviation Deviation Deviation Devia Deviation 
Mean from Mean from Mean from Mean from Mean tionfrom Mean from 
— ‘ Tecan mean mean Inean Ihean mean 
1 3 “4 (5 ‘ 7 ® ’ 10 11 12 
| 
80) Pract rrace Trace rrace 
I] 100 2 0.1 l 0.24) Prace Trace 
iW 12.7 j-4).7 w.5 ty 7 6.2 +0. 840.7 4 Oo. U4) S 2 “ 0.1 1.4 0. t-4) 
mW) 21 l 1.0 24.4 a .¢ 16.4 j .& .8 19. 1 1, 8-2. ¢ 11.4 Ol 12 ' 7 
2 “3 l 7 Am) s wt +3. 8-18 2? 9 2.4 12 l | +3 2.4 
ceed- pan 4.0 | +2.1-3.5 94 +2.4-2.0 M1 +2.8-5.0 52.8 | +4.8-2.3 73.4) + .2 70.3 | +2.8-3.7 
("UT 
eleht » These seven tests were made over a period of | week All other tests covered periods of 6 to 8 months 
and _ . . : +, . 
lhe data for nickel illustrate the agreement atmospheric conditions to which the powders were 
“ated ' . , . d ; : 
Bat may be expected in sieving tests made over exposed may have contributed to the variations 
1 the ; . ° ° ° 
short space of time as well as the lack of agree- noted in the sieve test results. A sampling error 
tant “~ . ry , - . ; 
, ent among tests made at different times. The built up by several repetitions of riffle cutting 
rom “ . ° . . i 
i 4 sults of seven tests with sieve Set | made over might also account for part of the differences ob- 
ma re . ' ; 
: period of about 1 week agree within 1 percent served. Further work will be required under 
} De- ore ° . . °° ° ° 
olumn 10). The three series of tests with sieve better controlled conditions, particularly with re- 
refiy ° » . . ° . ° 
an t Il, using 6, 4, and 6 samples respectively, spect to sampling and humidity, before the rela- 
riithe “> . . . ‘ 
ere made at different times over a period of tive effects of the two variables can be determined 
out 8 months. The values within each series with certainty. 
mage : ’ 
@ tests made over a period of a few days, agreed The average values of a large number of tests 
ron “1 . s . ° 
ithin 1 pereent, whereas differences between the of electrolytic iron, electrolytic copper, and nickel, 
de fF 2 ° ° - ° 
ries were as much as 6.5 percent. made with two different sets of certified sieves over 
Ss of ° ° ° . . ‘ A one 
Much of the data given in table 4 were obtained long periods of time, are given in table 4. The 
‘on } ° ° 
uring the spring and summer seasons of the year. differences between these averages (calculated to 
the ; er . ° ° . - 
t is possible, therefore, that variations in the the third decimal place) are compared in table 5. 
Phin 
the nie 5. Sieve tests of electrolytic iron, electrolytic copper, and nickel powders showing differences in results obtained with 
her | two different sets of certified sieves * 
e of {100-gram samples, sieved 30 minutes 
be- 
Weight percent retained 
age . . 
OUS Sponge iron Electrolytic iron Electrolytic copper Nickel 
‘ , Esti- - ” Esti- Mean, Mean, Esti : Mean Esti 
f ps oy Me a, | Differ-| mated Mean, Mean, Differ- mated Sieve Sieve | Difler- mated - = Sieve  Differ-| mated 
oY Ret I Set II ence of standard ‘Setl Set If ence of standard| Set I Set II | ence of standard Set | Set Il ence of standard 
| tests test 2 means) error ot (so. r sts 20 t .ete) means | error of (33 (48 (2means, error of (¢ pacts ii 2means error of 
Pen ; on difference ‘-” “ om difference tests tests difference \‘ tests difference 
uni 2 4 (5 (6) 7 (8 (¥ (10 ll 12 1S 4 15 (16 
ol 7 
race Prace rrace Trace Prace rrace rrace rrace 
ith ' ze 
Trace Pract 1.172 1. 260 0. O88 0.05 Trace Prace rrace Prace 
Ve 7. 267 1 40) 2. 867 0.35 12. 700 10. 495 2. 205 12 6. 152 4.204 1. S58 0.09 2. 286 1.431 0.855 0.10 
} 22. 600) 24. 400 1. 800 ‘5 21. 280 24. 380 + 100 22 16.352 19.065 2.713 23 11.429 12.484 1.055 0 
a 5 167 3. 000 1. 833 7 20. 260 23. 860 + 600 17 20.630 22.927 2.207 29 «12.514 15. 356 2842 77 
; Ss. 400) 7. SOT 5. ll $3. 956 9. 425 4.531 0) 56. 103 2. S21 3. 282 45 73. 371 70. 314 3. O57 SS 
on nickel with Sieve Set I were made over a period of 1 week, tests on sponge iron over a period of 4 weeks, all others over periods of 6 to & months 


> Tests of Metal Powders 499 














Included also are the averages of a few tests on 
sponge iron made with both sets of sieves within an 
interval of 4 weeks on samples from the same 
The maximum difference of the 
means for the two sets of sieves ranges from 2.9 to 


“sample supply”. 


4.5 percent and, although the same fractions are 
consistently larger or smaller for a given sieve 
mesh, the magnitude of difference varies with the 
different materials. 

The standard error or standard deviation of the 
difference of the means, op, given in columns 4, 8, 
12, and 16 was estimated according to the formula 
(ljn+1 S¥?/n+m—2)' 


op 


where 


AX =the arithmetic mean of the observations | 
sieve Set I at a given mesh 
Y =the arithmetic mean of the observations { 


sieve Set I] at a given mesh 


The small magnitude of the estimated standa 


error in comparison to the difference of the meay, 











indicates that the difference between sieve s op 
was significant in spite of the wide scatter in in ts ch 
vidual values for tests made at different times, ame 

Two new 325 mesh sieves were obtained fr eet 
the manufacturer of Sieve Set | for use in eit! vith S 
of the original sets if replacements should be nee Ine 
sury. These sieves were tested and certified | sting 
the National Bureau of Standards to conform | everal 


materi: 
The screen wire cloth was of monel metal inste; 


all specification requirements except 


luring 





, The 
n—the number of tests with sieve Set I , aol 
f th si a of bronze, the latter being unavailable at the ti n 325 
m==the number of tests with sieve Set , “ar : a asl set 
— of purchase. The differences in sieving resul) 44 
y Y y 2 n Y 2 . . a. i ~ . 7 7 l 7 > ] a“ 
- ~- A obtained when the new sieves were substituted nclusi 
Sy?=SY*—mFY? each of the original sets are shown in table (J...) 
~X*=the sum of the squares of the observations — The tests with the new sieves (columns 3 and} on, , 
i 
for sieve Set [at a given mesh were made on the same day. Those with th! 
S¥?=the sum of the squares of the observations — original Sieve Set I and Sieve Set Il (column Tar 
for sieve Set IL at a given mesh were made 2 weeks before and 2 weeks lat 
TABLE 6.-—Sreve tests of sponge tren showing variations due to use of different sieves in same set 
[The tests were made over a period of about 4 weeks in July Average of 3 tests on 100-gram samples sieved 30 minutes} 
Weight percent retained 
With original 325 With new 325 mesh With new 325 mesh 
mesh sieve sieve (NBS No. 8887) Difference sieve (NBS No. S888) Difference Difference 
Sieve from from between 
No original original new sieves 
aie Deviation 2 Deviation (3)—(1) - Deviation (6)—() (6) — (3) 
Mean | trom mean!) “8? | trom mean Mean | trom mean 
2 (3 (4 (5 (6) (7 (s @ ; 
; ' 
SIEVE SET I 
si) Prace Trace Trace 
Th) rrace Trace Trace 
i” 7.2 O12 7.0 +) 1-0. 2 os “9 +i) 1-4) “4 
aM) 22.6 I 22.4 +.2-.4 2 22. 6 +. 1 o 
25 1.2 4 3s. 2 +.1 +70 41.0 +.4-.2 +08 +258 
Pan 4 2-. 1 2.0 +. 2-.1 4 ay +.3 “5 ~3.1 
SIEVE SET Il 
sw Prace lrace Trace 
lon race rrace rrace 
iw 14 +06 5.2 +0. 1-0 +i s 5. 1 0 +07 
a) “44 » © 23.5 4 1-0 4 23.3 +1 —1.1 
2h ‘ 2 is 6 + 2 +5.6 wos + 1-0 +78 +22 
Pa 7.9 2 $2.2 17 2 1 7.7 2.2 
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All samples were taken from the 
ame “sample supply”. It will be noted that the 
tions larger than 200 mesh agreed within 1.2 


The new 


esp tively. 


rit 
oreent With any given sieve assembly. 
5 mesh sieves retained 7.0 and 9.8 percent more 
rial, respectively, than the sieve they replaced 
ve Set I. With Sieve Set I], the differences 
vere 5.6 percent and 7.8 percent with correspond- 
The amounts 


nal 


no 


ng changes in the pan fraction. 
tained by the two new sieves differed by 2.8 
wreent when used with Sieve Set [ and 2.2 percent 
vith Sieve Set Il. 

li connection with these differences it is inter- 
sting to compare the average openings of the 
everal sieves. These sieve openings, as measured 
luring the certification tests, are given in table 7 

The variation in the average opening permitted 
n 325 mesh sieves by ASTM Specification E-11 
s 44 microns+7 from 41 to 47 u 
nelusive. The average openings of the four 325 


percent, or 


nesh sieves range, in the order given in table 7, 


rom a value near the upper limit to one at the 


Average opening of sieves as measured in cer- 


tification tests 


Average opening @ 


NBs 
“Sieve No Certi‘ica 
tion No Bet ween Bet weer 
warp wires shoot wires 
ORIGINAL SIEVE SET I 

a Mo 
~ SSO; Ist} 174 
mo SAL) 143 14s 
iw SNOT 107 102 
a0 SSO 7s -F 
$25 be a 1) 44) 


ORIGINAL SIEVE SET I 


sD SARU 5 74 

100 ASM) i44 140 

140 RRO] 104 107 

an SAUD 77 7 

25 SAU 43 45 

NEW NO. 325 SIEVES 

2 SAAT 4} 4 

25 SASS 4} i2 
easturements were made under the supervision of L. V. Judson 
logy Division of the National Bureau of Standards who states 


surements were made with an accuracy sufficient to determing 


t the openings were within the limit ermitted by 


It is not likely that the errors are in eveess of about | a 
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lower limit. The weight of the fractions retained 
by these sieves vary in the reverse order, as would 
be expected, the sieve with the largest average 
the quantity, the 
sieve with the smallest average opening retaining 
The effect of seemingly 


opening retaining smallest 
the largest quantity. 
negligible differences in sieve aperture on sieve 
tests is illustrated by the results obtained with 
sieves S887 and S888. Although the sieves differ 
by only 1 w between the shoot wires, the fractions 
retained by Sieve 8888, with the smaller aperture, 
were 2.2 to 2.8 pereent greater than those retained 
by the other sieve. 

Variations of such magnitude in sieving results 
with different sieves are not peculiar to metal 
powders; similar variations observed by many 
workers with nonmetallic powders were men- 
tioned in the introduction to this paper. 


VI. Size Distribution within Fractions 


The particle size distribution of the fractions 
retained on the several sieves was determined by 
Xylene containing a 
few drops of cedar wood oil was used as a dispers- 
ing agent in making the slides. Measurements 
were made on the image of the particles projected 


microscopic measurement. 


upon the ground glass of a camera at a magnifica- 
tion of &50. The projected dimension measured 
was the one which, in the opinion of the observer, 
would control the passage of the particle through a 
sieve opening. 

The controlling dimension was easily determined 
when, as was usually the case, the projected image 
was circular, roughly rectangular or elliptical. 
However, when very irregular particles, such as 
the ‘‘tree-shaped” electrolytic copper particles 
shown in figure 5, were measured, selection of the 
controlling dimension was difficult. In such cases 
the branches were considered a part of the dimen- 
sion when they were small or close together, as at 
A and B (fig. 5). When the branches were widely 
separated, as at (’, the dimension was measured 
from the fork of the branch. 
mensions obtained for the three particles by this 


The controlling di- 
procedure (between 5 and 6mm at * 50 magnifica- 
tion) are in line with the dimensions of the other 
particles of the sample. 

The particle size distribution of the five sieve 
fractions of sponge iron larger than 325 mesh is 
shown in figure 23. The fractions retained by the 
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+80 MESH FRACTION 
PARTICLE COUNT—- 342 
AVERAGE PARTICLE SIZE - 267) 


-80+100 MESH FRACTION 
PARTICLE COUNT - 502 
AVERAGE PARTICLE SIZE - 226) 


-100 + 140 MESH FRACTION 
PARTICLE COUNT - IlI7 
AVERAGE PARTICLE SIZE — I67» 


-!140+200 MESH FRACTION 
PARTICLE COUNT - I734 
AVERAGE PARTICLE SIZE - !27p 





-200+325 MESH FRACTION 
PARTICLE COUNT - [748 
AVERAGE PARTICLE SIZE - 95y 








Ficure 23. 


X50 (accuracy 


80 and 100 mesh sieves contain particles ranging in 
The 
range in size distribution decreases considerably 
with the finer 

— 2004-325 


size is over 8Ou (from class mark 60 to class mark 


size from class mark 100 to class mark 340g. 


sieves. However, even in the 


mesh fraction the range in particle 


140u) with over 50 percent of the particles between 
90 and 110g. 

The average particle size of a sieve fraction is 
often described in the literature as the mean of the 
nominal openings of the sieves which retained, and 
which just The 
particle size of sieved fractions of the metal pow- 


passed the fraction. average 


ders listed in table 1, and the mean apertures of 
the passing-retaining sieves are given in table 8. 
Measurements were made at magnifications of « 50 
100, 


respect ively. 


+ Su, 
The average particle size of the 


and with accuracies of +10 and 
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Particle size distribution of sieve fractions of sponge iron measured microscopically at a magnification 


+ 10 microns). 


TARLE 8. 
powders as measured microscopically at magnifications 
<100 with accuracies of +10 and 


Average particle size of sieve fractions of met 


<50 and +5 microv 


respectively 


Mean 
aper 
ture of 
— Elect re 
ing re- Sponge 
taining iron 
Sie Ves 


Average particle size, controlling dimension 


Sieve 
fraction Electro 
lytic lytic 


iron copper 


Nickel Zine * 


>~—™) 
SO+ 100 
100+ 140 
140+ 200 
200+ 325 


* Measured at < 100 ler Measurements made at K 50 magnif 


fractions vary greatly with the different materials 
The values for the atomized powders, zinc, tin, and 


lead, composed of roughly spherical particles with 
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equal, are 
The 


<¢ or very irregularly shaped particles of 


ee dimensions approximately 
» the mean of the sieve apertures. 


«trolytie iron, sponge iron, copper, and nickel 
ve controlling dimensions considerably greater 
an the mean of the apertures. This tendency 
the projected dimensions of particles, particu- 
ly irregularly shaped particles, to exceed the 
rage aperture of the passing-retaining sieves has 
18, 19, 20}. 


en noted by other investigators [17, 


VII. Discussion of Results 


The results of the sieving tests indicate that the 
ving rates of metal powders vary considerably 
bth different sieves and with different materials, 
at appreciable variations in results may be ob- 


ined when the same powder is sieved at differ- 
{ times with the same sieves and that sharp or 


Diiplete separation into fractions cannot be at- 


ined by sieving, at least in a reasonable length 
time. Hence, selection of the proper sieving 
ne for a given material would appear difficult. 
owever, the results also show that reproducible 
ues can be obtained when certain sampling 
rors and other variables affecting the physical 
ndition of the powder, such as exposure to high 
imidities, can be eliminated or controlled. 
nder these conditions the scatter in values for 
mples of the same weight sieved for the same 
« apparently is independent of the time of 
ving (see figs. 17 to 22). Any time between 5 
inutes and 2 hours could be selected without 
riously affecting the reproducibility, provided 
As prolonged siev- 


e other conditions are fixed. 
g would increase the proportion of particles 


ussed by oversized openings, it would seem ad- 
sable to select a sieving time near the end of the 


‘riod of rapid sieving rate, that is, after the 


Pparation of the definitely undersized particles. 


same powder. 


he time required for sponge iron would be 15 to 
) minutes, depending upon sample size and 

her conditions. 
Two variables seemed to cause the most serious 
lferences observed in the results of sieve tests 
ut different times with the same sieves and 
One of these was a sampling 
roduced in riffle cutting, the other the result 
wsure of the powder to high humidity. 
mple error in riffle cutting can be attributed 
to the 


loss of fines as dust during the 
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operation. Evidently this is not a good method 
to use for sampling limited supplies of metal 
powders, particularly if the cutting operation is 
to be repeated, for then the sampling error be- 
cumulative. A 
possibly the use of a sample thief, might improve 


comes noncumulative method, 
the reproducibility. 

The effect of exposure to humid atmospheres on 
the sieving characteristics of some metal powders 
as demonstrated by the tests on sponge iron and 
by the results on electrolytic iron, 
nickel, difficult 
With wrought 
metals it is known that at ordinary temperatures 


suggested 
electrolytic copper, and poses 


problems of storage and handling. 


and low humidity corrosion is very slow or non- 
for 
instance, have been stored for many years in an 


existant. Polished samples of sheet steel, 
atmosphere of about 50 percent relative humidity 
of calcium nitrate) 
However, with finely 


(over a saturated solution 
without evidence of rust. 
divided metal powders oxidization or corrosion 
may proceed at much lower humidities. Relatively 
small temperature changes also may have an 
effect. For example, it has been noted that some 
very fine powders, notably electrolytic copper and 
tin (89 to 91 percent through 325 mesh) change 
color when heated at 110° C for drying. 

The problem posed by the differences in results 
with different certified 
There seems to be little hope that 


obtained sieves also is 
difficult. 
sieves can be manufactured, by present methods 
at least, to tolerances sufficiently close to eliminate 
the differences caused by minute variations in the 
dimensions of the average openings. The long 
continued and successful use of sieve analyses in 
the production and utilization of powdered sub- 
stances shows that present manufacturing toler- 
ances and calibration methods are adequate for 
control purposes at least. For research work, a 
more accurate method of evaluating the differences 
obtained with sieves that vary within specifica- 
tion tolerances would be desirable. 

It should be pointed out that the variations in 
results obtained with different sieves are not as 
great as is indicated by comparison of values on 
the basis of the nominal opening. For example, 
in figure 24, cumulative curves of the data given 
in table 6 are shown in which the percentage by 
weight passing the several sieves is plotted against 
the size of the average openings as determined 
will 


during the certification tests. It be noted 
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Figure 24. Cumulative weight distribution curves for 
sponge iron using the actual size of the average opening of 
the different certified sieves instead of the nominal opening 
as the bases «f comparison. 
100-gram samples sieved for 30 minutes 
that the differences between sieves are consider- 
ably less than they would be if the comparison 
were made on the basis of the nominal openings. 
Possibly still closer agreement would be found 
if the variations in size of sieve openings were 
taken into consideration after the metho of Weber 
and Moran [21]. This method consists in measur- 
ing random openings microscopically. 
are said to be identical if the average openings and 
The per- 


Two sieves 


the standard deviations are the same. 
centage standard deviation or coefficient of vari- 
ance V=100 o/X, where o is the standard devia- 
tion and X the average opening, is used as the 
An empirical for- 
mula based on the sieving time and the coefficient 


measure of sieve equivalence. 


of variation has been applied to evaluate the effec- 
tive openings of different sieves used in sieve 
analyses of soda ash and sodium bicarbonate. 

In the cement industry, the No. 200 sieves used 
in accordance with ASTM Standard C184—44 for 
fineness tests of hydraulic cement [10] are cali- 
brated against a standard sample of cement pre- 
pared by the National Bureau of Standards. 

Hateh [17], using crushed limestone, measured 
the dimensions of sieved particles microscopically 
and, by statistical methods, calculated the size 
distribution by weight. The calibrated size of a 
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sieve was taken as the median or geometri¢ ae 
size by weight of the material retained x ; 
sieve. 

Fagerholt [22], working with ground sand 4 
feldspar, found that the particle size to whic! 
result obtained by sieving in the time, f, should 
referred in determining the cumulative weig 
distribution is the average particle size of 
fraction passing the sieve by continued sie 
under the same conditions to the time 3f. 
particle size is determined by counting the num 
of particles in a weighed portion of the fract 
obtained during the interval ¢ to 3¢. 

It is possible that these or similar methods 
calibration could be applied to sieve analyses 
metal powders. 
in particle shape and particle size distribut 


Because of the wide variati: 


resulting from many different production method 


it is not to be expected that any one formula 
standard material or calibration factor can 
devised that will be equally applicable to all m 
powders. However, even .imited application 
powder from a given source or possibly to 
stricted classes of powders, such as iron made fr 


$ 


reduced mill scale or electrolytic copper or atow 


ized zine, would be of value. 
and limitations must be determined by furt! 
investigation that should include further stud) 
the effect of humidity. 


VII. Summary 


In sieve tests of sponge iron, electrolytic in 


electrolytic copper, and nickel, reproducible : 
sults were obtained only when certain varia)! 


were controlled. One factor was a cumulat 
sampling error resulting from repeated riffle « 
ting of limited powder supplies. Another varia! 
affecting the sieving characteristics of the powd: 


was exposure to different atmospheric conditiot 
The effects of these variables was demonstrat 


by tests on sponge iron. 


successive riffle cut samples were of the order o! 


percent of the original weight of the sample. T 
effects of humidity were considerably great 


differences of 3 to 6 percent being observed betwe 
sieve tests of dried samples and samples expo~ 


to high humidities. 


Variations of considerable magnitude were o'- 


served in the results of sieve tests of the sa! 


metal powder when sieved under compara! 
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ms with different sets of certified sieves. 
ponge iron, these differences were as high 
percent when comparisons were made on 
sis of the nominal sieve opening. When 


mpared on the basis of the actual size of the 


erage Opening as measured during the certi- 
ation tests, the differences were reduced con- 


erably. 
Further investigation of possible methods of 


iminating, controlling, or evaluating the effects 
all of these variables is desirable. 
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« preparation of graphs; and to Eleanor M. 
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Table and Mollier Chart of the Thermodynamic 
Properties of 1, 3-Butadiene’ 


By Cyril H. Meyers, Carl S. Cragoe, and Eugene F. Mueller 


This paper contains a working table and Mollier chart of the thermodynamic properties 


of 1, 3-butadiene in engineering units, for which the data were calculated from a set of 


empirical equations. 


These equations have been published elsewhere and shown to be 


consistent with exact thermodynamic relations and to represent the observed data within 


the limits of the experimental accuracy. 


fluid over the temperature range extending from the triple point 
305.6 


20 degrees of the critical temperature 


use of the table and Mollier chart (log P versus H). 


I. Introduction 


The extensive use of butadiene in the produc- 
tion of synthetic rubber led to a need for more 
complete and more accurate data on the proper- 
ties of this substance than were available when 
it xequired commercial and military importance. 
Letter Cireular LC710, a collection of the avail- 
able data at the time, was published in December 
1942. Further need tables of 
liquid densities of butadiene and other C, hydro- 


developed for 


carbons found in commercial C, mixtures, for the 
purpose of computing, from observed volumes, 
the quantities of materials bought or sold. This 
problem was brought to the National Bureau of 
Standards in July 1943, by representatives of the 
Letter Cireular LC 
736, containing tables of liquid densities of eleven 


Rubber Reserve Company. 


hvdrocarbons found in commercial C, mixtures, 
was prepared from data available in the literature 
November 1943. The need for 
tables of the thermodynamic 
properties of butadiene led to a program of ex- 
at the National Bureau of 
The results of 


and issued in 
more extensive 
perimental work 
Standards to supply the data. 
these measurements were published in 1945. in 
RP1640 and RP1661. The 


sent paper gives the data in the customary 


Research Papers 


neering units in the more useful form of an 
nsive table and a Mollier chart. 


was also published as a Technical Report to the Office of 
Reserve, June 30, 1947 


Thermodynamic Properties of 1, 3-Butadiene 


The table covers the properties of the saturated 


164.05 °F) to within 
Ff). Examples are solved to illustrate the 


Several other charts are also presented. 


II. Experimental Data 


The experimental data upon which these tables 
are based are those given in two recent publications 
ia, oe 
measurements made on butadiene at the National 
Bureau of Standards, compare them with the 
relatively small amount of data available from 
and develop empirical equa- 
with exact thermo- 


These publications describe the results of 


other laboratories, 
tions that are consistent 
dynamic relations and that represent the observed 
data within the limit of experimental accuracy, 
which for P-V-T data on the vapor is one to three 
parts per thousand, for density of the liquid a 
very few parts in ten thousand, and for the excess 
of the enthalpy or entropy over that for the solid 
at the triple point is approximately one part per 
thousand. This limit is probably as small as is 
practicable in view of the rather unstable charac- 
ter of the material that polymerizes to some extent, 
even while observations are in progress. 

Other publications |1, 2] were consulted, and the 
data they contained were used for comparison 
and verification, but not to any considerable 
The temper- 
ature of the triple point was taken as —164.05°F 
(—108.92°C), that of the boiling point as 24.06°F 
(441°C), and that of the point as 
305.6°F (152°C). The value for the critical point 


extent as source data for the tables. 


critical 


is believed to be reliable to one or two °C, Poly- 
Figures in brackets indicate the literature references at the end of this 

paper, 
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merization occurring before and during the obser- 
vations limited the accuracy of the determination 
of the critical temperature. 


III. Notation 


In general, the abbreviations approved as 
American Tentative Standard (ASA Z101-1932) 
have been used, with the exception that in giving 
numerica! the degree 


values of temperatures, 


sign (°) has been retained. The units an/1 symbols 


used are tisted in table 1. 


TABLE 1.—-l nits and notation 
Quantity Unit Symbol Formula 
Pressure Ib in P 
Specific volume ft Vib t 
Temperature F t 
Absolute tem perature RK T T =t+459.688 
Heat added Btu/Ib q 
Internal energy do " “du =dq—pdr 
Enthalpy (heat content do h h=u+pr 
Entropy Btu lb’ R * ds=dq/T 
Latent heat of vaporization Btu/Ib hve hig=h,—h 
SUBSCRIPTS 
Property of liquid 
Property of vapor gy 
Transition from liquid to vapor fg 


* The factor | (ft®) (bin 
tions In engineering units 


0.185052 Btu is necessary for numerical caleula- 


IV. Fundamental Units and Constants 
1. Temperature Scale 


The Fahrenheit seale used in these tables is 
derived from the international temperature scale 
[3] of 1927 by means of the relation: 


Fahrenheit temperature = 1.8 centigrade 
temperature + 32. 


It is not expected that any changes that may be 
made when the pending revision of the seale of 
1927 is adopted, will require any material change 
in the numerical values of temperatures as given 
in these tables. 

Temperatures on the absolute Fahrenheit seale 
were obtained by adding 459.688° to the tempera- 
tures as defined above. 








2. Heat Units 


































All of the calorimetric measurements were ma 
in terms of the international joule, which wy 
therefore the fundamental heat unit. 
upon the standards of resistance and electromotiy 
force maintained at the Bureau and upon the meg 
The relation between the interng. 
tional joule and the corresponding egs units, a. 
cording to the most recent data [4] is: 


It is dabsey 


joule 
} 


solar second. 
using 
olVel 


mult 


1 international joule= 1.00017 absolute joules 


wre th 


° ° ° iV be 
In the thermodynamic equations relations occw 


between quantities of energy, some of which wer 
measured in international joules, whereas other 
such as the product of pressure and volume cag! 
expressed in appropriate units, were measur 
directly in mechanical units, that is, absoluy 
joules. Although the difference might be cop. 
sidered of negligible importance, it was taken int 
account where it could have any appreciable effec 
on the results obtained. 

As secondary heat units the calorie and t) 
British thermal unit (Btu) have been used. Th 
calorie is the International Table (IT) calory 
widely used in engineering tables of the properties 
ol steam. It was officially adopted by the secon 
International Steam Table conference [5] held i 
London in 1929, and is defined as 3,600/860 (ap- 
proximately 4.1860) international joules. In this 
way the mechanical equivalent of heat, so sig- 
nificant in earlier tables, loses its importance as |! 
becomes a defined rather than an experimental] 


volun 
rela 
S5052 


value 


und \ 
lue of 
ye vali 
The v 


e stal 





ie low 


hes 0 


constant and therefore is in no way dependent| 
upon the results of past or future measurements 0! 
the properties of water. 

The Btu used is derivable from the calorie } 
use of relations that depend only upon the relativ: | 
sizes of the centigrade and Fahrenheit degrees 
and the relation between the kilogram and -" 
pound thus, 


1 Btu 


1054.866 int joules. This method of defining the” 
Btu retains the convenient relation 


o . 
453.592 4=251.996 IT calories | 


. 


1 calorie per gram=1.8 Btu per pound. 
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;. Relations Between Energy Units 5. Units of Volume and Density 


elation between the Btu defined above. The volumetric apparatus was calibrated either 
the standard foot-pound is 1 Btu=1054.866 — with water or with mereury. The volume unit 
: 10° 0.3937 778.160 ft Ib was defined by the relations: 
980.665 x 453.592 120° ™ , One cm? of water at 4° C under 1 atm pressure 
» the tables, quantities of heat or work that has a mass of 0.999973 ¢. One ml of water at 
nin Btu may be converted into foot-pounds 4° C under 1 atm pressure has a mass of 1.000000 
multiplying by the factor 778.160. In cases grams. One cm® of mercury at 0° C under 1 
wre the reverse conversion is desired, the relation atm pressure has a mass of 13.5951 ¢. Densities 
vy be written, measured in g/ml were converted by means of the 
relation [4], 
| ft-lb=-0.00128508 Btu. 
1 g ml =62.42658 Ib/ft 
volumes are given in ft® and pressures in lb/in.’, 
relation 1 (ft*) (Ib/in.2)= 144 (0.00128508) V. Description of Table 

985052 Btu becomes necessary wherever numeri- 
| values are given in an equation involving both Only the properties of the saturated fluid have 
ergy in Btu and the product pe in British engi- been tabulated. These are presented in table 2, 
ering units, which uses the notation described in the preceding 
section and in which the argument is temperature. 
4. Units of Pressure The table covers the whole i eh range 
In the British system the unit of force is the from the triple point 164.05° F to the critical 
ind weight, which becomes definite when the temperature 305.6° F. Values are given every 
lue of the acceleration due to gravity is specified. degree Fahrenheit in the temperature range 0° to 
vw value used here is 980.665 em/see~*.* 200° F and at less frequent intervals outside this 
The unit of pressure used in this publication is range. Such arrangement gives a table with an 
« standard pound weight per sq in. absolute. ample number of entries in the portion most used, 


« lower pressures have been expressed also in and keeps the table within a reasonable size. 


hes of mercury with use of the factor, The pressures are absolute and for the lower 
temperatures are given in two units, namely |b/in.* 

1 lb/in.2=2.03601 in Hg. and inches Hg at 32° F. The next two columns 

or are devoted to the specific volume of liquid and 


¢ International Committee of Weights and Measures recom, . , = . » 
idoption of an experimental value 980.665 cm/sec ? (32.174 ft/sec 2) vapor respectively in ft */lb. In that portion ot 
f gravity at latitude 45° and sea level [5 Since then, those who the table where onlv one column of pressures is 
nto use a standard value for g have been divided into two groups, - 


h uses the current “best” value for g at latitude 45° and sea level given, a column of densities of the vapor in lb, ft 
t value is about 980.62 cm/sec 2? Ihe other group has loyally fol has been inserted so that the table contains the 


ead of the international committee in using 980.665 cm/sec rhe 


3 


same number of columns throughout. 


m/sec ? has been used in computing these tables. 
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021117 
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o21279 
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O21 787 


021876 
021964 
022054 
22144 
022236 
022328 
2242) 
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022606 
22702 
22709 
022838 
O22S7s 
02217 


022957 


22007 
023037 
023077 
OZ31L1S 


O2Z315s 





OZS2s1 


(23325 


023364 


Or2s408 
(23448 
O2Z3491 
0723533 


23576 


(230819 


023662 


(23882 
ozs 
O2e7 1 


oO2401¢ 


mon 
. (24084 
O24107 
(24130 


024153 


Properties of saturate dl, 


volume 


Enthalp 


Vapor, r, | Liquid, hy | Evap.,h 


ft 3/ib 
5, 706 


4,44 
4, 310 
2, 461 
1,85) 
1, 406 


1, 079.6 
836.58 
i4. 6 
516.5 


411.0 


45.93 


Ptuilb 


a» 


122 


im 


126.7 
12s. 4 


131 


133, & 


135. 
138 
140. 
142 
145, 
147 
14y 
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165 


168, « 
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Ptu/b Ptu/lh 
282.7 41.8 
219. 2 M15 
218.3 42.7 
217.1 M35 
216.0 445.0 
214.5 346. 1 
213.7 47.3 
212.6 45.4 
211.5 340.6 
210.4 350.8 
2. 3 352. 0 
208. 2 353. 2 
7.1 354. 5 | 
aM. 855. 7 
25.0 a6. 9 
au8.4 458, 2 
aeY 350. 5 
ws 10.8 
200.8 362. 0 
199. 7 363. 3 
198.7 v4.7 
197.6 6. 0 
106. 6 67.3 
196. 2 367.8 
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144 you 5 
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194. 1 370 
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192. 0 373. 3 
191. ¢ 173.5 
191.2 v4.4 
190.7 74.9 
100.3 5 
sv. 1) 
189. 5 
sv. 1 1 
Iss 6 7 
TSA. 2 s78. 2 
187.5 17S. 5 
187.4 174.4 
186.9 174.9 
16. 5 wo) 
186. 1 ‘Ar 0 
ISS. R16 
185, 2 in? 2 
184.8 382.7 
184.4 sR. 3 
183.9 343.0 
183. 7 4.2 
183. 5 is4.4 
183. 3 84.7 
183. 1 mA. O 


3-butadiene 


y 


Vapor, 4, 














Entropy 


Liquid, sy 


Ptu/th°R 








6338 
t408 
_ 47S 


6546 
6613 


A680 
6745 
6sl0 
674 


6438 


Tun 
7124 
7iM 


7245 





7674 


THT 


~ eo 


THY 
71 
7TOH3 


Tus 
. aes 
. 8006 
SOLT 


sus 
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1. 3317 
1. 3317 
1. 3256 
1.3174 
1. 3006 
1. 3023 
1, 2053 
1, 2886 
1, 2823 
1. 2764 
1. 2707 
1. 2654 
1 

l 

1. 260 
1. 2406 
1, 2425 
1, 2386 
1. 2350 
1. 2315 
1, 2283 
1, 2252 


1, 2212 
1. 2201 
1. 2190 
1. 2180 
1. 2170 
1. 2160 
1. 2150 
1. 2141 
1. 2132 
1. 2123 
1.2114 
1, 2105 
1. 2007 
1, 2080 
1. 2081 
1. 274 
1. Die 
1. 2050 
1. 2052 
1, 204 
1. 2038 
1. 22 
1. 2025 
1 2e 
1. 2013 
1. 207 
1, 2002 
1. 1006 
1, 1991 
1. 1085 
1. 1083 
1. 1080 
1. 1978 
1. 197 
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TABLE 2 


Abs. pressure Specific volume 

















era 
p Liquid, t Vapor, r, | Liquid,hA 
I hin He ft2it ft 2/1 Ptu/lt 
9.5 9. 43 0. 024176 9.413 202. 44 
6 9. 770 19. 80 o24199 9, 208 202. 95 
7 10. 003 20. 37 024223 0. OOS 203. 4 
s 10. 241 20.5 024246 &. 814 203. 96 
’ 10. 482 21.34 024269 8. 625 204.4 
10 | 10. 728 21.54 024293 8. 441 24. 97 
11 0. 979 22. 35 024316 8. 261 205. 48 
2 11. 234 22. 87 024340 &. ORT 205. OS 
13 | 11. 404 23. 40 024363 7.917 206. 49 
14 11. 759 23. 4 024387 7.7 207. 00 
2. 029 24.49 v244il 7. 589 207. 0 
f 12. 308 25. 9 024435 7. 431 208. 91 
17 12. 5&2 25. 62 024459 7. 278 208, 52 
1s 12. 866 26. 2 024483 7. 128 200. 03 
If 13 5 26. 7S 024507 6. GS2 200. 4 
2» 45 27. 38 024531 6.5840 210.6. 
2 7 27.99 024555 6. 701 210. 57 
22 1.0 2s. ¢ 024579 6. 566 211.08 
2 14. 3 29.°24 024604 6. 444 11.59 
24 14. 6s 20. SS 24625 6. 308 212. 10 
9 x 0. 54 024653 6.179 212. 62 
2h 15. 32 31. 20 24677 6. O57 213.1 
27 Ld) s 937 213. ¢ 
~ l ” 2 ht 2 14.16 
24 16. 34 2t 706 214. 6S 
i) 97 24777 4 215.19 
1. 69 (24802 487 215. 71 
2 7.4 42 024827 SI 216. 23 
7. at 024852 277 216.7 
‘ 8 6.92 024877 ! 217. 27 
s he) 02490) O77 217.79 
t 5. & 8.4 (24928 4. OR] 18. 31 
2 ). 27 024954 4. S87 18.8 
s Os 40. ( 024979 4.79 219.3 
’ 20. OS 40.80 02500) 4.70 219. 87 
4 2. 49 | 2 025031 4 617 2). 40 
1 2.4 42 O25057 i l 
42 2 45.4 02.508 4.447 
‘ 2 f 4.0) 025109 4. 3 
14 22.19 45.19 0251 4. 28 222. 50) 
i 22. ¢ 46. OS 4. 7 223. 02 
i 25. US +; . iM 4. 150 22 
17 2 i 47. 92 4.{ 224. OS 
ss 24. is. ™ ms 24 ‘ 
iv 4.4 419.58 8. Vl 225. 1 
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values of enthalpy for the liquid have been 
to hundredths of a Btu/lb where they were 


TABLE 2.—Properties of saturated 1, 3-butadiene 
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s ipported by experimental data. 


of the rounding off of the values for the 
or and for evaporation, the latter may differ 


one in the last place from the difference between 


nthalpy of vapor and of liquid 


The values in the table and chart 


Density 
Vapor l/t, 


Liquid hy 


Rtu tbh 


In conse- 


were caleu- 


Enthalpy 


Evap. hye 
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Vapor &, 


thermodynamically consistent. 
In consideration of the large size of tables for 


rempera- 
ture, ¢ 


In consequence, the values given in the table are 


superheated vapor and of the fact that the data 
for 1,.3-butadiene ean be represented on a chart 


nearly within the accuracy of the observed data, 


the properties of the superheated vapor are pre- 


by methods deseribed in a previous publica- 


S|] after first converting the four empirical 


ms given in that publication to engineering 


mputing the figures in the table from the 


al equations, 


its bevond the number to be finally retained 


nces 


ept to a minimum, and exact thermody- 


relations have been used whenever possible. 


all ecaleulations were carried 
table. The values have been checked by 


the number of empirical equations used has 
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nates, 
coordinate as it is the property for which numerical 


sented only in graphical form. 


VI. Description of the Mollier Chart 


enthalpy 


in solving thermodynamic problems. 


chart such as first designed by Mollier [9] 


usually is 


values are of most interest 


chosen 


Five properties, namely, temperature, pressure, 
specific volume, enthalpy, and entropy are useful 


These prop- 


erties may be presented either in a table or on a 


Any 


two of the properties may be plotted as coordi- 


but 


ns one 


Any one of the other 
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four properties or a function thereof may be used 
for the other coordinate, depending upon the 
nature of the problems to be solved and the 
individual preference of the user. In the present 
instance, the use of the logarithm of pressure 
gives a constant percentage accuracy in reading 
the pressure and utilizes the rectangular paper 
efficiently. Usually all five properties are included, 
but for the sake of clarity, the constant volume 
lines have been omitted from the present Mollier 
chart. Values for the volume may be derived 
from figures 3, 4, and 5. 

The logarithms of pressures from 3 to 300 Ib/in.* 
have been plotted as ordinates Two scales 
along the abscissa have been used to represent 
enthalpies, the one at the left for the liquid, and 
the one at the right for the vapor. The curves 
representing the saturated vapor and the satu- 
rated to the 


left of the curve for saturated liquid represents 


liquid are so marked. The ares 
the subeooled liquid for which no data are avail- 
able. It should be noted that 
substances indicate that the constant temperature 


data for other 


lines are practically vertical over this area. The 
area between the curves for saturated liquid and 
for saturated capor represents mixtures of vapor 
and liquid. It has not been charted and the dis- 
tance between the two saturation lines has no 
physical significance. The area at the right rep- 
resents superheated vapor up to 300° F. 

The chart as originally drawn is 80 by 98 em 
(approximately 31‘ x 38% in The mid points 
of each coordinate line were located by means 
of a beam compass and an accurate seale, and 
the coordinate lines drawn in by means of a 
parallel straight edge. Various checks indicate 
that the coordinates are drawn with an accuracy 
of 0.22 or 0.3 mm, corresponding to about 0.05 
Btu/lIb and to about two parts in 1,000 of the 
pressure. 

Over 600 points were plotted so that each curve 
was drawn through a considerable number of 
suitably chosen points by means of a spline and 
weights. Errors in plotting larger than 0.2 or 
0.3 mm beeame evident due to the difficulty of 
fitting the spline to the points with a reasonable 
number of weights. This estimate of the accu- 
racy of plotting is supported by the fact that out 
of values read from 20 points on the original 
discrepancy between the 


drawing, the largest 


> 


chart and caleulated values was 0.05 Btu/lb (0.3 
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mm) and 0.0008 Btu/lb °R= in 


(0.3 mm). 


VII. Use of Tables and Charts 


The usefulness of the table and Mollier c, 
is based largely on three properties of the then 


entropy 


dynamic function enthalpy (sometimes called } 
content). The three properties referred to 
(1) In any process that eccurs at constant pre 
sure the increase in enthalpy is equal to the } 
in adiabatic compression— that is 
the work done by the e& 
pressor is equal to the change in enthalpy of | 


added; (2) 


constant entropy 


fluid between intake and discharge; and (3 
throttling process, such as occurs at an expats 
valve, the increase in enthalpy is equal to 
heat added. Therefore, when no heat is ad 
the enthalpy is constant. 

Two problems that illustrate the applicatior 
the three principles mentioned above, through 
of the chart, will now be discussed. 

Problem 1. 
tillation unit as liquid under an absolute press 
of about 80 Ib/in? and at a temperature of 70 
It is to be converted under the same pressur 
gas at 120° F. 


pound of butadiene? 


Butadiene is admitted into a 


How much heat must be added | 





Due to lack of experimental data, the chart d 
not give values of enthalpy for the subcool 


liquid; but from data on other substances it n 
be inferred that the change of enthalpy \ 
pressure for subcooled liquid butadiene at consta 
temperature is so small that for all practi 
purposes the value at 80 Ib/in.? absolute and 7 
F is the same as the value 236.5 Btu/lb for sat 


rated liquid at 70° F. The chart indicates th 


vapor leaving the still at 80 > Ib/in? abso! 
pressure and 120° F is practically saturated vap 
and that its enthalpy is 418.2 Btu/lb. 7! 


181.7 Btu per Ib mus: 
The answer sho 


difference 418.2—236.5 
be supplied to the butadiene. 
be affected only slightly by the presence of appre 
able amounts of other Cy hydrocarbons, as th 
have properties very similar to butadiene, but 
dissimilar compounds render 

The heat 


solution is involved in evaporation of mixtu 


presence of 


answer unreliable in such problems. 


that for dissimilar substances may differ appre 
ably from the heat of evaporation for butadic! 
Although correct answers can be obtained « 


when these data are applied to the pure mono! 


omkaie 
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of | 3-butadiene, yet in many of the problems to 
lved, mixtures of the monomer with polymer 
or with other compounds are encountered. In 
come cases if the composition of the mixture is 
known, a judicious use of the chart may lead to a 


sufficiently approximate answer. 

In the following problem, which involves such 
a mixture, the results of two methods of calcula- 
tion with the use of the chart are compared with 
the results of a third method of calculation in 
which data are used from tables or charts appro- 
priate for each component of the mixture. 

The correct answer to this problem could be 
obtained only from measurements made on the 
mixture itself. Such data are rarely available, 
but the results obtained by the third method are 
believed to approach the correct answer most 
closely. 

Problem 2. 
saturated with water vapor and containing about 
00 percent of butadiene and 8 percent of am- 
monia is encountered. The initial temperature 
and pressure of this mixture will be assumed to 
be 78° F and 25 |b/in.’, respectively. In actual 
practive the pressure of this mixture is increased 
only a few pounds per square inch through com- 
pression. We will assume for the present prob- 
lem, however, a compression ratio of two, i. e. a 
discharge pressure of 50 Ib/in.*, for the sake of 
sizable differences on the Mollier chart that may 
be read with sufficient accuracy for the compari- 
son of three methods of calculation. The prob- 
lem will consist in finding (1) the work performed 
by the compressor (2) the temperature of the dis- 
charged gas and (3) the increase in work of com- 
pression caused by pressure drop through valves 


In some processes a gaseous mixture 


or other constrictions. 

In the first method, the procedure is to caleu- 
late the results for pure butadiene as a rough 
approximation to the results for the mixture. 

The initial condition is represented on the 
Mollier chart by the intersection between the 
pressure ordinate 25 Ibn” and the constant 
temperature 78° F interpolated between the 
dotted curves for 70° F and for 80° F. This in- 
tersection corresponds to an enthalpy of 408.0 
Biu/lb and an entropy of 1.209. The intersection 
of this value for entropy with the pressure ordi- 
nate 50 Ib/in2 is at an enthalpy of 421.5 Btu/lb 
and a temperature of 120° F. Then for pure 
butadiene (1) the work of compression is 421.5— 
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408.0=13.5 Btu/lb and (2) the discharge tem- 
perature is 120° F, 

A second method, which gives a better estimate 
of the work of compression, uses the chart by tak- 
ing account of the fact that for equal compression 
ratios with the same inlet temperature, the work 
of compression per mole varies only 10 or 20 per- 
cent from substance to substance, whereas the 
work per unit mass varies over a wide range. On 
the molar basis the calculation is as follows: 


0.01665 Ib mole. 


0.00471 Ib mole. 
0.00111 Ib mole. 


0.9 lb butadiene 0.9/54 
0.08 Ib ammonia =0.08/17 
0.02 Ib water 0.02/18 


Total 0.02247 Ib mole. 


1 Ib butadiene= 1/54=0.01852 mole. 


Work of compression for 1 Ib of mixture is 


13.5 0.02247 /0.01852= 16.4 Btu/Ib. 

The third and probably the best method for 
estimating the work of compression is to sum up 
the work for compressing the various components 
separately through the same compression ratio 
from their initial partial pressures. Fortunately, 
such procedure is possible for the mixture presented 
in this problem, as tables of the preperties of am- 
monia [10] and of steam [11] are both available. 

The partial pressures are for butadiene 22.5 
Ib/in, for ammonia 2.0 Ib/in2,and for water the 
saturation pressure at 78° F, which is given by 
Keenan and Keyes in their table 1 as 0.4747 Ib/in.? 
Using the same procedure as before for butadiene 
at 25 |b/in? initial pressure, and the same initial 
temperature 78° F and compression ratio, two, 
one finds the same values for the work and the 
discharge temperature, for example, 13.5 Btu/lb 
and 120° F, respectively. 

The ammonia tables mentioned do not include 
pressures below 5 lb/in2 but the work of compres- 
sion for a compression ratio of 2 and initial tem- 
perature 78° F is nearly independent of initial 
pressure. For initial pressures of 5 and10 Ib/in’, 
the work of compression is respectively 46.8 and 
46.7 Btu/lb, the discharge temperature in’ both 
cases being 170.5° F. One may assume the work 
of compression to be 46.8 Btu/lb for an initial 
pressure of 2 lb/in2, and the discharge tempera- 
ture 170° F. 

The steam tables by Keenan and Keyes give 
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for the enthalpy and entropy of saturated steam 
at 78° F, respectively 1,095.8 Btu/lb and 2.0416 
Btu lb °R. The discharge pressure is 2 x 0.4747 

0.95 lb/in’. For convenience we will use the 
Ib/in 
temperature corresponding to an entropy of 2.0416 
and a pressure of 1 lb/in® are 1144.1 Btu/lb and 
186° F, The work is then 1144.1 

LO95.8 


rounded pressure | The enthalpy and 


respectively 
48.3 Btu/lb. 

The work of compression for the mixture is 
then, 


0.9 13.5 0.08 16.S 0.02 S35 


16.8 Btu lb, 
und the discharge temperature ts 


0.9 120-78 0.08 170—78 


125° F 


0.02 IS6—7S 7S 


We now have 
work of compression, namely 
13.5 Btu/lb 
using the chart for buta- 


three values calculated for the 


] For pure butadi he 
2) For the mixture 
diene together with the assumption that the work 
ol compression pet mole is the same for different 


16.4 Btu lb 


3) For the mixture using 


substances, 
the appropriate chart 
or table for each component of the mixture, 16.8 
Btu lb 
This last 


nearest to the correct answel 


method would be « X pec ted to give the 
but cannot alwavs 
ly used as data mav be lacking on some of the 
components. 

The discharge temperature 120° F) for pure 
butadiene is in this case lower than that (125° F 
for the mixture. This is due to the facet that the 
molar specific heats of ammonia and steam are 
less than that for butadiene 

The additional work that the compressor must 
perform in consequence of a pressure drop due to 
throttling through valves or other constrictions 
may be estimated by emploving the third property 
previously mentioned for enthalpy, 1. e. the con- 
ol enthalpy throttling process 


when no heat is added o 


staney aqduring a 


rt moved 


The loss is 
very nearly proportional to the pressure drop and 
may be caleulated by multiplying the work pet 


if on the Mother 


chart we follow from the point representing intake 


Ib in by the pressure drop 


downward along 


conditions (78° F and 25 |b/in 
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the constant enthalpy abscissa (408.0 Btu/lh) 4 
24 Ib/in.? the new entropy is 1.210. If. this 
constant entropy line is followed back to 25 Ib in 
an enthalpy of (409.0 Btu/lb) is obtained. Th 
additional work that the compressor must perforn 
in consequence of a 1 lb/in. pressure drop in th, 
suction line is 409—408=1 Btu/lb. By a simila 
process one finds that 1 lb/in? pressure drop oy 
the discharge side where the pressure is 50 Ib in 
causes the compressor to perform extra wort 
In the case where thy 
total compression amounts to only a few pounds 


amounting to 0.5 Btu lb. 


losses due to pressure drop may become an im- 


portant portion of the whole. 


Vill. Graphs of Certain Properties of 


1, 3-Butadiene 


A number of smaller graphs of the prope 
of butadiene are presented in this publication (1 


to furnish data not readable directly from tl 
Mollier chart and (2 
picture of the basis for, and relationships betweer 


Btu ll 


is especially in th 


to give the reader a clea 
these data. Figure 1 in which enthalpy 
is plotted vs temperature (°F 


latter class. 
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A. will be noted in this figure, the value for the 
enthalpy of the crystalline solid has been taken as 
nwro. This has been done for two reasons, 
namely, (1) exeept in problems involving chemical 
reactions no negative numbers are encountered 
with such choice (2) as many of the theoretical 
calculations of thermodynamics are based on the 
theorem that the entropy of crystalline substances 
is zero at the absolute zero of temperature, the 
zero for entropy was chosen for this state, and it 
seemed appropriate to assume the same state for 
the zero of enthalpy. 

This choice of zero state is satisfactory for all 
the problems for which this publication was 
designed. A value for the heat of formation from 
the elements [12] may be added to the enthalpies 
considered 

As illus- 
trated by the continuous curve in the lower left 
the enthalpy for solid 


ised here, but its assignment was 


beyond the scope of this publication. 


hand corner of figure 1, 
butadiene gradually increases with temperature 
ip to the triple point (—-164.05° F) then increases 
through fusion to liquid without any change in 
When the liquid is warmed, the 
enthalpy again changes with temperature at a 
rate that gradually The normal 
boiling point is reached at 24.06° F, in other words 


temperature. 
increases. 


the vapor pressure has increased to one standard 
With further warming of the liquid, 
the enthalpy continued to increase at a rate that 
temperature 


atmosphere. 
becomes infinite at the critical 
305.6° F, although the enthalpy itself is finite 
about 421 Btu/Ib). 

The heat of vaporization is zero at the critical 
temperature, but increases to finite values as the 
emperature is decreased (see dotted curve fig. 1). 


The rate of increase at the critical is infinite but 


At the triple point 
the heat of sublimation exceeds the heat of vapor- 


finite at lower temperatures. 


zation by an amount equal to the heat of fusion. 
\t lower temperatures the heat of sublimation 
hereases to a maximum and then decreases. As 
eas can be deduced from values of C,° from 
s formulation [8] together with the assumption 
hat ©C,° for the vapor approaches asymptoti- 


cally toward 4R at absolute zero, the enthalpy of 


‘vapor and consequently the heat of sublima- 
The 


thalpy of the saturated vapor is the sum of the 


m at absolute zero is about 287 Btu/Ib. 


thalpy for the solid and heat of sublimation, or 
higher temperatures the sum of the enthalpy for 


Thermodynamic Properties of 1, 3-Butadiene 


the liquid and the heat of vaporization. The 
enthalpy of saturated vapor increases’ with 
temperature without any break at the triple point, 
passes through a maximum at about 265° F and 
joins up with the enthalpy of the liquid at the 
critical temperature. The enthalpy of the vapor 
at zero pressure, which is the same as that at the 
standard state (ideal gas at 1 atmos), is not dis- 
tinguishable from the enthalpy of the saturated 
vapor at low temperatures, but at higher tem- 
peratures the difference calculated from the equa- 
tion of state increases so that the enthalpy at 
p= increases steadily with temperature without 
any peculiar turn at the critical temperature as 
exhibited for the saturated vapor. 
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2 Svecific heat at constant pressure of qaseous 


1, 3-hutadiene. 


FiGguri 


Figure 2 is a plot of C, the specific heat of the 
vapor at various pressures from 0 absolute to 
200 Ib in abs as a function of temperature ex- 
pressed in degrees F. The values of C, at zero 
pressure are ealeulable from spectroscopic data 
but were actually calculated from an empirical 
equation [8] and the values at other pressures can 
be deduced from these by using the equation of 
state. The temperatures at the bottom of the 
graph apply to the lower set of curves, those at the 
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top to the upper set, with the two sets overlapping 
in the range 120° to 140° F. The curve for 200 
lb in.? shows a minimum in C, for the vapor, and 
if the curves were computed for higher pressures, 
these minimums would be more pronounced. 
Figure 3 shows values of the quantity Pe/RT 
which is unity for an ideal gas) as a function of 
pressure for saturated vapor and for superheated 
vapor at various temperatures. It may be noted 
that at the highest pressure plotted, 120 |bs/in® 
the value of Pe for saturated vapor is more than 
16 percent less than the value for an ideal gas, but 
that at low pressure, the values of Pr approach 
RT, 


temperatures, 


the approach being closer for the higher 


Pe/RT 


Figure 4 isa plot of the quantity Pike a 


10° 





Py/RT 
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“oO © 2 © 4 SO 6 7 8 9% WO WM 120 
PRESSURE. LB/IN* 

hi Ri 3 Pressure-to ime-lemp atu 


which is equal to zero for an ideal gas) for the 
vapor, as a function of the absolute pressure for 
various temperatures from — 40° to 180° F. The 


quantity plotted never reaches the value for an 


ideal gas, but approaches it as the pressure is 


decreased or as the temperature is raised 
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Ficure 4 Pressure-volume-te mperature relations fo 


gaseous 1, 3-butadiene 
For use at higher pressures and temperatures 
the coordinates used in figure 4 require too large 


a page for accurate presentation of the data 


This difficulty is overcome in figure 5 by plotting | 


) as ordinate vs pres- 


the product 1 PuRT (_T 
ie produc P Ib in. ( 100 


sure as abscissa. Through the use of these coor- 
linates accurate values of Pe/RT may be obtained 
over the whole pressure and temperature rang: 
although 


more than 


calculation is involved 
using figure 4 
well as the saturation curve at low pressures an 
temperatures. In this region the isotherms hav 
been omitted and the use of figure 4 will be less 


ambiguous. 


An acknowledgment with thanks is given 
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appearing in the table and to Robert E. MeCosk: 
for calculating the data needed for the Molli 
chart and for drawing this chart. 
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Heats, Equilibrium Constants, and Free Energies 
of Formation of the Alkylcyclopentanes and 
Alkylcyclohexanes ° 


By John E. Kilpatrick,” Helene G. Werner,’ Charles W. Beckett,* Kenneth S. Pitzer,‘ and 
Frederick D. Rossini ° 


For ecyclopentane, cyclohexane, the seven dimethyleyclohexanes, and the normal 
alkvlevelopentanes and alkyleyclohexanes, values are presented for the following thermo- 
dynamic properties to 1,000° or 1,500°K: The heat of formation from the elements; the 
free energy of formation from the elements; and the logarithm of the equilibrium constant of 
formation from the elements. For cyclopentane, cyclohexane, and the normal alkvicyclo- 
pentanes and alkyleyclohexanes, values are also given to 1,500°K, for the following proper- 
ties: The heat-content function; the free-energy function; the entropy; the heat content; 
and the heat capacity. 

Equilibrium constants and concentrations are given in tabular and graphical form for 


some reactions of isomerization, hydrogenation, and cyclization. 


I. Introduction stants of a number of reactions involving these 
compounds have been calculated. 
{s a part of the work of the American Petroleum 
Institute Research Project 44 at the National II. Fundamental Constants 


Bureau of Standards and the University of Cali- The unit of energy, atomic weights, and values 
fornia, values have been calculated for the heat- of the fundamental constants used in this report 
content function, free-energy function, entropy, are the same as those previously used [5, 13]°. 
heat content and heat capacity of the normal alkyl- 


eyclopentanes and the normal alkyleyclohexanes III. Heat-Content Function, Free-Energy 
above ethyleyclohexane, in the gaseous state to Function, Entropy, Heat Content, and 
1.500° K. These data, together with previously Heat Capacity 

published data on cyclopentane, cyclohexane, 
methyleyclohexane, and ethyleyclohexane, have 1. Method and Data Used in the Calculations 


been combined with values of the heats of forma- For cyclopentane, cyclohexane, methvlevelo- 


4 hexane, ethvleyvclohexane, and the seven dimethyl- 


tion at 25° C to calculate values of the heats, free 
nergies, and equilibrium constants of formation cyclohexanes, the values were taken from refer- 
17 alkvleyclopentanes and 24 alkvleyelohex- ences [1, 2]. 


The free energies and equilibrium con- For methyleyclopentane, the thermodynamic 


snd Inorganic Chemistry of functions were calculated by means of the follow- 


Y., Sept. 1947 ing relation: 
the work of the . 


nal Bureau [methyleyelopentane] = [cyclopentane] + [methyl] 
+ [restricted rotation of ring puckering]+ (2? In 10). 

(1) 

The terms written in brackets represent a given 


thermodynamic function of the indieated sub- 
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stance. In ea 1, the term [methyl] represents the 
contribution (vibration and restricted rotation 
of the methyl group and was calculated by the 


following relation: 


methyl] — [methyvleyclohexane]— [steric factor] 


Rin 6 2 


In eq ae the term |steric factor] is the contribution 


evclohexane] 


of the equatorial-polar tautomerism of methyl- 


cyclohexane to Its thermodynamic functions 
The term 


tween the functions of methylevelohexane and 


methyl] is therefore the difference be- 


evclohexane, without the inclusion of the effects of 


tautomerism and of symmetry number. The 
svmmetry number corrections, which are written 
in parentheses, are to be used for the negative of 


function, F—T1,)/T, and for 


. but not for the other functions. 


the free energy 
the entropy, S 

The term [restricted rotation of ring puckering| 
is the contribution arising from a restricting bar- 
rier of 750 cal/mole from the asymmetry induced 
by the presence of the methyl group on the ring 
See reference [1] for details A classical partition 
function was used 

The thermodynamic functions of ethylevelo- 
calculated from the following 


pentane were 


relation 
ethyvlevclopentane] = [methyleyclopentane] 
2-methyl butane] —[2-methvlpropane] 


Rin 3 >) 


For normal props leyclopentane and each of the 
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higher normal alkyleyclopentanes, the follow ing 


relation, with the appropriate number of CH 
groups, was used to calculate the thermodyna ni 
functions: 


fnormal alkvleyelopentane] lethvleyvelopenta ( 
n (CHL. 1 


The thermodynamic functions of normal propy 


evclohexane were calculated from the following 


relat Ion: 


jnormal propylevclohexane] = [ethvlevclohexane 


[2-methyvlpentane]— [2-methylbutane]. 5 


In the case of normal butyleyclohexane and thy 
higher normal alkvlevclohexanes, a relation simila; 


to eq 4 was used: 


fnormal alkvleyclohexane] = [normal propyl- 


nN [(CH,]. 0 


cvclohexane| 


The thermodynamic functions of the parafli: 
hvdrocarbons and of the CH, increment that ar 
used in eq 3 to 6 are taken from reference [14] 


2. Results 


The resulting values of the thermodynamy 
properties for the normal alkyleyclopentanes and 
alkvleyclohexanes are presented in tables 1 to 10 
which give values of the following properties 
Heat-content function, (77 Hy) i: 
funetion, (F°—//,)/7T; entropy, S°; heat con- 


tent, // /7,; and heat capacity, C, 


free-energy 
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Heat 





BLE | Values of the heat content function, for the ideal gaseous state, to 1,500° K, for the normal alkylcyclopentanes 











\ ’ 
‘ i remperature in °K 
ann 
Compound For ( 208.16 00 4100 Mn) 600 700 SOO) O00 1.000 1,100 1,200 1.300 1,400 1 
ga ila 
fT] Heat content function (H°-H lr, in cal/deg mol 
} 
wntane Cs5H ) 12. 07 12. 12 15.12 18.52) 21.98) 25.29 28.40 31. 29 $3. O7 6. 46 a8. 75 40). 87 42.83 44. 65 
cyclopentane C,H 0 16.01 16. 09 19.88 24.03) 28.17 32.09 5. 75 wy. 15 42.28) 45.18) 47.87) 50.34) 52.62) 54.74 
| vclopentane CrHy 0 18. 91 19. 01 23.70 28.66) 33.55 38.16 42.43 46. 38 00) 53.36 56.47 50.32) 61.96 64.40 
“ vicyelopentane CsHy 0 | 22.34 | 22. 4€ 27.9 33.6 39. 1 44.4 49.2 53.7 57.8 61.6 65. 1 68.3 71.3 74.1 
cyclopentane Col 0 | 25.77 | 25.90 32.1 38.5 44.7 “6 6.0 61.0 65. € 69.8 73.7 77.3 80.6 &3.7 
vieyclopentane CwH 0 | 29.20) 29.35 6.3 43.4 3 56.8 62.8 68.3 73.3 78.0 82.3 86.3 90.0 93.4 
Icyclopentane C),Hy 0 $2. 43 32.80 40.5 48.3 55.9 63.0 #9. 5 75.6 Si.1 86.2 91.0 95.3 90.4 103.0 
i¢ vicyclopentane ( HH. 0 16. 06 6. 24 44.8 53.2 61.5 69.2 76.3 S28 ako 04.4 oo 6 14.3 108. 7 112.7 
cyclopentane ( H 0 su. 4Y 34. 69 49.0 58. 2 67.1 75.4 SS. 1 vO. 1 06. 6 102. 6 108. 2 113.3 118.0 122. 4 
yelopentan C yHo 0 42.92 43.14 3.2 63.1 72.7 R1L¢ 80.9 97.4 104.4 110.8 | 116.8 | 122.3 | 127.4 132.0 
cyclopentan CH 0 46.35 46.59 57.4 68.0 78. 2 87.9 WH. ¢ 104. 7 112.2 | 119.1 125. 4 131.3 | 136.7 141.7 
I cyleyclopentane CH 0 | 49.78 | 50.03 61.6 72.9 83.8 94.1 | 103.4 | 112.0 | 120.0 | 127.3 | 134.0 | 140.3) 146.1 | 151 
ri evicyclopentane ( HH ) 3. 21 53. 48 65.8 77.8 20.4 100.3 110.2 119.3 127.7 13. 142.7 149.3 155. 4 161.0 
vleyclopentane Call 0 | 4.64 |) 56.93 70.0 a2.8 05.0 | 106 117.0 | 126.6) 135.5 | 143.7) 1 3) «158.3 | 164.8 | 170.7 
idecvleyvclopentane CH 0 60.07 60.37 2 87.7 100.6 112.7 123.8 133.9 143. 151.9 159.9 167.3 174.1 180.5 
idecyleyclopentan CoH, 0 | 63.50 | 63.82 78.4 92.6 | 106.2 | 118.9 | 130.5 | 141.2 | 151.0 | 160.1 | 168 176.3 | 183.4 | 190.0 
leevievel ntam ( H 0 66. 93 67. 27 82.6 97.5 111.8 125.1 137.3 148.5 ISS 8 168.3 177.2 185.3 192.8 100. ¢ 
t per CH, group 0 ,. 430 $447 4. 21 4.92 ». 59 6. 21 6.78 7.29 rears 8. 2) &. 62 00 9, 35 0. Ot 
{ 
if] 
ABLE 2 Values of the heat content function, for the ideal gaseous state, to 1,500° K, for the normal alkylcyclohexanes 
remperature in °K 
‘ } il ' Jus. 16 uM 4M LD 600 Too siM) GOO 1.000 1. lx 1, M0) 1, 300 1, 400) uM 
i 
I 
Heat content function (/7°-H T, in cal/deg 
tithe 
iv " CoH 4.2 1.28 8.38 22.8 7. 34 1. 4 67 39.40 42.8 45. 5.9 t 4.0 “ 
ies hexa C°H 7. tit 22.84 28.2 s s 18. It 17.42 4.9 8.2 61.2 4. (it 
hexa Cu o Ww4 ws 2.9 so 44 19.9 1.7 ’ t iis 70.2 7 76.2 
rey lohexa CoH 0 23.66 23.8 0. ¢ 7.7 14 wt 4 61.8 6.7 75.4 79, 2 82.6 85.9 
on hexane CwH ( 27.00 27. 2 4.8 $2.6 w.0 iS 4&3. 2 69.1 74 7¥ S40 SA. 2 92.0 oe) 
yclohexar C,H 0. 52 0. 70 0 ‘7 63.0 70.0 76.4 82.2 87.7 2.7 7.2 101 105, 2 
hexa CoH o 11 43.2 2 t iv. 2 6.7 83. 7 w{ "0 0 2 10.7 14.8 
hexa Cyl s “i 7.4 7.4 ti. ¢ 7 s 1.0 17.8 4.2 OY 2 20.4 a4 
hex CH 0 4.8 11.0 é 2 22) 81.7 ”) s O05. 24 118 244.2 120.4 134.2 
hexa CisH ( 4.24 44.49 ’ 67.2 77.8 87.9 7.1 05. ¢ Ww. 27.1 2 is 43.8 
hexa CoH 7.¢ 47.04 ti 72 S34 "4 3.8 29 2 5 ‘ 42.2 is 
‘ hexane Cc ' 4 77.0 NY. ( T) vf a mY 0 4.4 2 { 1th 
cyclohexane CisH 4 1.8 is 82.0 Mt in 7.4 27.4 Yi 45.2 0 2 ‘6.8 172.8 
c hexa Cw ” 8. 28 2.7 86.9 00. 2 2.7 24.2 4.7 404 $ I6le iv. 2 7 Is 
yl hex CoH ( ‘ >| 61.7 6.9 18 05.8 54 0.9 42 2.2 61.7 70.2 8.2 x 2 
adlecy hexane ( Hi ( 4 SZ th s mM. 7 125. 2 7.4 tt] “O.0 “4.4 73.4 87.2 4.8 Dis 
cyl exa CoH is 2 S. 62 s 11. ¢ 16.9 4 44 1M. ¢ 67.7 7s S7 196.2 24.2 2 ‘ 
per CH» grou 4 7 12 4. 92 ” 6.2 6. 78 20 7 8.2 8. (2 00 ’ 1, i 
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, PaABLe 5 Values of the en opy, for the deal gaseous slate, lo 1,500 K, fo ti normal all yic jclo pe nianes 
Temperature in © K 
= . = * und I Js. If wn iu vn un 70 sin) PLD on in 1, aK 1, 1, 40 ‘Ti 
Entre py, ! " te t 
Zz 
ntam C.H 0 70.00 70. 12 77.00 84.14 G1.28 98.24 104.92 111.33 17.45 123.30 128.88 134.18 139.23 144.07 
cyclopentane CoH 0 81.24 81.42 9.33 90.36 108.24 116.81 125.01 132.83 140, 25 172. 37 
= lopentane CH 0 ”). 62 GO. 84 101.60 112.38 122.89 133.05 142.73 151.91 160.62 198, 19 
aA cyvclopentane CH ' GSO 100. Of 112.7 125. 2 137.3 149.0 10.1 170.7 1s. 6 293. ¢ 
. cyclopentane ColHyse ( Os. 09 109. 20 123.8 138.0 151.8 165.0 177.5 180.4 An), 7 249. { 
icyclopentane CwH { 18.17 118. 51 134.8 150.8 166. 2 181.0 194.9 206.2 2007 274 
yelopentane Cy @ 127.35 | 127.74 145.9 163.6 180.6) 196.9 212.4 226.9 240.7 . 
3 vieyclopentane CoH 0 i. 54 136. 96 157.0 176.4 195.0 212.9 220.8 | 245.7 | 200.8 {25.3 
" yelopentane Cyl { 45.72 146.18 6s. 1 189. 2 228.8 | 247.2 | 264.4 | 280.8 3h). 8 
cvclopentane ( Hos ( 4. oO 155. 41 2.0 244.8 24. 6 2R3. 2 Lim. $76.2 
-* eyclopentane ( Hi 0 64.08 164. 63 2 214.8 260.8 282.0 “1.9 120.9 lf 
t= ¢ cevieyclopentane Cy 0 173.27 173.86 3 | 227. ¢ 76.7 209.4 320.7 340.9 360.0 78.2 | 305.4) 411.7 | 427.0 
i vieyclopentane C,H ( R245 183, OR 4 | 240.4 | 287.2 | 202.7 316.8 330.4 0.9 BR1.3 400.¢ 418.9 436.2 452 
vicyclopentane CH ( 91.63 192.30 223.4 253.2 281 ¢ 308.7 334.2 358.2 381.0 402 123.0 442.4 400.6 477.9 
» = lecvicyclopentane Cll 0 200.82 WIS 234.5 26.0 206.0 $24. 6 S51. 6 77.0 wo 423. 7 $45.3 465.8 485. 1 M3. 3 
= udecyvieyclopentane CoH 0 210.00 210.75 245.6 278.8 310.5 +40), € 169.0 105.7 $21.0 $45.0 467.7 $80.3 “Nn. ¢ 5oR. 7 
“ lecvicyclopentane CoH 0 219.18 | 219.98 256.7 291. ¢ $24.9 $5 ahh. 4 $14.4 441.0 466.2 490.1 12.7 4.1 4. 2 
t per CH), grouy { ¥. 183 9. 224 11.08 12. Ss) 14.43 7.40 S.7 20. 0 21. 24 22. 38 23. 4 24.48 25. 43 
PARLE 6 Values of the ent opy, fo the ideal gascous state, to, 1500 K, for the normal alk sleyclohexane 8 
" Temperature in ° K 
- 7 . “ . : f Js. If SLU wi “nH man 70) si) wn ooo 1, 100 AK) 1, 300 1, 40 wm 
--* 
Entropy, S n cal deg ‘ 
-% 
vane CH 28 71.44 SO.IN 80.24 W300 107.14 115. ¢ 23. SO 0 139.0 4.0 152.7) 159.0 165.1 
t C;H { S2. (Mt 82. 28 %.22 104.32 21 125.77 85. Be 5.46 154.57 163.2 71.4 1702 186. 193. 7 
! CH 0 1.44 91.70 104 117.3 129.9 142.0 153.¢ 14 74.9 | 184.8 | 194.2 | 203.0 211.4 219 
'~ 2 cyclohexane CoHys 0 WwO1l 100. 43 115.0 129.6 143.7 157.3 170, 4 182.7 194.3 205, 5 215.9 226.0 235.3 244.4 
* yelohexane Croll 0 109.29 109.66 126. 1 142.4 158.1 173.3 IS7.8 201.5 214.3 226.7 238.3 249.5 250.8 269.8 
cyclohexane CyH { 8.48 118.87 7.2 2 | 172. 89.2 | 205.2 | 220.2 | 234.4 248.0 260 72.9 | 284.3 295.3 
evclohexane { iH ( 27. 128. 10 148.2 148.0 187.0 205. 2 222 ¢ 239. 0 254.4 249. 2 283.0 206 4 wR. 7 i200. 7 
~ yelohexane Cyl 0 136.84 137.33 150.3) 180.8 201.4 | 221.2 240.0 0 257.7 274.4 | 200.4 305.4 319.8 | 333.2 | 346.1 
= yclohexane CyH ) 146.038 146, of 170.4 U3. ¢ 215.8 | 237.1 257.4 276 244 11.7 327.8 343.3 | 357.7 7 
1 cyclohexane Coll » 1 21 5. 77 Isl.5 206.4 230 253 274.8 | 295.2 $14.5 $32.9 $0.2 366.8 382.2 397.0 
yelohexane ( H ( (4. 3Y 165. 00 1v2 219. 2 244.7 2. 1 202. 2 14.0 34.5 4. 2 372.6 ww. 2 106. 7 422.4 
1 cyleyclohexane H 0 8. ST 74.23 203.¢ 232.0 259.1 285.0 sO. f $2.7 34 375.4 304.9 $13.7 451.1 447.5 
= vicyclohexane CywH i S2. 7f ISS. 44 214.7 244.5 273. ¢ 1.0 27.0 s51 74.6 sh. f 417.3 437.1 $55. ¢ 473.2 
1 vieyclohexane CywH 1.94 192.67 225.8 257.6 QRS. ( 7.0 344.4 70.2 4. 117.9 439.7 460. is).1 498.7 
lecyicyclohexane CoH oO 2 2 Die 236.8 270.4 wr. 4 32.0 61.8 INU. 414. ¢ $39. 1 Wiz is4. 1 i. € 524. 1 
adie j lecyleyclohexane CyH 210.31 211.11 247.9 283.2 6.9 348.9 3702 40 134.7 460.3 484.4 507.5 520.1 549 
} evi hexa CeH 219.49 220.34 259 204 ‘4.9 M.7 426 4.7 isi. f iS 1.0 ; Th.0 
I CH i 41S 224 s 2M 14 ; i s 2.03 21.24 22.38 23.4 24.48 25.43 
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TaBLe 7.— Values of the heat content, for the ideal gaseous state, to 1,500° K, for the normal alkyleyclopentanes 














T 
Temperature in ° K 
Compo ind For- ; > | } : a | een 
(gas) mula |? | 298.16 300 400 500 600 700 SOO 900 1,000 1,100 1,200 1 1,300) 1,400 1, 
Heat content (/7°—H°), in cal/mole 

C yelopentane CsHy 0 3,509 3,6 6.048 9,280 13,188 17,703 22,720 28,161 33,970 40,110 46,500 53,130 59,960 & 
Methyleyclopentane CoHy 0 4,774 | 4,827 7,952 12,015 16,902) 22,463 28,600 35,235 42,280 49,700) 57,440 65,440 73.670 82 Byclopents 
Ethyleyclopentane CrHy 0 5, 638 5, 708 9, 480 14,330 20,100 26,700 33,900, 41,700 50,000 58,700 67,800) 77,100 86,700 & BR erhvieve 
n-Propyleyclopentane Coy 0 6,661 6, 737 11,160 16,790 23,500 31, 100 39,400 48,300 57,800 67,700 78,100 88,800) 99,800 11) BBrhylcyele 
n-Butyleyclopentane Cols O| 7.683 | 7,77 12.850 19,250 26,800 35,400 44,800 54,900 65,500 76,800) 88, 400 100, 500 112, 900 125 BR propyley 
n-Pentyleyclopentane CyHy 0 8,706 | 8,805 14,530 21,710 30,200 39,800) 50,200) 61,400 75,300 85,800 98, 800 112, 200 126, 000 140 Bb Butyley« 
n-Hexyleyclopentane CnuHa 0 9,729 9,839 16,220 24,160 33,500 44,100 55, 600 68,000) 81, 100 94, 800 109, 100 123, 900 139, 100 15: Pentyley 
n-Heptyleyclopentane Col, 010,752 10,874 17,900) 26, 620 36,900 48,500 61,100 74,600 88, 900 108, 800 119, 500 135, 600 152, 200 1 BB Hexyleye 
n-Oct yleyclopentane CyHe OUL774 (11.908 19,580 29,080 40,200 52,800 66,500 81, 100) 96, 600 112, 900 129, 800 147, 300 165, 300 Is. Be Heptyles 
n-Nonyleyclopentane C Hos 012,797) 12,942 21, 270) 31, 40) 43, 600) 57, 200) 71,900) 87, 700 104, 400 121, 900 140, 100 159, 000 178, 400 1% Bb Oct yleye 
n-Decyleyclopentane CisHye 013,820 13,976 22,950 34.000 46.900 61,500 77,300 94, 300 112, 200 130, 900 150, 500 170, 700 191, 500 212 He Nonyley 
n-Undecyleyclopentane CyHe 014,842 15,010 24,40) 36,460) 50, 300) 65,900) 82, 700 100, 800 120, 000 140, 000 160, 800 182, 400 204, 600 22 BB Decyley« 
n-Dodecyleyclopentane Cry 015,865 (16,044 = 26,320 38,920 53.600 70,200) 88, 200 107, 400 127, 700 149, 000 171, 200 194, 100 217, 600 24) HU ndecyl 
n-Tridecylevclopentane CrsHas 0 16,888 17,078 28,000 41,370) 57,000) 74,600) 93, 600 113, 900 135, 500 158, 000 181, 500,205, 800 230, 700 2% Dodecyh 
n-Tetradecyleyclopentane CigHys O17,910 (18,112 29,690 43,830 60,400 7S, 900) 99, 000 120, 500 143, 300 167, 000 191, 800 217, 500 243, 800 270 Tridecy] 
n-Pentadecyleyclopentane CoyHy 018,933 19, 146 31,370) 46, 280) 63,700) 83, 200 104, 400 127, 100 151, 000 176, 100 202, 200 229, 200 256, 900 2k HD Tetradec 
n-Hexadecyleyclopentane CyHe 019,956 20180 33,060 48.750) 67, 100) 87, 600 109, 900 133, 600 158, 800 185, 100 212. 500 240, 900 270, 000 2 HD Pentades 
Increment per CH, group 0 1,022.7 1,084.1 1,68) 2.458 3,350 4.350 5.420 6.560 7.770 9.080 10,340 11,700 13,000 1) fi Hexadec 
crement 
TABLE 8.— Values of the heat content, for the ideal gaseous state, to 1,500° K, for the normal alkyleycloheranes I 

Temperature in ° K 
Compound For- | » | 28.16] 300 0 = 500 s| 800 | | «1,000 «1,100 1,200 1,300 | 1, 400 
vas mula | 
Heat content (/7°-—-H°%>), in cal/mole 
C yelohexane Cy 4, 284 7.352 11,435 16, 404 28, 536 460 42,850 50,600 58,680 67,080 75,609 9 i yelohexs 
Methyleyclohexane Coll 5, 208 9136 14,130 DL 1IY 26.971 34,528 42,678 51,330 60.300 60.840 70,500 SO GO & lethyley 
Ethylevelohexane CH, 6,174 10, 680 16,450 23,300 31,200 39,900 49,200 59,100, 69,400 80,200 91, 300 102, 600 1 AP hylevel 
n-Propyleyelohexane Cols 7,148 12,240) 18,850) 26,600 35,400 45,100 55,600) 66,700 78,400 90, 500 103,000 115, 700 12 Propyle 
n-Butyleyclohevane CH» A177) 13,920) 21, 310 2,200 74, 500) 87,500 100, 800 114, 700 128, 800 14 Butyley 
r-Pentyleyclohexane Cyl 9, 211 15,610 23, 770 38, 700) 82,200) G6, 500-111, 200 126, 400 141, G00 E> Pentyle 
n-Hexyleyelohexane Coll 19,245 17,280) Oe Be 5, 300) 0, 000.105, 500 121, 500 138, 100 155, 000 17, Hexyley 
n-Heptyleyclohe xan Cooly 11,279 = 18, SO, wn!) 52,800 66,800 81,900 97,800 114, 600 131, 900 149, 800 168, 100 ISN Heptyle 
n-Oct yleyvelohexine CyHs 12,313 20,660 31,140 43,300) 57,200) 72,200) 88, 400 105, 600 123, 600 112, 200 161, 500 181, 200 2 Ocetyley 
n-Nonyleyelohexane CoH 348 22.340 33,600 46,700 61,500 77,700 95,000 113, 300 152, 600 152, 500 173, 200 194, 200 21 y-Nonyle: 
n-Deeyleyclohexane Cully 14, 382 24,050) 36,060) 50,0005 65, 900) 83, 100 101, 600 121, 100 141, 600 162, G00 184, G00 207, 300 29 Deevley 
n-U ndeevleyclohexane Cyl 15,416 25, 710 38,520) 53,400 70,200 88, 500 108, 100 128, 900 150, 700 173, 200 196, 600 220, 400 24 Undeey 
n-Dodeeyleyclohexane CH» 16,450 27,400 40,980 56,700) 74,600) 43, G00 114, 700 136, 700 159, 700 183, 600 208, 300 233, 500 2° Dodecy 
n-Trideeyleyclohexane CoH 282 17,484 20,080) 43,440 60,100 78,900 99, 300 121, 300 144, 400.168, 700. 198, 900 220, 000 246, 600 27 Trideey 
n-Tetradecyleyclohexane CoH O 18,304 18, 518 3), 7H) 45,890 68, 500) 83, 300 104, $00 127, 800 152, 200 177, 800 204, 200 231, 700 250, TAO 2. Tetrade 
n-Ventadeeyleyclohevane Calle 019,327 19, 552 32.450) 48,3500 66.8000 87,600 110, 200 134, 400 160,000 186, 800 214, 600 243, 100 272, S00 we Pentad 
n-Hexadeeylevclohe vin Colt 020,350 2.586 34,130 50,810) 70, 200 99,000 115, 400 141,000 167, 700. 195, 800 224, G00 255, 100 285, goo ay: MP iexade 
Increment per CH, croup 0 1,022.7 1.0341 168 2458 3.35) 4.350 5.420 6.560 7.770 9.030 10,340 11,700 13,090 1 glneremen 
For 
For 
_ 

rhe 
deseri 
For 
eveloy 
in this 
reu 
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TABLE 9. Values of the heat capacity, for the ideal gaseous state, to 1,500 
Temperature in ° K 
a ~ = 0 | 208. 16 300) 400 0) 600 700 800 900 | 1,000 | 1,100 | 1,200 | 1,300 1,400) 1, 500 
Heat capacity, C,, in cal/dec mole 
entane CH 0 19.82 19. US 2k, 24 35. 86 42. 36 47.81 52. 44 56. 37 0.75 62.68 65, 18 67. 36 69, 24 70. 89 
cyclopentane CoH 0 26.24 26. 46 $6. 11 44.94 52.43 58.68 64.00 68.53 72.44 75.82 78.72 81.24 85. 35 
yclopentane C;H 0 31.93 32. 18 43.39 53.55 62.00 69.24 75.31 80.48 84.94) SSSI 92.12 95.00 99. 69 
Propyleyclopentane CH 0 37.69 37. 99 50. 5 62.0 71.6 79.6 86.4 92.3 97.3 | 101.7 105.4) 108.7 113.9 
jutyleyclopentane Coys 0 43.45 43. 80 7.7 70.4 81.0 90.0 97.6) 104.1 | 109.7) 114.5 118.7) 122.3 128. 2 
ntyleyeclopentane CyH 0 49.21 49. 61 64.8 78.9 100.3 108.7) 115.9 122.0) 127.4 132.0 136.0 142.4 
Hea vieyclopentane CyH 0 54.97 5. 42 71.9 87.3 110.7 119.9 27.6 | 134.4 140.2 145.3 149.7 1M. 7 
feptyleyclopentane CyHy 0 &.73 61. 23 79. 1 95.7 121.1 131.0) 139.4 146.8 153.1 158.6 163.4 171.0 
vtvleyciopentane CyH 0 66.49 67.04 86.2 104.2 131.4 142.2) 151.2 | 150.2 | 166.0 171.9 177.0 185, 2 
Nonyvleyelopentane Cy Hos O 72.25 72. 85 93.3 112.6 141.8 153.3 163.0 | 171.5 | 178.8 185.2 190.7 199. 5 
ecvieyclopentane CiysH » 78.01 78.66 100.4 121.0 152.2. 164.4 | 174.8 | 183.9 | 191.7 198.4 204.4 213.7 
ndecyleyclopentane CywH 0 83.77 84.47 107.6 129.5 162.5 175.6 186.6 196.3 204.6 211.7 218.0 228. 0 
odecyleyclopentane Cy 0 89.53 90.28 114.7 137.9 172.9 186.7 | 198.4 | 208.6 | 217.4 225.0 231.7 242.2 
lecyleyclopentane CH 0 95.29 96.09 121.8 146.4 183.3 197.8 210 221.0 | 230.3 238.3 245.4 Lith 
radecyleyclopentane CywH 0 101.05 101.90 120.0 154.8 193.6 209.0 | 222.0 | 233.4 243.1 251.6 259.0 270. 7 
Pentadecyleyclopentane CyeH 0 106.81 107.71 136. 1 163. 2 24.0 220.1 233.8 | 245.8 256.0 264.9 272.7 ‘4 285.0 
lexadecyleyclopentane CoH 0 112. 57 113.52 143.2 171.7 214.4 231.3 | 245.6 | 258.1 268.8 278.2 | 286.4 293.3 | 200.2 
ent per CH» group 0 5.760 5.810 13 844 10.37) 11.14 11.79) 12.37) 12.86 13.29) 13.67) 13.99 14.25 
TABLE 10. Values of the heat capacity, for the ideal gaseous state to 1.500 K, for the normal alkylcycloheranes 
Temperature in ° K 
‘ na gs ba. 0 Us. 16 400) 100 OO) 600 TOO S00 900 1.000 1.100 20) 1.300 1,400 1,500 
Heat capacity, ¢ in cal/deg mole 
exane C.l Tt) 2A. 40 25. 5s s2 45.47 ‘Se WO. ST i. Tt 71. 68 75. 80 7Y¥ a2. 2 4.7 SHS AA. f 
vieyclohexane C:H ( 2.27 $2. 51 44.35 55.21 4.46 72.23 78.74 §4.20 SS.79 «92.7 Ow. ( MS 101.2 108.2 
velohexam Cl 7. 4) 51. i. 74.1 S28 wi w.2 101 105. 7 1. 4 112. 4 If 117. 
vievelohexane Coll 0 4 2] 45.80 ms. f 72.1 be] 03.1 11.1 107.9 113.6 118. ¢ 122.7 12.2 120.34 131.9 
itvleyelohexane CoH ( 9.3 19.70 65.7 8) 5 @.0 108 12.2) 119 146.0 131.4) 126.0) 139.8 143.3) M61 
tvlevclohexane ( iH ( 5A. 11 58. 5) 12.9 SY. 0) 12. 4 113.8 123.4 131 13S i144 14yu ; Lav 4 
Hexvleyelohexane Cyl 0 6O.87 61.32 S00 v7.4 111.9 | 124.2 | 134.5 1458 1M) 157.1 12. ¢ 167.2 71.2 174.6 
Heptvleyclohexane ‘ i tii, 67.1 S71 105, 8 121.4 134.6 145.7 155.1 i. 1 170.0 175.8 ISOS 185. 2 ISS. Y 
yieyelohexane CoH 0 72.29 4 4.3 #1143 DOS 144.9) ISS) 16.9 75.4 | 182.9 | 189.1 144 2.1 
nvlevelohexane Coll 0 7s 8.75) (109.4 122.7 | 140.3 | 155.3) 167.9 | 1787 | 187.8 6.7 | 202.4 | Aw.2 217.4 
wevieyelohexan Cyl x M6 Os 31.2) 1408 165.7 179.1 wv) 20. 2 iS. 215.7 | 221.9 227.2 | 231.6 
lecyileyclohexane Coll SY uF WO. 37 115.7 139. 6 1.2 «176.0 19.2 22.2 |) 212.6 21.4 | 220.4 2 241.2 M59 
lecvleyclohexane Cull 0 5.4 68. 1s 122. 8 148.0 168. 7 Ist. 4 i.4 214.0 224.9 24 242 244.2 255. 2 a0 1 
leevleyclohexane ( H © W119 1. wy 4.4 lit 178.2 1H. Ss 212 247.2 a, # 729 2.2 74.4 
decyleyclohexane ( Hi 0 106.9 107. SO 137.1 14.9 187. € D7. 1 273. ¢ 2.0 aU vt 283. 1 2s. ( 
tadecyleyclohexane Coll u 12.7 113. 61 144.2 17 197.1 217.5 24.8 272 9 | 262.2 w.2 27.1 wr2. Y 
vadecyviey clohexane ( iH 0 118.47 11%. 4z Sl ISLS WME 227.9 245.0 221.2 2744 m5. 7 205, 5 tM ’ 11 17.1 
ent per CH» grouy ( 5. 700 s10 ! 44 47 10.37) 11.04) 11.79) 12:37) 12M BD Bw ow OM 
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The same method of calculation was used 


described in section IV, 1 of reference [5]. 


the formation at 25° C of the 


lopentane and cyclohexane derivatives treated 
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For heats of 


his report, the values given in references [8, 
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properties for the 
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the normal alkvleyclopentanes 
hexanes, in the gaseous state, from the elements 
carbon (solid, graphite) and hydrogen (gaseous 

are presented in tables 11 to 19, which give values 
of the following properties, to 1,000° K or 1,500 
K: heat of formation, A//f°; free energy of forma- 
tion, AF/°; and the logarithm of the equilibrium 
constant of formation, log)/v/ 
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\BLE 13 Values of the heat of formation, for the ideal gaseous state, to 1,500° K, for the ? dimethylcycloheranes 








4 
remperature in ° K 
= ( pou For 7 
iz . Jus. lt sO) un) un) ooo Tuo Sun) CLD 1, 000 1100 1. MK we) ws 1 
1 i mula 
Heat rmat Al n keal/mok 
Z 
7 ! eve t i Cl 4 4 ” 4 4.4 1<_0 | Ls »4 2 24 ] 4 140.4 is 4 
etl cve ex CH is. 4 41.1 1.0 4 Wi. 7 is. 4 ivf “2 * w4 wo 44 ist 1.4 wY 
> Dimethyvlevelohexane C.H ) 3. 02 { 7 Wi 1 is 4 1) ; 2 ; l “9 “ 49.2 i8. 2 
x 
irme lohexanes¢ Cl 202 HI 4.21 i7 19.7 1.4 2 ! ‘ 3. 2 2s 2 2 v4 it) 
Dimethvlevclohex Cl 1). 12.2) $2.2 i ‘ iv w.7 l.¢ 2 ¢ 9.4 iv is 
t hexane CH s 2. 22 42. 27 45.4 47.7 iv , ‘ 7 t I wie iv. ¥ iv is. 2 
» j-1) etl hexal Cu " 44.12 4.17 i; uF 2 0 ° “ J H) 4 iu 
” } tal 
RLE 14 Values of the free enerqy of formation, for the ideal qaseous state, to 1,500° K. for the normal alkuyle clopentanes 
rem peratu ! kK 
( pound I . 
im 1f uM 100) uM) int iT sit wn “nw “) 1M " TT " 
‘ 
» ¥ Free em i \/ kca 
tar CH . +o 2 40 Yu. i 29. 2 Oo. 78 4 2 ‘ rT 2.3 x ts “4.4 m4 i4 7. 4s is 4 
- \ al CH ‘ 5 7 x) ; ow 4 s x 7u j s ms. 4 2 . x4 as M48 
jopentar C-H 5 ) m4 2 30. 70. 6 si. 20 2 “f 6S 140. % 65.20 18 Wis 
IM ‘ Cll s 2 4 2 s4 m4 4.7 64 ws . ; 7 ; 2.4 s » 
CoH j 1.69 9 rt . “4 ‘ ; ws «(219 oN f om 
pent ( | t x w.7 a3. 4 in » 9 f Mw) ¢ 2") 246.9 » 4 et 
ent ( Hi j a. 78 a { f f 2.4 s 4 2 2 222. 2 Jis 774.2 2 fi. 2 
| ta { t s ; 0) s 2 - j 4 74.6 in { s s y y s 14 Zita. | = ; ~ 
‘ i CoH 2 12 22.88 | 2 m s 7 12.4 2 24.2 2 4 20 sf ; w)S 
r ‘ " ( iH 24 25. 47 H.4 aS. 4 21.2 { s »” 254.4 JAS 22 2 s a) 12 
ent ( tl 0. 4 4. 97 7. f f) ; 4 1) ru Pd 238 274 1. ts 414.2 j 
pentane ( HT 4 29. 2 2» ¢ ‘ 4 2 4 8, 2 216. ¢ 2 4 » $10.2 $40 is - 
enta Coll ‘ 07 1 69.8 (Ww. 2 v4 0.2) 2 272. 4 f iv i. 2 i ‘ i7s.s a 
Cc [™ ( I! is 2 -. { ‘ su 2 245. ¢ mu 4 i). Wid ‘ 9 
Fi CoH j ° 35. 2 a8 2 6s 214 » 4 2s "2 “4 j ) - x48 
lee clopentane CoH s 2 wa S 7 ~H | 74 ¢ 7.4 24 42 {70.2 ) 6s ( 
i ( 1 4 a. rt a7_ 8 “9 a7. 2 28 ony " 2 44 144.8 if 17s 4y 
~ rcH I I 2 is 20%" 1. 4s 4.4 4 i 4s * ” 2 24 fF 9 »~ ~ 7f 9 
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TABLE 15.— Values of the free energy of formation, for the ideal gaseous state, to 1,500° K, for the normal alkylcyclohex: 
| ! 
remperature in ° K 
Compound For . . . - 
0 208. 16 300 My aw ow 700 SOO) oun 1, 000 1, 100 1, 200 1, 300 1, 400 1. we 
(gas mula 
Free energy of formation, AF/°, in keal mole 
C yelohexane CHy 20.01 +7. 59 7.81 20.66 34.07 47.86 61.85 75.96 90.13 104.30 118.42) 132.58 146.69 160.73 174° 
Methyleyclohexane CHy 26.30 +6. 52 6.79 21.84 37.51 53.55 69.80 86.16 102.59 119.03 135.43 151.83 168.13 184.43 20 
Ethyleyclohexane C.H 2.8 #+9.38 89<£.69 27.12 45.19 63.72 82.50 101.37 120.28 139.23 158.17 177.11 195.93 214.71 23 
n-Propyleyclohexane CoH 32.79 +11.33 11.70 31.5 52.2 73.2) 94.6 | 116.0 | 137.6 | 150.1 | 180.7 | 202.3 | 223.6 | 245.0 26 
n-Butyleyclohexane CywH 36.29 +13.55 13.97 36. 2 59.3 82.8 106.7 130.7) 154.8) 178.8 | 203.0 | 227.1 251.0 275.0 20s 
n-Pentyleyclohexane Cua 9.96 +1560 16.06 0.7 | 66.2 92.3) 118.7 | 145.2 171.8 | 198.4 | 225.1 | 251.8 | 278.2 | 304.8 | 331.2 
n-Hexyleyclohexane Cys 43.64 +17.65 18.15) 45.2) 73.2) 101.7 | 130.6 | 150.6 | IRBS | 218.0 0 247.2 | 27 305.4 334.5 3¢ 
n-Heptyleyclohexane CyHy 47.31 +19.69 20.24 499.6 SO1 TLD) 142.6 174.1 205.8 237.6 20.4 BL] 332.6 364.3 BH 
n-Octyleyclohexane C Hos 4). OS 74 | 22.33 M1 87.0 120.6) 154.5 | 188.6 | 222.9 257.2 201.5 325.8 350.8 304.0 4202 
n-Nonyleyclohexane CsHw 4. 66 79 «24.42 58. 6 04.0 130.0 166.5 208.1 | 239.9 276.7 313.6 350.5 387.0 423.8 4 
n-Decyleyclohexane CH 58.33 +25.84 26. 51 63.1 100.9 139.4 178.4 217.6 26.9 206.3 335.7 375.1 414.3 453.6 492 
n-Undecyleyclohexane CyrHy 62.00 +27.89 2.60 67.6 107.8 WAS 190.4 232.1 274.0 315.9 357.9 300.8 441.5 483.3 52 
n-Dodecylicyclohexane C ysHy 65.67 +29.93 30.69 72.0 114.8 |) 158.3 | 202.3 246.6) 201.0 335.5 380.0 424.5 | 468.7 513.1 | SAT 
n-Tridecyleyclohexane CH —60.35 +31.98 32.78 76.5 121.7 167.7 214.3 261.0 308.0 355.0 402.1 449.1 495.9 342.9 SRO 
n-Tetradecyleyclohexane CoHw —73.02 +34.03 | 34.87) 81.0) 128.6) 177.1 | 206.2 275.5 | 325.0 | 374.6 424.2 | 473.8 | 523.1 | 572.6 622.2 
n-Pentadecy leyclohexane CyHay 76.69 +36.08 36.96 85.5 135.5 | 186.6 238.2, 200.0 342.1 394.2 446.4 498.5 550.3 G024 465 
n-Hexadecyleyclohexane CeHy —S80.37 +3813 39.05 90.0 142.5 196.0 250.1 304.5 | 350.1 413.8 468.5 | 523.1 | 577.6 | 632.2 GRAS “ 
Increment per CH, group 3. 673) +2.048) 2.000 $48 6£.93) 9£.43 11.95 14.48 17.08 19.58 22.12 24.67 29.76; 32 
TABLE 16.— Values of the free energy of formation, for the ideal gaseous state, to 1,500° K, for the seven dimethylcyclohexan x 
remperature in ° K “i 
‘ a Lo 0 208. 16 00) 00 rT) 600) 7H) S00 wo | 1.000) 1,100) 1,200 | 1,300) 1,400) 1% 
Free energy of formation, A F/°, in keal mol 
1,1-Dimethyleyclohexane CH ¥). 93 +8. 42 8.74 26. 6 45.2 1 83.4 102.7 122.1 141.4 160.8 | 180.2 |) 198.5 | 218.4 237 
s-1,2-Dimethyleyclohexane C3H 28.95 +9.85 10.17 27.8 46.1 (4.8 S38 102.9 22.1 141.2 160.4 179.5 198.6 217.4 2 
trans-1,2-Dimethyleyclohexane CH Ww. 91 +8. 24 &. AS 44.6 63.4 82.4 101.5 120.6 139.9 150.1 1781 197.3 216.1 244 
cis-1,3-Dimethyleyelohexane * CH 22.02 +7.13 7.4 43.6 62.4 S14 100.6 119.8 139.0 158.3 177.5 196.6 215.6 234 
frans-1.3-Dimethyvieyvclohexane* CcH +). Of +S. 6S 09 44.5 63.6 82.5 ‘O15 120.6 139.8 158.9 178.0 197.2 216.0 234 
cis-1,4- Dimethyleyclohexane CH 30, ON +0. 07 0. 38 45.5 (4.3 83.4 102.6 121.8 141.2 160.4 179.7 190.0 217.8 238 
frana-1,4-Dimethycyclohexane C.H 31.00 +7. 58 7.89 4.3 63.2 82. 4 101.6 121.0 140.3 159.6) 179.0 198.2 217.2 23 
*See footnotes a and 6 of table 22 
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Figures 1 and 2 show the thermodynamic 
bility of the normal alkyleyclopentanes and 
vievelohexanes in the gaseous state as a func- 
» of the temperature, in the form of a plot of 
standard free energy of formation, per carbon 
m, divided by the absolute 


AFT " 7. 


temperature 


| Free Energies and Equilibria of Some 
Reactions of Isomerization, Cyclization, 
and Hydrogenation Involving Alky]l- 
cyclopentanes and Alkylcyclohexanes 


In figures 3 and 4 are plotted, as a function of 


perature, values of the logarithm of the 


iilibrium constant for the reaction of cyeliza- 
20 Free energies and eq silibrium constants for the 


erizalion of ethyulcyucloherane to the isomeric dimethyl- 


ohexanes in the ideal gaseous state to 1,500 K 
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tion of a given normal paraffin to a normal alkyl- 
cyclopentane or a normal alkyleyclohexane. 

In figure 5 are plotted, as a function of tempera- 
ture, values of the logarithm of the equilibrium 
constant for the reaction of isomerization of a 
given normal alkylevclohexane to a normal 
alkvleyclopentane. 

Values of the logarithm of the equilibrium 
constant for the reactions of hvdrogenation of the 
normal alkylbenzenes [17] to normal alkyleyelo- 
pentanes and alkyleyclohexanes, as a function of 
temperatures, are plotted in figures 6 and 7. 

Figures 8 and 9 show the logarithms of the 
equilibrium constants of the reactions of eveliza- 
tion of normal l-alkenes [18] to normal alkyleyelo- 
pentanes and alkylevclohexanes, as a function of 
the temperature. 

In figure 10 are plotted, as a function of the 
temperature, the values of AF 
zation of ethyleyclohexane into each of the seven 


T for the isomeri- 
dimethyvlevclohexanes. The corresponding  nu- 
merical values are given in table 20 

In figure 11 are plotted, as a function of the 
temperature, for the CH), alkvlevclohexanes, the 
amounts, in mole fraction, of each of the isomers 
present at equilibrium with its other isomers in 
phase. The 


values are given In table 21 


the gas corresponding numerical 


raBLe 21 Equilibrium concentrations of the CoH, alkyl- 
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Figure 12 shows, for the equilibrium between 
cyclohexane and methyleyelopentane, a = com- 
parison of the values derived from this report with 
values reported experimentally by Schuit, Hoog, 
and Verheus [9], Glasebrook and Lovell [10], Shell 
Development Co. [11], and Mizusima, Morino, 
and Fujisiro {12}. 

Table 22 gives, for the relative amounts of the 
cis and trans isomers of 1,2-dimethyleyclohexane, 
1,3-dimethyleyelohexane, and 1,4-dimethyleyelo- 
hexane, a comparison of the values derived from 
this report with experimental values reported by 
Boord, Greenlee, et al. [15]. 

TABLE 22 


( omparison of the calculate d and er pe ri- 


mental values for the relative amounts at equilibrium of 


sone pairs of twomeric cis and trans dime th yle iclo- 
heranes 
Relative amounts 
it equilibrium 
rem 
Compound pera 
ture Caleu Experi 
lated mental ¢ 
liquid liquid 
( 
cis-1,2- Dimethyleyclohexane 250) 21 35 
trans-1,2-Dimethyleyclohexane 7Y 5 
cis-1,3-Dimethyleyclohexane * 225 ri 7s 


trans-1,3-Dimethyleyclohexane 
cis-1,4-Dimethyleyclohexane 


trans-1,4-Dimethyleyclohexane 7s 63 


* This isomer, formerly labeled “trans”, has the following properties (16) 
Boiling point at | atm, 120.09° C; refractive index, ap at 25° C, 1.4206; density 
at 25° C, 07620 g/ml 


This isomer, formerly labeled “‘cis’’, has the following properties [16] 


Boiling point at 1 atm, 124.45° C; refractive index, 1D at 25° C, 1.4284; density, 


at 25° C, 0.7806 g/ml, 


¢ From reference {15 
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actors Affecting Operation of Apparatus for Counting 
Alpha Particles in an Ion Counting Chamber 


By F. J. Davis 


The resolving times of two amplifier circuits for the counting of alpha particles in a 


eviindrical ion chamber is discussed 


stage of the amplifier, gave a resolving time of 2.5 10 


using a l-megohm grid resistor on the first 
= cond 


The effect on the counting rate of the 


and HCl in the ion chamber was found to be negligible. 


however, reduces the counting efficiency by 


The original circuit, using a floating grid on the first 


2 ‘second. An improved circuit 


stage lowered the resolving time to 1.4 10-* 


presence of small amounts of CO, CO»), H.O 


The presence of 1 percent of air. 


20 percent. 


A method whereby the collection time of radon in radium solutions may be of the order 


of 's hour is given. 


ml radium solutions for 17 minutes removes 99.99 percent of the radon. 


of radon by small neopre 


mi 100 en’ 


I. Introduction 


In the operation of the alpha particle counting 
paratus for the determination of radium by the 
don method ! iwo problems that have been in- 
stigated are (1) the resolution of the electronic 
lses resulting from the alpha particles and (2 
quantitative transfer of the radon representing 
radium. 
Factors contributing to the first problem are 
resolving time of the electronic circuit and the 
mposition of the gas. Pulses occuring within 
interval less than the resolving time are re- 
rded as a single pulse. This affects the calibra- 
m constant so that it is a function of counting 
te as well as resolving time. The electronic 
ses are essentially produced by the collection 
lectrons on the center electrode of the ion 
imber. In the presence of gases such as oxygen, 
ch have a high affinity for electrons, the elee- 
ms are captured with the formation of heavy 
As a result, the 


in height so 


itive ions of low mobility 


are broadened and reduced 


t many pulses fail to reach the threshold pulse 


Curt und F. J. Da J. Research NBS Sl (194 RI 


ne meteorlogical balloons was found to be approximatly 


Experiments show that boiling and bubbling nitrogen through 200 


The absorpt lon 
0.03° 


height that will just activate the pulse discrim- 
inator and thus fail to be recorded. 

Factors involved in the second problem (trans- 
fer of radon) are (1) deemanation time, (2) collee- 
The 


deemanation time or the time required to extract 


tion time, and (3) absorption of radon. 
the radon from the solution should be sufficient to 
transfer all radon from the solution to the cham- 
ber, and it will depend on the volume of the ion 
chamber as well as on the amount of solution. 
The collection time, that is, the time allowed for 


the growth of radon in the radium solution, 
should be sufficiently well defined so that the 
amount of radon produced may serve as a measure 
of the radium present. This is especially critical 
when this time interval is short. The absorption 
of radon in rubber is important in the taking of 
breath samples that are usually obtained by the 


inflation of a balloon. 


II. Resolving Time 


The resolving time used in this paper is defined 
bv the relation N=N,e7", where N is the count 


per hour per 107" curie, ¢ the resolving time, / the 
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actual count per hour and NV, is the count 


per 
*“curte for /=0. 

constant NV 
the larger the resolving time. 


hour per 10 
becomes smaller 
The immediate 
advantage of a short resolving time, seen from 


The calibration 


the formula, is that larger samples may be meas- 
ured without having to apply a large correction 
As the statistical 
error of a random count is determined by the 


to allow for unresolved pulses. 


square root of its total count, a further advantage 
gained by a higher counting rate is that greater 
accuracy is obtained in a given time of measure- 
ment. 

If the apparatus used is an ion chamber with 
an amplifying device, such as a FP-54 circuit 
with a galvanometer, the resolving time may be 
of the order of 1 second, thus limiting the radon 
from the sample measured to the order of 107“ 
curie. With this limitation, samples that must 
be fused in a vacuum furnace may require quanti- 
ties of the order of milligrams or iess, thereby 
introducing sampling errors. 

In figure 1 is shown the calibration or count 
" different) counting 


per hour per 10>" curie for 
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Ficure tl. Variation of calibration count per hour per 
10>" curie with counting rate for the original circuil. 


Resolving time caleulated from slope of curve equals 2.5X l0-* see 


rates for the apparatus as described in a previous 
paper? The per hour per 10°" 
plotted as the ordinate on a logarithmic scale 


count curie is 
and the counting rate as the abscissa on a linear 


seale. It is to be noced that the calibration con- 


stant .V drops from 257 counts/hr per 10~™ curie 
for small counts, to 239 for a counting rate of 
100,000 counts/hr, which is the equivalent of 
4.2<10~-" curie. 

Calculating ¢ from the slope of the line in figure 
1 gives =2.5%10~ see. Although this resolving 
time is adequate for breath samples and room air 


2 See footnote | 
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samples where the amount of radon rarely exc eo 
5x<X10°"" requires the use of the ca! 
bration curve for samples over this value. 

In the amplifier used to secure the data show 
in figure 1, no grid leak was used on the fips 


curie, It 


amplification stage (type 38 tube). To shorty 
the resolving time for this cireuit, the obvioy 
method was the placing of a 1l-megohm gr) 
leak between the grid of the 38 tube and grown 


The effect of this is shown in figure 2. T| 
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Fie. 2. Variation of calibraiion count per hour per It 
curie with the counting raie for short resolving 
circuil, 


Resolving time calculated from slope of curve equals 14x 10-¢ se 


circles represent independent samples of radon The crosses repre 


different recounts of a single sample while decaying during a perio 
1 week 


calibration count per hour per 107" curie drop 
from 252 for small counting rates to 238 for 
counting rate of 1,400,000 counts/hr. It is 

be noted that the counting rate in figure 2 cove: 
a tenfold range of that shown in figure 1. Ca 
culating the reselving time from the slope of th 
1.4 1074 see. 

These two curves are for the same cham) 
38 tube and amplifier. For the first curve, thy 
for the 38 tube and amplifier was 4 
Waddell In our amplifier 
coupling condenser of 500 wuf was used for th 
plate of the 38 tube to the grid of the followin 
For the second curve the changes mad 


line gives ¢ 


cireuit 


suggested — by 


stage. 
in the cireuit of the 38 tube were (1) replacing th 
500-uuf coupling condenser on the plate with 
100-uuf condenser, (2) connecting a 1-megohn 
grid leak to the ground, (3) limiting the eathod 
emission by placing a resistor of about 30 ohm 
in parallel with the heater, and (4) the 38 tub 
was connected as a triode, that is, with the sere: 
grid and plate connected together. 

The original circuit with the longer resolvi 
time has a greater efficiency of counting at slow 


‘KR. C. Waddell, Rev. Sci. Instruments 1@, 311 (1939 
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arc 


as shown by comparing figures 1 and 2. 


s increased efficiency is probably due to a 
signal to noise ratio, which is somewhat 
As the number of 


ses that are below the threshold pulse height 


cher for the original circuit 


either cireuit is small, the difference in the 


ficiency of counting by the two circuits amounts 
only 2 perent. 

In figure 2 the calibration curve for the short 
shows both circles and 


time cireuit 


The cireles are independent fillings of the 


psolving 
unbers from two different standards, one con- 
ng 10 
other 1.077 107° ¢ Ra for the higher counting 


'o Ra for the lower counting rates and 


The crosses are all from a single filling and 
esent different recounts, beginning on the right 
| proceeding to the left, as the radon decayed in 
for A, the decay 
The decay of 


chamber. The value used 
stant of radon was 0.00755/hr. 
radon was followed for 1 week, and each cross 
hours 


s the average of a 12-hour run with 24 


tween the beginning of each run, except  be- 
en the third and fourth crosses from the right, 
en it was 48 hours. 

\ side tube containing water was sealed on a 
umber and the temperature of the water con- 
No effect on the counting rate was noted 


led 


temperatures below room temperature. At 
mm temperature spurious pulses occurred, prob- 
across insulators in the 


v due to leakage 


umber. A single drying trap using CaCl, re- 
wes water vapor sufficiently for satisfactory 
eration 

Hvdrogen chloride of the amount driven off 
m deemanating a 5-percent HCI solution has no 
fect on the counting rate. However, if strong 
mixtures 
used, the counting rate is reduced. A trap 


d with Asearite followed by a drying trap can 


ICl solutions such as constant boiling 


ised to remove the hydrogen chloride. 

The CO and CO, were taken from commercial 
s cylinders and passed through the purification 
stem described in a previous paper,* which re- 
oved the oxygen and water vapor. In the case 
CO, an additional trap of Ascarite plus a 
ving agent was used to remove any CQ), present. 
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Ill. Effect of Composition of Gas 


the 
presence of carbon monoxide, carbon dioxide, sir, 


The effect on the calibration constant of 
and water vapor was determined. Considerable 
quantities of CO are formed by a carbon-crucible 
vacuum-furnace when oxygen is liberated from the 
samples fused. The effect of using CO in place of 
nitrogen in the ion chamber was to lower the 
calibration constant by only 1 percent 

A small amount of CO, seems to decrease the 
resolving time of the apparatus. The effeet of 
the addition of CQO, is shown in figure 3. The 
ordinate represents the resolving time in 10~¢ 
sec, and the abscissa represents the percent by 
The 


time falls to a minimum near 10 percent of CO, 


volume of CO, in the chamber. resolving 
whereas at 20 percent of CO, it is about the same 
as pure nitrogen and at greater percentages the 

The effect of CO, on the 
for concentrations less than 


resolving time rises. 
calibration constant 
20 percent at a counting rate of 250,000 counts/br 


1 


is less than 's percent. 
The effect of air in the chambers is vy ery marked; 
1 percent of air reduces the counting efficieney at 


30,000 counts/hr by approximately 20 percent. 


x10 


SECONDS 
b 








o 





he circles represent observations made t Ising the short resolving time 
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IV. Deemanation Time 


The method used for deemanating solutions is 
to boil the solutions under reduced pressure of 
approximately 1/2 atmosphere in a reflux con- 
denser while bubbling N, through it as the transfer 
agent to carry the radon into the chamber. As 
the chamber is filled to atmospheric pressure, the 
pressure and boiling temperature in the condenser 
rises accordingly. The time required to fill the 
chamber is usually 16 to 18 minutes in order to 
efficiently remove O,. 

A test was made to measure the efficiency of 
transferring the radon from a 200-ml solution for 
various deemanation times into a 3.6-liter cham- 
ber. This was accomplished by filling a chamber 
in a specified time and immediately on completion 
reboiling the same solution and filling a second 
chamber to measure any radon remaining in the 
For a 17-minute deemanation time, no 
radon was found in the solution in excess of that 


solution. 


amount that would grow from the radium in the 
solution in the short time between the two runs. 
The efficiency of the above test was such that 
one-hundredth of a percent 
detected. For a 10.5-minute run, 0.020 
of the radon remained in the 


have been 
+ 0.003 
solution. 


would 


percent 
For a run of 6 minutes 1.01 percent of the radon 
remained, 

The results of a slightly different experiment 
to determine the «amount of radon transferred 
during a normal run of 16 minutes is shown in 
figure 4, where the amount of radon removed from 
io-* 


is plotted against the time from the start of the 


a solution containing 1.62 > curie of radon 


deemanation. In this case, the flow of nitrogen 
through the solution was adjusted to a rate that 
would normally fill a chamber to atmospheric 
pressure in 16 minutes. It was then allowed to 
for the first number 
one; for the second minute into chamber number 


flow minute into chamber 
two; for the next 2 minutes into chamber number 


three; the next 4 minutes into chamber number 
four, and the remaining S minutes into chamber 
number five. After the deemanation was com- 
pleted, the chambers were filled to atmospheric 
pressure with No. The amount of radon trans- 
ferred at the termination of the filling of any one 
chamber would then be the sum of the amount of 
radon in the chamber and the preceding chambers 


548 














x1o7!0 1 _ ome 
- a 
15 
.. 
a 
’~ 10 
ie 4 
~ 
oO 
o5 
s 
. at 
oO 5 10 5 2 
MINUTES 
Figure 4.—Curve showing radon fransfe rred from refi 


condenser to ion chamber during deemanation. 


lotal amount of radon transferred 1.62 10°" curie 


the first minute of deemanation th 


are being adjusted, consequently t! 
first point may vary considerably for differe: 


During 
stopcocks 


runs. According to the curve in figure 4 th 
radon transferred in the first 6 minutes is 
percent. As only a portion of the nitrogen ha: 


passed through the solution to fill a chamber, les 
radon is transferred than when all the nitrog: 
passes through the solution to fill a chamber as 
the earlier experiment. 


V. Collection Time 


The time allowed for the collection of the rado: 
from a sample is the interval between a_ pre- 
liminary deemanation the filling of th 
chamber by a second deemanation process. Thi 
interval is measured from the instant when tly 
stopcocks are closed on the reflux condenser at th 


and 


end of the first deemanation to the exact insta! 
during the subsequent filling of the ion chamber 
when all the radon destined to enter the cham! 
The colli tn 


tion time terminates before the ionization chambe! 


has been removed from the solution. 


stopeock is closed, as it takes some time for th 
radon to pass through the purification system 
The effective point of termination of the collectio 
is fixed according to the following consideratio 
first that the 
filled to atmospheric pressure and the purifying 


Suppose ionization chamber \ 
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wsiem evacuated. If then the stopeock 73 in 

re 5 is opened, the mercury manometer will 
ome to a level B depending on the relative 
olumes of the ionization chamber and the purify- 

system. Vice versa, in our actual operation, 
vhen the manometer reaches the level B, we know 
hat all the radon the purifying 
.ystem and that radon alone will be swept into the 
is filled up te 


contained in 


mization chamber when this 
tmospheric pressure. All that is necessary to 
termine the time of the second deemanation is to 
cord the time when the mercury crosses the 
wrk B while filling the chamber. The radon 
ccumulated in the additional time of filling fails 
o reach the chamber. 

Using this interval for short time collections 
woved to be accurate within the statistical error 


f O03 percent, even for a collection time of 30 


nutes. 
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VI. Absorption of Radon 


The absorption of radon by small neoprene 


eteorological balloons, such as are sometimes 
ised in transferring samples, is shown in figure 6 
The percent absorption of radon is plotted against 
time of contact with the neoprene. The area of 


he balloon used was about 2,500 em?*. Using 
this value, the amount of absorption is approxi- 
mately (0.03°%/min)/100 em*. The 


vas measured by placing a known amount of 


absorption 


vion in a balloon, allowing it to stand for the 
then 
The amount of 


iod of time under consideration, and 


Measuring the remaining radon. 
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Area of balloon 2.500 em?’. 
radon lost under these conditions would be due 
not only to absorption but also to transmission 


through the neoprene. 


VII. Additional Remarks 


A satisfactory arrangement when a sufficient 
number of chambers are available is to reserve 
some chambers for large samples and use the short 
resolving time circuit. The remaining chambers, 
if reserved for small samples, can be used either 
with the short or long resolving time cireuit. A 
sample of 10 curie was used to calibrate one 
chamber connected to a short resolving time cir- 
cuit and its value was found to fall on the same 
line of calibration points used for smaller rates of 
counting with this arrangement. This indicates 
that quantities of this order can be readily meas- 
ured within the same limits of error as smaller 


samples. 
The writer acknowledges the helpful criticism 
of L. F. 


ments by Christine Shaffer and Sarah Culpepper 


Curtiss and the assistance in measure- 


WasHinaton, July 28, 1947. 
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A Statistical Analysis of Some Mechanical Properties 
of Manila Rope 


By Sanford B. Newman and J. H. Curtiss 


Breaking strength, weight, and circumference are three important qualities that are 


determined when fiber ropes are submitted for test. 


These properties are all subject to some variation because of differences in fiber quality , 


method of fabrication, and previous treatment. The results of tests on more than 800 


samples of 3-strand manila ropes ranging in size from *j, in. to 3 in. in diameter are tabulated 


and analyzed by the methods of modern mathematical statisties. 


Considerable emphasis is 


placed on the rationale and details of the methods of analysis, as they are considered to be 


applicable to a broad variety of similar projects 


I. Introduction 


Strength, weight, and size are usually deter- 
mined for samples of fiber rope submitted to the 
These properties 


Bureau for acceptance tests. 


have been found to be subject to some variation 

manila rope by Stang and Strickenberg [1].' 
This variation would appear to be caused by differ- 
neces in fiber quality and in the mode of fabrica- 
ion of the rope. 

In many applications of fiber rope where even 
small economies in weight are important, 1. e., in 
use of ropes on cargo airplanes, a knowledge of 
the probable strength of a set of ropes can result 

an increase of the pay load. For other applica- 
tions Where ropes pass over sheaves and through 
eves, the probable range of size is an important 
design consideration 

The manila ropes discussed in this paper were 
submitted for test by a Government agency from 
IN3S to 1941, 
by rope works and rope contractors in the 3 years 


They represent material supplied 


st prior to the loss of sources of manila fiber in 
World War II 


neous sizes of the samples available for the various 


This accounts for the heteroge- 


nominal diameters, as the smaller sizes and the 


tegral multiples of \, in. in diameter are ordered 


tore frequently than other sizes. Although some 


ires in brackets indicate the literature references at the end of this 


Properties of Manila Rope 


data manila 
they were not of sufficient quantity to allow statis- 


tical 


were available for 4-strand rope, 


treatment and consequently are not dis- 
cussed here. In general, the requirements for 
4-strand rope in the Federal Specification [2] will 
yield a basis for comparison with the 3-strand 
ropes treated in this paper 

The methods of test described are those found 
in the Federal Specification for manila rope [2], 
but the data are applicable to many situations 
where the conditions may be somewhat different. 
Whittemore [3] found that the speed of the moving 
head of the testing machine between 1 and 4 in. 
mio had little effect on the observed strength of 
the rope. It has also been noted [4] that measure- 
ments of circumference with increasing loads on 
the specimen yielded decreasing changes in cir- 
cumference for equal increments of load up to 
loads equal in pounds to 300 times the diameter 
in inches squared. 


II. Methods of Test 


1. Circumference and Weight 


Both weight and circumference of manila ropes 
were determined on an unspliced sample that had 
been conditioned in an atmosphere of 654 
+ 2° F for at least 
The samples were 


2 per- 
cent relative humidity and 70° 
72 hours preceding the test. 
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long enough to provide the minimum free length 
specified in table 1. 


TaBLe 1*,—Length of weight—circumference Specimens 
; Minimum 
( nference ' 
reum me free length 
! Ft 
Below 3 a 
ito 6, inel 5 
Over 6 
Values from Table IIL of Federal Specification for Rope, Manila 


horizontal testing 
200 JF |b (where 
DP) is the nominal diameter of the rope in inches) 


The rope was placed in a 
machine and a tensile load ?P 
was applied. A single fiber was passed snugly 
around the rope, near the middle of the sample, 
The cut length of 


fiber was measured and the circumference recorded 
1 


and cut where it overlapped. 


to the nearest '\¢ in. 
With the load, P, 


length as specified in table 1 was marked off on 


still applied to the rope, a 


the free length, and the load was then removed, 
The marked length was cut from the sample and 
weighed, and the weight in pounds per foot was 
computed, 

2. Breaking Strength 


A breaking strength sample had an eye splice 
at each end and measured from 5 to 6 ft between 
the inner ends of the splices. The samples were 
conditioned in the same fashion as the cireum- 
The then 
removed from the conditioning room and_ the 


ference-weight samples. ropes were 
splices were immersed in water for 15 minutes in 
order to minimize the possibility of a failure in 
the splice. 

Ropes with breaking strengths over 2,000 Ib 
were loaded to failure in a horizontal hydraulic 
testing machine, and all others were tested using 
a horizontal, pendulum, serew testing machine. 
The ropes were loaded by means of 3-in. diameter 
steel pins passed through the eyes at each end 
of the sample. The speed of the moving head 
of the testing machine was 3 in./min during the 
The until at 


the maximum load the rope failed in one or more 


tests. tensile load was increased, 


strands. 
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III. Methods of Statistical Analysis 


Four fundamental relationships were studiod 
(1) Circumference as a function of nominal diam- 
eter, (2) weight as a function of nominal diame er 
(3) strength as a function of nominal diamete 
and (4) strength as a function of weight. Thy 
first three are useful in preparing specifications 
and predicting the characteristics of individual 
ropes purchased under a contract or order. |) 
was therefore considered desirable in these cases 
to furnish careful predictions of the dispersions 
encountered in practice, as well as of the mean 
values, so that realistic tolerances for individual 
The fourth relationship is 
important as a measure of innate characteristics 
of the material. 


ropes can be set up. 


Central tendeney, rather than 
the dispersion, seemed to be of paramount interest 
in this case. 

Observations on 863 ropes were available for 
this study. With truly random sampling and a 
clearly defined and homogenous universe or popu- 
lation, such a large sample, if properly handled, 
should provide close estimates of the underlying 
frequency distributions and relationships of the 
variables. In the 
empirical equations containing, if necessary, sev- 


such circumstances use of 


eral parameters would be justified. In the present 


instance, no direct control could be exercised 
over the sampling, and the universe, which pre- 
sumably consists of the entire National output 
of manila rope of the relevant nominal sizes during 
the period from 1938 to 1941, from 


homogeneous. The data, as might be expected, 


was far 


contained certain anomalies that are more or less 
directly attributable to the composite nature of 
this universe, or to the nonrandom sampling. 

In spite of these shortcomings, the observations 
appeared to exhibit a sufficient amount of internal 
consistency to warrant the use of mathematical 
methods, provided that the description of the 
underlying causal situation attempted thereby 
should not be too elaborate. 

The mean values of the dependent variable were 
represented in each case by a curve of the genera! 
form 


~~ 


eX”. (1 
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hore Y designates the (arithmetic) mean value of 
he dependent variable and Y the independent 
ariable. The parameters & and 6 were adjusted 
othe data, except in the case of the circumference- 
liameter relationship, where 6 was arbitrarily 
aken as unity. 

The curves were fitted by the following general 
wocess: All data were first transformed to loga- 
ithms (to the base ten), and then a regression 
quation of the form 

log Y=a+b log NX, 

least 
At the 


same time the standard error of estimate s (that is, 


vas fitted DY the standard unweighted 


sjuares method to the logarithmic data. 


the root-mean-square deviation of the data about 
The 


esulting equation (2) yielded an optimum esti- 


the regression line) was also determined. 
5 


mate of the mean log Y for each log Y, provided 
that it can that the 
between mean log Y and log Y is exactly linear 


be assumed true relation 
ind that the standard deviation of the variable log 
Y for a given log XY is a constant independent of 
Y. As the antilog of the mean of a set of loga- 
rithms of Y is not the arithmetic mean of Y itself,’ 
t was necessary to add a correction of some sort 
The 
1.151298", 
which is based on the further assumption that log 
) is normally distributed for each value of log X 
see [5], pp. 120-121). 
tion used in 


0 (2) in transforming back to the form (1). 


correction chosen in this case was 


Thus in terms of the nota- 
connection with (2), formula (1) 


becomes 


10 ari 1512982 YD 


The calculations mvolved in the curve-fitting 


process were carried out almost entirely on 


machinery at the Computation 
Bureau. The logarithmic 
was accomplished automatically 


punched-card 
laboratory of the 
transformation 
by the use of master logarithm cards in conjunc- 
tion with a collater that simultaneously punched 
both the logarithms and their squares onto the 
data cards. Cross products and cumulative sums 
were then obtained in the standard way on multi- 
plying punches and tabulators 

\ few 


remarks on the underlying 


ationale of this method of curve fitting and the 


—_—_ — 


general 


s the geometric mean of Y, however. 
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choice of the type of curve are in order, as the 
situation about to be described is typical of that 
frequently encountered in materials. 
The most elementary physical theory of the data 
would suggest that mean observed circumference 


tests of 


should be directly proportional to the nominal 
diameter (which is actually defined to be one- 
third of the that 
weight then be 
directly proportional to the square of the observed 
This would ordinarily imply that 
the mean values of circumference, weight, and 


nominal circumference), and 


observed and strength should 
circumference. 


strength, as functions of nominal diameter, could 
be most appropriately estimated by fitting linear 
trends by the classical least squares method to 
the observed values, respectively, of circumfer- 
ence, the square root of weight, and the square 
root of strength. 

Unfortunately, as is so often the case, the data 
for all of these variables reveal an unmistakable 
dispersion of the 
(measured, say, by their standard deviation) to 
This would necessitate 
The 


simple unweighted least squares solution would 


tendency for the readings 
increase with their mean. 
the use of a weighted least squares solution. 


give an inefficient estimate of the mean and would 
not yield any over-all estimates at all for the 
But 
the difficulty of determining a posteriori the 


variances within nominal diameter groups. 


proper individual weights to use in such cases has 
led in recent years to the widespread use of trans- 
formations of the scale, such as the logarithmic 
effect, at least 
theoretically, a stabilization of variance over the 


transformation here used, to 
Such trans- 
formations under certain circumstances have the 
additional property of rendering the distribution 
of the transformed variable more nearly normal 


whole range of values of the mean. 


or Gaussian. 

In the present case, as far as could be deter- 
mined from the data, the standard deviation of the 
observed circumference seemed to be roughly a 
linear function of the mean. The transformation 
ordinarily employed in that event is a logarithmic 
transformation [5]. Preliminary exploration §re- 
vealed that this transformation indeed seemed to 
stabilize the variance of the circumference meas- 
It also seemed to stabilize the variance 
This 


would follow mathematically if the transforma- 


urements. 
of the weight and strength measurements. 


tion really did stabilize the circumference variance 
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and if mean weight and strength really were 
direetly proportional to some power of the circum- 
ference. Thus additional evidence was furnished 
thereby as to the correctness of the choice of the 
logarithmic transformation for the circumference 
data. 

As for the choice of. the curve, it is natural to 
choose a type that contains as few parameters as 
possible and that can easily be handled in conjunce- 
tion with the logarithm transformation. 
and Stickenberg [1] chose to use a one-parameter 


Stang 


curve of the form 
Yok N (N41) =k N44 N, 


to represent the strength-diameter relationship. 
Although this equation is readily adaptable to the 
logarithmic least squares approach (it seems to 
have been fitted empirically in {1]), nevertheless 
the first degree term is not easy to explain in terms 
of the physical theory. Exponential equations 
of the type (1), on the other hand, not only can be 
conveniently handled after a logarithmic trans- 
formation, but also accord very well with a slight 
extension of the simplest physical theory of the 
tests. They imply that weights and strengths of 
a series of ropes of different sizes whose successive 
nominal diameters are in a constant ratio will 
themselves be in a constant ratio. 

As a check on the adequacy of (1), three-pa- 
rameter curves of the type 


were actually fitted to the weight-diameter rela- 
tions, and in spite of the additional parameter, 
the results obtained were approximately equivalent 
to those obtained with (1). 


TABLE 2. 
Dependent variable Independent variable 


Sym- 


Meaning uo Unit Meaning bol Unit 
Circumference Cc 1/16 in Nominal diameter D 1/16 in 
Weight Ww Thy ft Nominal diameter D 1/16 in 
Strength S Ib Nominal diameter D 1/16 in 
Strength - S Ib Weight Ww sib/ft 


IV. Results and Discussion 


The various equations for mean values, obtaiiie 
as described in section ILL, together with the cor. 
responding values of the regression constants a. 4 
and s associated with eq 2 of section IIL are a! 
given in table 2. 

In interpreting the equations, it should be und» 
stood that they give the estimated mean values o| 
the dependent variable for each fixed, predeter- 
mined value of the independent variable. Thus 
the last equation, which gives the estimated mea 
of Sas a function of W, cannot be obtained by 
merely eliminating D between the second an 
third equations, because the distribution of values 
of S corresponding to a fixed D, say D=Dy an 
the distribution of values of S corresponding to 9 
fixed W, are not in general quite the same, even 
if the fixed W was determined by setting D=/) 
in the equation. (Actually the 
methods of deriving the fourth equation in table 2 
happen to agree to two significant figures in the 
exponent and three in the constant factor, owing 
to the relatively small seatter of the data about 
the various curves.) 


second two 


The standard errors s, and s, of a and 6 are also 
presented in table 2. It will be noted that these 
standard errors turn out to be exceedingly small 
This is partly due to the large number of items 
entering into the sample, and partly due to the 
relatively small size of s in each case. Statistica! 
theory would state that if the sampling had been 
truly random and if curves of the type used ex- 
actly described the mean relationships in the 
sampled universe, then ranges of 
b+3s, would, with high probability, contain the 
corresponding “‘true’’ values of a and 6; that is, 
the values for the universe. Due mainly to the 


a+3s, and 


Regression equations and statistics of distributions of circumferences, weights, and strengths of 3-strand manila roy 


Regression statistics * 


Equation of means 


a } s s 
C=#3.119D 0. 49353 (b) 0. 02026 0. 00089 
We 0.001447 Diet —2. 84023 1. SSU6S O2541 O03 O. 
S=70481 Dit LMag9 1. 828193 | . 042619 . 0064 i" 
Sa 4027S W' ver 4. 602511 0. 9880 | . O47059 002s | (mt 





* For explanation of a, 6, and 4, see discussion of eq. 2 in sec, IIT, 


tion that for each value of J’, the data constitute a random sample from a universe of such data, 


b Not adjusted to data, 
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The symbols «, and s, denote the standard errors of a and 6, computed under the assum) 


(See see IV.) 
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that the method of sampling was not rigor- 
y random, such a statistical interpretation of 
nd s, is apparently unwarranted. It is be- 
d, however, that these standard errors afford 


ile basis for the following important general 


ements about the coefficients and exponents in 
equations in table 2 


1) In 


elationship, the coefficient of D is significantly 


the case of the circumference-diameter 


rver than 3, which is the value suggested by the 
lefinition of nominal diameter 
2) In the weight-diameter and strength-diam- 
ter relationships, the exponent of PD is signifi- 
antly less than the value 2 suggested by the 
simplest physical theory 

}) In the 
xponent of Wis significantly different from unity. 


strength-weight relationship, the 


If “true”? values of a and 6 may be postulated, 
hen it follows that there exists a ‘“‘true’’ value of 
he ordinate of the equation of the mean for each 
riven value of the independent variable. — Investi- 


ration of the appropriate standard errors reveals 


that theoretically with high probability the mean 
alues of C, WY and 
responding “true” 


S do not deviate from the 
values by more than 2 per- 
ent in the case of C and W, and 4 percent in the 
ase ofS 


ith much caution, 


But this statement must be interpreted 
not only because of the non- 
indomness of the sampling, but also because even 

the curves had been fitted in some analogous 
nner to the universe instead of to the sample, 

is possible that deviations of this order from 
the actual true means of the universe might be 
observed, because of the fact that the simple type 
for fitting 
lescribe the real functional relation between the 


f curve chosen may not accurately 
true means and the independent variables. 

The preceding discussion of standard errors may 
be summarized by stating that from the viewpoint 
f theoretical statistical analysis, the curves have 
been fitted with a considerable amount of preci- 
sion; but in default of exact knewledge as to the 
nechanism real 
weuracy of the curves as a description of the 


of the sampling method, the 


rational manila rope technological situation during 

the data period must be taken largely on faith. 
The information in table 2 has been tabulated 

numerically in table 3, and represented graphically 


n figures 1, 2, 3, and 4. The ‘Estimates of 
\lean’’ columns in table 3 and the central curves 


in all of the graphs were obtained by straightfor- 
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ward substitution into the equations in table 2. 
The small that graphs 
represent the means of the observations for the 
independent. 


circles appear on the 
respective indicated values of the 
variable. 

The information as to dispersion of individual 
observations about the curve of means is presented 
in table 2 in the form of the standard error of 


estimate s. For convenience in the applications, 


this dispersion information is presented in table 3 
in terms of tolerance limits obtained in the follow- 
ing manner: 

Reverting to the notation of equation 2. of 
section IL], if it be assumed that the value of log 
Y given by eq 2 is the true mean of values of log 
and s is the true standard 
ts, log V4 ts] 
will bracket a fixed proportion p (f) of the under- 


) for each value of NX, 
deviation, then the interval {log Y 


lving distribution, the proportion being dependent 
only on the value of ¢f. (The trivially small size 
of all standard errors involved in the present case 
makes these assumptions entirely tenable from the 
theoretical point of view. In the present case, 
the value of f was so chosen that if the distribution 
of individual value of log ) were normal or Gaus- 
sian, the value of pt would be 0.95. Specifically ; 
this value of f to 6 decimal places is 1.959964. 
Thus the formulas for the tolerance limits appear- 


ing in table 3 were: 


10 19500648 4 
] ()! 9599645 i” 


Lower tolerance limit 
Upper tolerance limit 


In tabulating the numerical value of the toler- 
ance limits given by these formulas, the general 
practice was to round off to the number of signif- 
icant figures appearing in the raw data. However, 
a number of exceptions were made, chiefly in the 
direction of retaining one extra significant figure, 
to conform with the conventions of tabulation. 
Ambiguous cases were always rounded outwards. 

The tolerance limits are plotted on the graphs 
as the outer curves. It should be emphasized 
that these outer curves pertain to individual measure- 
ments, and not to means, such as those re prese nted 
by the circles on the graphs 

In general, the analytic representation of the 
863 observations is remarkably faithful, as can be 
seen from a glance at the figures. (In interpreting 
deviations of the circles from the central line, the 
varying number of observations represented by 
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TABLE 3 Circumference, weight, and strength of 3-strand menila rope 
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each cirele must be taken into account.) The 
representation of dispersion is apparently as sue- 
As a partial 
check it might be noted that, for example, 34 
observations on S fell outside the tolerance limits 
for S given in table 3; this is 4 percent of the total 


cessful as that of the mean values. 


number of observations on S, a proportion that 
compares favorably with the specified theoretical 
The fact that the number of 
observetions that were too high was about equal 


value of 5 percent. 


to the number out on the low side gives added 
credence to the validity of the logarithmic trans- 
formation. 

of Wand C. 

The following irregularities in the data should 
be explicitly noted: 

(1) The observed distributions of C, W, and S 
for nominal diameters of ® in., 1, in., and 1% in., 
not found in the Federal Specification [2], were 
very similar to the corresponding distributions for 
the nominal diameters of '\\. in., 16 in., and 1"), 


Similar situations obtain in the cases 


in., respectively. In each case, the means and 
tolerance limits given in table 2 for the larger 
size represented the mean and range of the obser- 
vation for the smaller size (as well as for the larger 
size) much better than the mean and tolerance 
limits given in table 2 for the smaller size. Ap- 
parently, rope of the appropriate one of the three 
larger sizes is supplied whenever one of the three 
smaller sizes is stipulated in a purchase. 

(2) The observed values of C, W, and S for a 
fixed value of D exhibited some tendency to occur 
in clusters having smaller individual dispersion 
than that represented by the values of s in table 
2 or the tolerance limits in table 3. This phe- 
nomenon was undoubtedly due in part to nop- 
random sampling. 

The chief result 
noted in (1) above, and of the clustering noted in 
(2), is to increase the values of s (and thus in- 


of misclassification such as 


crease the spread of the tolerance limits) over the 
values that would have been obtained if such 
irregularities had been absent. 

(3) In the case cf the *\. in. nominal diameter, 
the observed circumference and weights of the six 
ropes tested were not properly represented by the 
analytical expression. In view of the various 
testing, rounding off, sampling, and classification 
errors involved in the measurements, it may well 


be that extrapolation from the analytic represe 
tation based on over 800 measurements shou! 
furnish more reliable information than a samp) 
of six possibly anomolous observations. Howe , 
for completeness, the means and ranges of th: 
and W observations for ))=*\. are given in tabk 


TABLE 4. Summar " of observed circumfe rences and we 


for %\@ in. 3-strand manila rope 
[Based on 6 observations on each variabk 


Circum- 
ference 


ie am 
Arithmetic mean 11.2 
Minimum 10 
Maximum 12 


Figure 5 exhibits a comparison between t! 
results on mean strength obtained in the prese 
paper, and the curve fitted by Stang and Stricke 
berg [1] to the data on strength which they ob 
tained in 1921. An examination of the closenes 
of fit of Stang and Strickenberg’s curve and of th 
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Figure 5.—Comparison of present sample with that 
reference |1| 
The solid line is the present sample. The broken line is the sample of 
ence [1] 
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spersion of their data, indicates that the theo- 


repre 
s she 
1 Samp 
lowe 
of th 
tabl 
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ud we erage 


it one-half as many observations. 


breaking strength of 


tical accuracy of their curve should be roughly 
mparable with that of the corresponding curve 
the present paper when due allowance is made 


the fact that their curve is based on only 


It follows 


for sizes of 1 in. diameter and greater, the 


rope in 1921 as 


presented by Stang and Strickenberg’s curve is 


nificantly lower than that of the rope discussed 


the present paper. 


en tl 
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Voltage Anomalies of the Glass Electrode and the 


Chemical Durability of the Glass 
By Donald Hubbard and Gerald F. Rynders 


Additional work has been done on the correspondence between voltage departures 
errors) of electrodes prepared from Corning 015 glass and the change in the chemical dura- 
bilitv of the glass The results confirm and extend the earlier findings that voltage depar- 
tures of the electrode are accompanied by durability changes of the glass. The pronounced 
voltage departures in the alkaline region and in hydrofluoric acid solutions are accompanied 
by conspicuous attack of the glass. Dilute solutions of acids that do not form soluble 
compounds of silica cause swelling of the glass. This swelling is repressed as the concen- 
tration of the solutions is increased. As is to be expected from the distribution law and 
Donnan membrane considerations, this repression of swelling was also brought about by 
high concentrations of organic acids and salts as well as by the so-called “strong’’ acids 
The attack of Corning 015 glass by alkaline solutions of different ions was found to be in 
the order, Na* Li* K* Ra* NH,*, which is the same order dictated for these ions 


by the voltage departure data. 


CaO 6 percent (Corning 015) [1, 7]. Although 
this glass combines the two necessary properties 


I. Introduction 


eee lel —— ee we rt 
Che ability of a glass to function satisfactoril; of adequate hygroscopicity and uniform durability 


s an indicator of the hydrogen ion activity of 
queous solutions appears to be determined 
irgely by two properties of the glass, namely, 
he hygroscopicity (water sorption) and the uni- 
ormity of the chemical durability over an 
xtended pH range. Glasses of inadequate hygro- 
wopicity produce electrodes of high resistance 
hose pH responses fall appreciably below the 
heoretical dictates of the Nernst equation: 


AE=0.000198TApH [1, p. 78, 2, 3]. 


On the other hand, solutions that cause detect- 
ible change in attack of the glass [4] also bring 
bout voltage departures (errors) from the straight- 
ine relation of the Nernst equation [1, 5, 6}. 

The glass whose hygroscopicity and durability 
haracteristics have most nearly fulfilled the 
quirements for pH measurements is one of 
mposition SiO, 72 percent, Na,O 22 percent, 


ure n bracket neicate the terature references ai the end of 
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Voltage Anomalies of the Glass Electrode 








to a considerable degree, electrodes prepared from 
it exhibit voltage departures in the “super acid” 
range [5] and even more conspicuous errors in 
alkaline solutions above pH 9 [1]. Such voltage 
departures ih alkaline solutions covering a limited 
range of conditions and pH values have been 
reported as being roughly directly proportional to 
the amount of attack [8, 9]. 

The present investigation was undertaken to 
ascertain if electrodes prepared from Corning 015 
glass always have voltage departures accompany- 
ing changes in the chemical durability of the glass 
and also to determine to what exten. the magni- 
tude of these voltage departures is directly pro- 
portional to the change in the attack of the glass 


by the solutions 


II. Experimental Procedure 


The glass electrodes were prepared by blowing 


a thin bulb of Corning 015 glass on the end of 
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The inner electrical connection 
was made by filling the bulb and tube with 
mereury [10]. Such metal-filled are 
simple to prepare and have the advantage that 
an electrode broken in service does not seriously 


soft glass tubing. 


electrodes 


contaminate or alter the pH of the solution under 
observation. All voltage and pH measurements 
were obtained at room temperature by means of 


a Beckman pH model G. 


Unless otherwise stated, the hydrogen electrode 


meter, laboratory 
was used as the reference electrode, and the value 
in 0.02 N hydrochloric acid as the zero reference 
This solution was chosen to avoid the 
the equilibrium the 
hydrogen electrode and M/20 potassium 
phthalate [11]. Further, the acid solution served 
to restore the glass electrode after severe abuse to 
its initial condition {12}. 

The chemical durability (reported as attack) of 
the glass was determined by the interferometer 
method [4, 13]. This procedure offered a rapid 
means of determining the amount of solution or 


voltage. 
uncertainty of between 


acid 


swelling of the glass over a wide range of condi- 
tions under circumstances in which the conven- 
tional analytical methods would be time-consum- 
ing, unreliable, or impossible. ? In_ this 
method the magnitude of the attack was deter- 


even 


mined by observing the displacement of the inter- 
ference fringes when the specimen was placed 
under an optical flat. The optical flat of fused 
silica or other glass must have both surfaces 
polished “true plane’? and adjusted to make an 
angle of 20’ with each other so that interference 
bands will not be formed by reflection between the 
upper and lower surfaces. 

The specimens of glass for chemical durability 
measurements were small rectangular strips ap- 
proximately 3 em long with one surface piteh- 
polished, and were sufficiently flat to show inter- 
ference bands when placed under an optical flat. 
The exposures were made by immersing these 
specimens to approximately one-third of their 
length in the desired solution. The surface of the 
solution in most instances was covered throughout 
the test period by a thin layer of liquid petrolatum 


to prevent evaporation. For the experiments 





? An attack of less than 0.1 of an interference fringe can be readily detected 
For a glass of density 2.5, an attack of 0.1 fringe (filtered helium source 
amounts to approximately: 


S87 ™ 10 ' ie = => 
ox = xix? 0.073 10—*g or 0.0073 me/em? 


with high percentages of ethanol (C,H,OH), +}, 
density of the solutions was too low to float th 
oil layer. In these experiments it was necess ry 
to make the exposures in a closed system by sus 
pending the specimen from the cork stopper of a 
Erlenmeyer flask.* 

The desired temperatures were maintained con 
stant to +0.2° C by means of an electricall 
heated thermostated water bath. All chemica’ 
durability data were obtained at 80° C except fo 
the ethanol and hydrofluoric acid solutions. Fo 
the ethanol solutions the temperature was lowere: 
to 78° C to insure that the surface remained wn. 
disturbed by boiling during the period of exposur: 
In the hydrofluoric acid solutions it was necessary 
to make the durability exposures at room temper. 
ature because of the vigorous attack on the glass 

As the purpose of the investigation was to fu 
nish a comparison between the voltage departur 
of the glass electrode and the chemical durability 
of the glass, the solutions were adjusted to approx. 
imate concentrations only. Whenever possible 
the voltage departure of the glass electrode was 
observed for a given solution, followed immediate 
by a determination of the durability of the glass 
in the same solution, 


Ill. Voltage Departures (Errors) of the 
Glass Electrode and the Chemical Du- 
rability of the Glass 


1. Sodium Hydroxide Solutions 


One of the most conspicuous voltage departures 
for the glass electrode occurs in the alkaline range 
beginning near pH 9 and increasing rapidly with 
increased alkalinity [1]. It has been shown that 
this region of voltage departure is accompanied 
by changing chemical durability of the electrode 
glass, the rate of attack increasing with increasing 
pH |4]. Published data covering a limited range 
of conditions indicate that this relation is linear 
[S, 9]. 
make as simple an experiment as possible, an 
aqueous solution of 1.0 N NaOH was prepared 
from 
ranging in concentration from 0.015 N to 1.0 N, 
with each successive member of the series doubling 


To check this observation further and to 


which a series of solutions was obtained 


' The use of a closed system v is tried at the suggestion of A. W. Helv. of 
the Portland Cement Association Fellowship. It proved very effective 
greatly increases the range of usefulness of the interferometer method 
durability measurements. 
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with the 


fferent solutions 


glass electrode chemical durability of the glass in 


NH,OH, Ba(OH 


hours at SO° ( 


NaOH + Nat ( 


Dural t\ \posure 6 


‘20H solutions; B KOH, LIOH 


ind NaOH 


lation between voltage departure and attack 
The 


aken in order beginning with the weakest solu- 


vithin reasonable limits emf data were 


on. The individual readings were usually steady ; 
iwowever, In accordance with published data, the 
eproducibility was found to depend greatly on 


|8, 9, 12). Nevertheless, the data are sufficiently 
reliable to support the trends indicated by the 


curves, 


2. Sodium Hydroxide Solutions Plus Sodium 
Chloride 


It is known that the voltage departures of the 
glass electrode in alkaline solutions are increased 
by the addition of neutral sodium salts such as 
NaCl [1 p. 129, 14, 15, 16]. Chemical durability 
data have shown qualitatively that an increase in 
attack on the glass also accompanied this increased 
voltage departure [13]. To obtain quantitative 
data on this point, the same NaOH series as 
deseribed in the prev ious section was prepared and 
sufficient NaCl added to make all of these solutions 
Once again the 
individual emf readings were steady. The results 
B of table 1 and 

A pronounced 


normal with respect to Na* ions 


obtained are listed in section 

plotted as curve B of figure 1 
increase in both voltage departure and attack was 
brought about by the increased! |Na*]. However, 
the voltage departure exhibited a greater increase 
than did the attack, and the data do not coincide 
with the previous straight line A. From this it is 
evident that factors other than the durability shift 
of the glass contribute to the voltage departure of 


the electrode. 


3. Solutions Containing Na*, Li*, K*, Ba**, and 
NH,* ions 

The magnitude of the voltage departure in alka- 

line solutions is different for different ions [1, p. 

130]; for example, the departures caused by cations 

the Nat >Li* >K*>Bat*>NH,?*. 


are in order, 


he immediate previous history of the electrode Karlier work suggested that this was the order 
BLE I Comparison of the voltage departures (errors) of the glass electrode with the chemical durability of the glass n solu- 
tion of NaOH, NaOH plus NaCl, and other hydroxides 
Durability exposure 6 hours at 80° C 
\ B Cc 
Solution Depar Attack Solution Depar- Attack (fringes Solutio Depar Attack (fringes) 
ture ure ture 
ur uw 
e HC] 0 Zio (swelling Dilute HC} 0 o (swell Dilute HC] 0 249 (swelling) 
N NaOH 3 attack L015 N NaOH+Nal nv 1 (attack N NH,OH attack 
N i is | O31 N do 3 5 do Ba(OH), satu 20 Loy... che 
N i 61 i O62 N i 107 pi, do rated 
N do S 2 i 125 N do 122 3 ck N KOH 9 24 do 
N do 105 | 284 lo 25 N | 135 sche N LiOH 78 | 2 lo 
N do 133 tly i WN i 51 do N NaOH 140 | 
‘N do 140 Sty i 
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of decreasing attack on the electrode glass by 
these ions [4, 13]. To check this point, normal 
solutions were prepared for each of these hydrox- 
ides except Ba(OH),. For Ba(OH), a saturated 
solution was used. The voltage departure of the 
glass electrode was determined in each of these 
solutions, followed by a durability test of the 
glass. As the concentration of saturated Ba(OH), 
at room temperature is less than 1 N, and at 
80° C is greater than 1 N, the voltage departure 
value is probably too low and _ the 
attack value is too high. The results given in 
section C of table 1 and plotted as curve C in 
figure | show that the relative attack on the glass 


obtained 


by these solutions is in the same order as for, the 
voltage departure, but the two properties do not 
plot as straight line nor does the resulting curve C 
coincide with either curve A or curve B. These 
results again indicate that some factor other than 
durability shift of the glass contributes to the 
errors of the glass electrode. However, before one 
accepts these data as an indication of a specific 
equilibrative response of the glass electrode to 
these ions, it must be emphasized that the glass 
electrode gives little or no specific response to any 
of them except 
durability of the glass. 


under conditions of changing 


4. Buffer Solutions Containing NH,* ions 


The voltage departures of the glass electrodes 
are much less in alkaline buffer solutions contain- 
ing NH,* ions than in basie solutions prepared 
with Na* ions [17]. A typical set of data illustrat- 
ing this point and making a comparison of voltage 
departure and chemical durability is given in 
table 2 and plotted in figure 2. The voltage 
obtained on buffer 
prepared from the Britton-Robinson 
buffer mixture [18] in which ammonia was sub- 
stituted for the NaOH to obtain the desired pH 
range, the highest pH value being that of full 


departures were solutions 


universal 


strength ammonia solution.’ 

Most of the voltage results are qualitatively 
similar to those for buffers containing Na* ions, 
the departure remaining nearly constant in the 
acid range, then increasing with increasing alka- 
However, at 


linity beyond pli s. the higher 


alkalinities beyond pH 11, the rate of departure 


————_————$<$— 


‘In the durability tests no difficulty was encountered from loss of NH 


through the oil layer except at pH 12 und for the full strength solution 
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Figure 2. 
glass electrode with the chemical durability of the glass i 


Comparison of the voltage de partures o 


buffe r solutions containing N H,* ions. 


In this and other figures the s 
the circles does not indicate the probable error of the individual px 


Durability expsoure 6 hours at 80° C 


for durability, but merely emphasice that the data must not be consice 
highly precise nor accurate 


TABLE 2. vollage 


of the glass electrode 


Comparison of the de partures (errors 


with the chemical durability of U 
glass in buffer solutions containing NUH,* ions 


[Durability exposure 6 hours at 80° ¢ 


Surface cut 


pi oe Attack (fringes oil-solutic 
interface 

2.9 0 eo Swelling None 
a) 0 219 Swelling Do 
4.9 0 Detectabk 
a 2 i 
9.7 1 bie Attack Lio 
10.3 7 ‘io Attack ! 
11.0 s » Attack ‘ 
12.3 7) Attack None 
12.9 full strength NH S5to0 Detectabl Slight 


NH, escaped through the oil layer during the early minutes of the att 
period, 

? The temperature of the full-strength NH, fell sharply upon bubbh 
hydrogen through the solution 


falls off and reversal for full 
strength NH, even yielding a negative 
the latter This value for full 
strength NH, bas been previously reported [4] bu! 
not adequately The  hydroge 
stream bubbling through the full strength an 


goes through a 
value in 


case. negative 


investigated. 
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«nia causes a large drop in the temperature of 
i solution. When adequate precautions were 
aken to minimize this temperature change, the 
oltage departure for the glass electrode approxi- 
wited zero. 

The trend of the departure curve over the pH 
uige investigated was in complete accord with the 
irability curve for the glass in these same buffers. 
n the acid range*the glass exhibited a slight 

welling, and between pH 6.9 and pH 11.0 exces- 
ive surface cuts appeared at the oil-solution 
nterface, similar to the results reported for buffer 
«lutions containing Na* ions [4]. An increase in 
Hat the oil-solution interface brought stout by 
referential orientation of ions could account for 
he surface cuts in this pH range. In full- 
strength NH, the durability 
roved in accordance with that indicated by the 


of the glass im- 
oltage departure shift. 


5. Electrode Prepared from a Glass of Low 
Hygroscopicity 


The voltage departures for Corning 015 glass 
lectrodes that appear under conditions of chang- 
ng chemical durability of the glass [4] must not 
« confused with the voltage errors that are ex- 
ibited by electrodes prepared from glasses of low 
vgroscopicity [1, 2,3,]. In the latter, departures 
ppear over the entire pH range, with any depar- 
ires accompanying a change in durability being 
wrely superimposed on the general departure 
irve [3]. A typical example of the behavior of 
ectrodes prepared from glasses of low and high 
iwgroscopicity, Pyrex and Corning 015, respec- 
vely, is illustrated in table 3 and figure 3. These 
lata were obtained in Britton-Robinson universal 


BLE 3 Comparison of the voltage departure of an elec 
ule prepared from a glass of low hygroscopicity, Pyrex, 
departure of an electrode Corning O14 glass, 


} yd ogen elect ode as clerence electrode 


Voltage Anomalies of the Glass Electrode 


Qcevporture, millivolts 
8 
T 
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Figure 3. 


trode prepared from a glass of low hygroscopicity (Pyrea 


Comparison of the vollage de parture of an elec- 


with the de parture of an electrode of Corning 016 glass. 


hese data were taken over a pH range for which both glasses exhibit un 


form durability 


buffers over a pH range in which both glasses 
exhibit uniform durability [13], in order to elim- 
inate any voltage effects attributable to change in 
the durability of either glass. 


6. Hydrofluoric Acid Solutions 


Electrodes of Corning 015 glass exhibit very 


pronounced voltage departures in hydrofluoric 
acid (HF) solutions, the departures increasing 
with increased concentration [4, 16]. The data 
given in table 4 and plotted in figure 4 show that 
these increased voltage departures are accom- 
panied by increased attack on the glass. How- 
ever, the relation between voltage departure and 
Although the attack 


on the glass was very severe, the emf readings 


attack is not a straight line. 


were surprisingly steady. The attack was so 
rapid that to keep the durability measurements 
within the range of the interferometer, the time 
of exposure had to be reduced to 5 minutes at 
room further shortened to 1 
minute, 30 seconds, and 10 seconds for the 2.0 N, 
3.0 N and 5.0 N solutions, respectively. Hence, 
the results shown in figure 4 are only qualitative, 


temperature and 


but oby iously the reneral shape of the curve will 
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Fieure 4. 


Comparison of the voltage departures of the glass electrode with the chemical durability of the glass in 


aqueous solutions of hydrofluoric acid. 


Durability exposure 5 minutes at room temperature 


not be altered even though the individual values 
may possess large cumulative errors. 


TaBLe 4. 
of the 


Corning 015 glass in aqueous solutions of hydrofluoric 


Comparison of the voltage departures (errors) 


glass electrode with the 


acid 


{Attack reported for 5 minutes exposure at room temperature] 


Solution Depat Attack (fringes 
ture 
me 
Dilute HCI 0 Not detectable. 
oooOLN HF 3 Do 
orN HF 12. Detectable 
WN HE 71 by 
SN HF 1m) 2 
LON HE Is! tha 
20 N HF a4 11h 
LON HF 21s 21 
SON HF —208 145 


' Values based on attacks observed for | minute, 30 seconds, and 10 seconds 
for concentrations 2 N, 3. N, and 5 N, respectively. 


7. Ethanol and Acid Solutions 


The voltage departures of the glass electrode 
that appear in nonaqueous such as 
ethanol (C,H,OH) are often compared with the 
characteristic departures in strong acid solutions 
(1, 5, 19, 20, 21]. The increasing negative de- 


solutions 


parture of the glass electrode in acid solutions 
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chemical durability of 


of increasing concentration is accompanied by a 
repression of swelling of the glass [4]. Similar 
durability experiments carried out for ethano! 
solutions indicated that a repression of swelling 
with increased ehtanol concentration also accom- 
the negative voltage departure. How- 
been convineing 


panied 
ever, the latter data had not 
because of the difficulty of applying the inter- 
ferometer procedure for durability, as the low 
density of the ethanol solutions prevented the 
floating of the oil layer. To obtain more con- 
clusive data demonstrating that a durability 
shift of the glass does accompany the voltage 
departure in ethanol solutions, durability ex- 
posures were made in a closed system that elimi- 
nates the need for a layer of oil. 

Table 5 and figure 5 give the results obtaine:! 
for the repression of swelling of Corning 015 glass 
and the accompanying voltage departures of the 
glass electrode in water-ethanol solutions buffered 
at pH 4 with Britton-Robinson universe! bufter 
mixtures [18]. The performance of this glass in 
ethanol appears to be similar to its performance 
in strong acid and salt solutions. This can lx 
seen from inspection of table 5 and figures 6, 7, 8 
and 9, which give similar data for sulfuric [5| 
formic, and acetic acids [21] and magnesium 
sulfate [22]. All of these have been plotted on a 
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percentage basis in order to make an easy inter- 
comparison and also to avoid the confusion of 
thought caused by the abnormally large activities 
for high concentrations of strong acids |23}. 


TaBLe 5 Comparison ol the voltage de parture o the glass 
electrode with the chemical durability of Corning 015 glass 
in solutions of C,H;,OH, H,SO,, HCOOH, CH,COOH 
and MgSO, 


Solutior Departure Attack (swelling 
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Fictre 9 Comparison of the voltage departure of the gla 
electrode with the chemical durability of the glass repres- 
sion of swelling) in acid solutions of increased MgSO 
content, 


Durability exposure 48 hour. at 80° C 


The repression of swelling by the higher concen- 
trations of organic and inorganic acids, as well as 
for the magnesium sulphate ° is readily predictabl 
from Donnan membrane considerations and_ thy 
distribution law [24]. The accompanying voltag: 
departure can be attributed to the chemical chang: 
of the exposed glass surface [4], to a reduction of 
the water activity of the solution [1, 19, 20], or to 
a loss of electrode function because of induced 
decrease in water content of the glass [1, 2, 3 
Dehydration of the glass appears to be the most 
likely cause, as the major voltage departure is 
taking place over a concentration range in which 
no further detectable change in durability occurs 

For alkaline solutions of ethanol, the data 
obtained were far from satisfactory or conclusive, 
table 6 and [1 p. 81, 25]. Solutions of increased 
C,H,OH content, when prepared immediately 
before use, gave voltage departures that plotted 
as a straight line for 0.1 N and 0.01 N NaOH 
curves A and B, figure 10. However, afier the 
attack data, curves C and D, were obtained, thy 
voltage-departure curves E and F for these samy 
solutions were entirely different from curves A 
and B. The response of the hydrogen-gla 
electrode cell was sluggish, giving the impression 


A drop of H»SO, was added to the MgSO, solutions to insure t 


nain in the acid pH rang 
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that the hydrogen electrode was poisoned. How- 
ever, checking the combination in the dilute HCI 
solution used to establish the zero reference poten- 
tial showed that the initial voltage was duplicated 
to within 2 mv. This indicated that the electrodes 
had not been affected permanently and that some 
alteration had taken place in the solutions during 
the 6-hour period at the high temperature. Fur- 
ther data taken on 0.1 N and 0.01 N NaCH solu- 
tions prepared from a stock solution of NaOH in 
C,H;,OH gave the departure curves G and H, 
respectively. The latter curves might be rational- 
ized as compatible with the durability curves C 
and D, but it seems unlikely that any of the emf 
values are reliable. At least, it seems improbable 
that all of the voltage anomalies observed can be 
accredited to the glass electrode. However, it 
must be borne in mind that the glass is being 
vigorously attacked by some of these solutions, so 
that reliable performance of the glass electrode is 
not to be expected. 


IV. Conclusions 


Voltage alkaline solutions of 


electrodes prepared from Corning 015 glass are 


departures in 


accompanied by changes in the chemical dura- 
bility of the glass. The positions in a voltage 
departure series occupied by different metal ions 
are in the same order as the attack on the electrode 
glass by these same ions. However, the previous 


suggestion that the magnitude of all voltage 
departures might be directly proportional to the 
extent of the change in attack on the electrode 
glass is not universally true; in fact, the data 
obtained show that such linear relations occur for 
only a limited number of conditions. Further- 


more, in acid solutions of high concentration 


with the exception of HF), ethanol, and also 
salts of high ionie strength; dehydration of the 
glass is probably responsible for voltage depar 
tures, 
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Instability of Simply Supported Square Plate With 
™ Reinforced Circular Hole in Edge Compression 


mt to 


“4 By Samuel Levy, Ruth M. Woolley, and Wilhelmina D. Kroll 


Bool 
l 
" A method is presented for computing the compressive buckling load of a simply supported 
elastic rectangular plate having a central circular hole reinforced by a circular doubler 
plate. Numerical results are presented for six square plates having hole diameters up to 
one-half of the plate length. Comparison of these results with those computed for plates 
without holes shows that an unreinforced circular hole causes a relatively small reduction in 
buckling load, and reinforcement of a circular hole by a doubler plate causes a substantial 
increase in buckling load 
I. Introduction around a circular hole in an infinite sheet under 
- uniform compression in one direction or under 
The stressed skin cover of airplane wings and 
shear. 
iselages has to be perforated by holes at a number on ica Sn =e 
‘ The second phase of this investigation consisted 
if places to give access to the interior of the wing ; 
fe of a check of the theoretical analysis by tests on 
r fuselage. These holes must be reinforced to . . . 
: ' ; plates having reinforced circular holes. It was 
wevent weakening of the entire structure by the é‘ 
found that the analysis gave an adequate descrip- 
Ore, ° ° ° 
, ' tion of the stresses and displacements provided 
{n ideal remforcement would be one that . 
. . . that the doubler plate was fastened to the sheet 
onfines the disturbance in the stress flow to the 
; , , by at least two concentric rows of rivets. 
mmediate neighborhood of the hole and which, at at “ake nt nee , 
, , The third phase of this investigation is described 
the same time, adds minimum weight to the ; air ange 0 
in this paper. An analysis is given for the stabiliz- 
structure, : . ; 
rh . . . ing effect of the doubler plate for square plates 
Too little is known about reinforcements to 
— with central circular holes when the plates are 
ipproximate this ideal in practice. In the 2 Cite 
subjected to edge compression in one direction 
ibsence of a rational design procedure, it has been 
ustomary to reinforce holes by circular doubler , 
a - a II. Method of Analysis 
plates riveted to one side of the sheet, thereby 
nereasing the effective thickness and rigidity of An energy method for determining the buckling 
the sheet at the edge of the hole. load of rectangular plates of constant thickness 
The Bureau of Aeronautics, Navy Department, under compressive loads is presented by Timo- 
nitiated a study of reinforcements around holes shenko.' A review of this derivation shows that 
it the National Bureau of Standards to provide a the method is also applicable to plates of variable 
etter understanding of reinforcements and to thickness. In this case, Timoshenko’s integrals 
ndicate practical improvements in their design. I, and J, are 
lhe first phase of this investigation was a plane PEE 
a ‘ S. Timoshenko, Theory of elastic stability pp. 307 to 326 (McGraw-Hill 
stress analysis of a doubler plate reinforcement Book Co., New York, N. Y.. 193¢ 
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where 
h=plate thickness (function of + and y). 
r, y=rectangular coordinates with origin at 
center of plate and z-axis in direction 
of load. 
lateral deflection of plate. 
D= Eh 12(1—p*), flexural rigidity of plate 
(function of ¢ and y). 


~ 


u—0.3, Poisson’s ratio. 

tensile stress in x direction. 
, = tensile stress in y direction. 
Ty ~shear stress. 


S= tensile stress in « direction far from hole. 


The stress ratios o,/S, o,/S, and 7,,/S just prior 
to buckling may be obtained as shown by Gur- 
ney. The stresses o,, o, and +,» referred to 
polar coordinates r, @ are given, with some change 
in notation, as 


o,—S' F+ KR?*/r?— (A+ 3CR'/ r+ 2DR?/r*) cos 20} 
og S| F—KR?/r? + (A+ 6Br*/ R? + 3CR'/r') cos 20 
t~a=S{ (A+3Br/R?—3CR'/r' DR: r?) sin 26) 


(2) 


where 


r, @=polar coordinates with origin at center of 
hole, 
@=angle between radius r 
tion of load, 
R radius of hole. 


point and diree- 


Fk, Kk, A, B,C, D=coeflicients with different 
values in sheet and in 
doubler plate. 

The values of F, A, A, B, C, and D may be deter- 
mined from the dimensions for a_ particular 
doubler plate by using figure 1. The coefficients 
from “‘1"’ region, figure 1, 
computing the stresses in the unreinforced sheet 
the coefficients from ‘‘2” 


are used in eq 2 when 


region are correspond- 


G. Gurney, Brit. Rep. Memo. No. 1834 (Feb. 1938). 
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ingly used when computing stresses in the r. iy 
forced sheet. The abscissa in figure 1, (7 
(6?/a?-1), is the ratio of the volume of materia) | 
the reinforcement to the volume of materis 
removed from the sheet to make the hole. 

The stresses in rectangular coordinates may |) 
computed from the stresses in polar coordinat, 


as given in eq 2 by using the conversion formula. 


74 sin 26 

Te Sin 20 
Co, 0a ° 

Try ( = ) sin 20+ 7.9 cos 26. 


o,= 6, cos*6+ a9 sin’6- 
o, =, SIN*O+- a9 COs*O4 


The lateral deflection w of a simply supporte: 
rectangular plate of length a and width 6 will } 
approximated by the first terms in the trigono- 
metric series 


wr ry wr Sry 
cos | (y, COS cos . 
a b Fis a b + 


wea), COs 


r ry Sarr Sry 
COS of +d COs cos ;°+ .. (4 
b a bh 


3m 
ay, COs 
a 


Substitution of eq 4 in eq 1 and integratio: 
(see next section for method of numerical in 
tegration) reduces integrals 7; and /, to quadraty 
expressions in the coefficients a@)), ays 
The critical value of S at which buckling of th, 
plate occurs is obtained from these expressions 
according to Timoshenko, as that value of 4 
which reduces to zero the determinant of th: 
coefficients of a,;, a, ete. in the set of simultaneous 
equations, 


ol, , of, 
‘ 0 

Oa, On); 

ol, _ ol, 5 
S =. ==() 

Od); Ol) ; 





The sign before S in eq 5 is plus, and that before 4 
in the corresponding eq 212 given by Timoshenko 
Where S is considered positive fot 
tensile stress, as is done in this report, the signs 1) 


is minus. 


eq 5 are correct. 
iS. Timoshenko, Theory of elastic stability (McGraw-Hill Book | 
New York, N. Y., 1936 
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re 4 FiGure 1 Values of coefficients F, K, A, B, C, and D to be used in eq 2 for computing the stress distribution in 
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Ill. Numerical Integration 


The evaluation of the integrals in eq 1 over the 
surface of the plate is made difficult by having a 
circular inner and a rectangular outer boundary 
and by involving a stress that is a complicated 
function when expressed in rectangular coordi- 
hates. 

The integral /,;, eq 1, was evaluated as the 
resultant of three integrations: 


L=lfhetlp—l, (6) 
where 

I,, ~The integral J, for the surface enclosed 
by the outer rectangular boundary of the 
plate taking Y—D),, the value of DP in 

the unreinforced portion of the plate. 
/,,—The integral /, for the surface enclosed 
by the outer circular boundary of the 
D),, where 
D, is the value of ) in the reinforcement. 


reinforcement taking J)= J), 


/ The integral /, for the surface enclosed 
by the cireular boundary of the hole 


taking J)— 1),. 


The double integration necessary to evaluate /,, 
could in every case be done directly. However, 
this was not possible in evaluating /,, and /,.. 
The integrals /,, and /,, were obtained by integrat- 
ing directly in respect to one variable and using 
Gauss’ method of numerical integration,’ for in- 
tegration with respect to the other variable. In 
each case where numerical integration was used to 
evaluate /,, the number of Gauss points was in- 
creased to the place where the addition of another 
point caused less than | percent change in the 
integral. In general, five Gauss points were 
enough for this purpose. 

The integral 7, was evaluated by considering 
the surface of the plate in three portions. The 
first portion, A, figure 2, was taken as the cireular 
disk of the reinforced area; the second portion, 
B, figure 2, was taken as the circular disk between 
the outer edge of the reinforcement and the largest 
inseribed circle in the plate; and the third portion, 
(’, figure 2, was taken as the remainder of the 
plate. The integral for each cireular portion was 
determined, using Gauss’ method of numerical 
integration, by first integrating numerically in a 
— 


‘FE. T. Whittaker and G. Robinson, Calculus of observations, p 
(Blackie and Sen, Ltd., London, 144 
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Figure 2.——Subdivision of plate into three portions A, B, C, 


in evaluating I>. 


circumferential direction and then integrating 
re. : 

The integral 

for the remainer of the plate was obtained also by 

using Gauss’ method, first integrating numerically 


numerically in a radial direction. 


in the direction of the load and then at right angles 
to the direction of the load. Only one quadrant 
of the plate had to be considered because of sym- 
metry. In each ease, a sufficient number of 
Gauss points was used to reduce the estimated 
error to less than 1 percent. Twenty-three points 
in each quadrant were used for cases 1, 3, 4, and 
5, table 1, thirty-three points for case 2, and 
twenty-eight points for case 6. 


TABLE 1.-- Dimensions for plates investigated 
Width rhickness 
Outer 
Cast Radius of radius of 

. Direc- rrans- hole, R reinforce- | Rein- | Unrein- 
tion of verse to ment forced forced 
load loud region | region 

l a a 0. 1254 ! ( h 

2 a a OO25a 1 i h 

a a 2ha ! h 

4 a a 125a 0. 25a 2. 0h h 

a a 125a 25a 1. 5h h 

6 1 a 125a . IS75a 3. 4h h 


' Hole in plate not reinforced 


A typical distribution of Gauss points for 
computing /, is shown in figure 3 for one quadrant 
of a square plate with a central hole. For this 
distribution, the integral of a function F' over the 
plate surface is 
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FiGuRE 3 Distribution of Gauss points for computing |. in case 
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«(0.0753 16F, + 0.068366 F, + 


SS FAA 


surface 

0.068366F),+ 0.041605F,-4 
0.084048 F', + 0.099898 F;, + 
0.084048 F,+- 0.041605F 
0.028176F,+ 0.0569 19F; 
0.067653 F;, + 0.0569 19F,, 4 
0.028176F),+0.01040LF;, + 
0.021012 F,, + 0.024975F;,4 
0.021012F,,+- 0.010401 F\, 4 
0.007044 F),, + 0.014230F,, 4 
0.016913 F,, + 0.014230F,, 4 
0.007044 F | (7) 

where F,, Fy... Fo are the values of the 


function F at the positions indicated in figure 3. 

An approximate check on the over-all adequacy 
of the numerical integration methods used was 
obtained by determining the buckling load of a 
square plate with no hole both by the exact 
method,’ and by numerical integration with 23 
points for /,. The differed by only 
0.7 pereent. In the case of a plate with a rein- 
forced hole, it is likley that the more complicated 
stress distribution and the greater prominence 


results 


of higher order terms in the series used for the 
deflection causes the error to be somewhat higher. 


IV. Convergence of Deflection Function 


The correct value of S for buckling of the plate 
would be that value of S which reduces the deter- 
minant of the coefficients of ay,, ay.. in the 
infinite set of eq 5 to zero. In order to limit the 
work of computing Sto a finite amount, preliminary 
computations were made to see which terms in 
the deflection function, eq 4, were most important 
in determining the value of S for buckling. These 
computations were made for the square plate of 
ease 4, isble 4b. 
10.33 Eh?/a? using only the first term of eq 4, 


The compressive stress (—S) was 


wr ry 
we ay, COS Cos 
a 4 
Using only terms | and 2 of eq 4, 


ws Ony 
+ Ay, COS cos me 
a b 


wr ry 

w=), COS cos 7 

a b 

(—S) was 9.48 Eh? a’, 
Using only terms | and 3, 


a decrease of 8.6 percent. 


$8. Timoshenko, Theory of elastic stability, p. 331 (MeGraw-Hill Book Co., 
New York, N. Y., 1936 
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ws ry 3rr ry 
w=), COS COs > dg, COS cos "> 
a b a b 
(—S) was 8.80 Eh? /a’, 
Using only terms | and 4, 


a decrease of 14.7 pereeit 


I ory 
w=), COS cos -% 


wr ry on 
q COS pf tas C08 
(—S) was 10.32 Eh? /a?, 
Had additional terms of the trigonometric series 
been used with term 1, the decrease in the buckling 
stress would probably have been proportionate!) 


a decrease of 0.1 percent 


smaller. It is believed, on this basis, that thy 
first three terms of eq 4 approximated the latera 
deflection of the square plates with simply sup- 
ported edges with sufficient accuracy to give thy 
buckling stresses within 5 pereent. Accordingly 
only the terms 


ry wr Sry 


mr 
wd), COs cos +), COS > COS + 
" a b ' a b 


wi] 


Srsr 
(ls, COS cos 
a b 


were used in the remaining computations for 


square plates. 
V. Results 


The compressive buckling load was computed 
for six square plates with reinforced and unre- 
inforced holes, figure 4. The dimensions for thi 
plates are given in table 1. 

The analysis gave the compressive buckling 
stress (—S) far from the hole corresponding to 
the stress distribution derived in reference {2} 
Values of ( 


pressive stress ( 


S) are given in table 2. The com 
S) is somewhat larger than the 
average stress in the direction of the load along 
A better value of the 


average buckling stress o,, along the edge of the 


the edge of the plate. 


plate was obtained from 


where o, is the average compressive stress on the 
loaded edge of the plate and o, the average com 
pressive stress on the center cross section, obtained 
by dividing the load on that section by ah. The 
value of o; was obtained by numerical integration 
and the value of o by direct integration of the 
Values of (—S), o;, 
The average stresses 


stresses o, given by eq 3. 


and o,, are given in table 2. 
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Figure 4 Loading of square plates with simply supported 


edges. 


Dimensions in table | 


and o, differed from each other by less than 3 
percent, indicating that the equilibrium condition 
is substantially satisfied. 


PARLE 2 Critical compressive stress for plates investigated 
Average tre Average stress 
0 © on at @) correspond- Buckling stress 
Case mn _ ng to loud at | oer, average of 
or role a a a ‘ " t 
far from hob Se eg 0 divided og, ande 
é é é by ah 
! SOL Eh a 8.454 Eh $400 Ehta 8472 heya 
2 8.620 Eh?a 8.578 Eh? a 501 Pha t.5M4 Eeh2a 
8.75) Eh? a? L076 Fh 8.164 Fh a2 $120 Eh? at 
i 8.021 Ehta 7.084 Eh? a S000 Fh? at 8.027 Ehta 
5.700 Eh2/a2 S578 Rhea 5146 Eha LHl2 Rhea 
6 7.945 Eh? a 7.SSI Eh 7.84 Phat 7.932 Eh 
The buckling stress of a simply supported 


square plate of constant thickness without a hole 
is given on p. 331 of reference [3] (taking Poisson's 
ratio as 0.3) as 


o---3.615 Eh?/a’. (9) 
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Comparison of this value with the values computed 
for cases 1 to 3 of table 2 shows that unreinforced 
central holes, having a diameter of \, 4 and ' of 
the length of a side of a square plate, reduce its 
buckling load by only 1, 4, and 14 percent, 
respectively. Comparison of the value in eq 9 
with the values for cases 4 to 6 of table 2 shows 
that the reinforcements increased the buckling 
stress over that for the plate without a hole by 55 
to 122 percent. 

The effect of the shape of the reinforcement. on 
the buckling load is indicated by comparing cases 
4 and 6, table 2. 
volume of reinforcing material are identical but 
the reinforcing material is concentrated nearer the 
The com- 


For these cases, the hole size and 


edge of the hole in case 6 (see fig. 4). 
puted buckling stresses are the same within 1.2 
pereent, 


VI. Conclusions 


A numerical procedure was developed for 
estimating the buckling stress of simply supported 
rectangular plates with circular holes and doubler 
plate reinforcement. The procedure provides a 
convenient method for solving the integrals for 
the energy stored in the plate. 

The computations showed that the buckling 
stress of square plates is reduced only a small 
amount by the presence of unreinforced holes. For 
the cases considered, the greatest reduction, 14 
percent, corresponded to case 3, asquare plate with 
a hole diameter 0.5 times the length of the plate. 

Reinforcements of the doubler-plate type are 
the computations to cause marked 


For case 4 with 


shown by 
increases in the buckling stress. 
reinforcement, the computed buckling stress was 
2.3 times that for case 1 with no reinforcement. 
The buckling stress seems to be insensitive to 
shape of 
-_ 
“he computed buckling stresses were 


change in reinforcements of given 
volume. 
the same within 1.2 percent for cases 4 and 6, 
which differed in having doubler plates of different 


thickness and radius, but of the same volume. 


The authors thank the Bureau of Aeronautics, 
Navy Department, for their 
carrying out this investigation and for permission 
to publish the results. 
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